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Co3aaH B WeJIAxX MHPopMupoBaHHA YHTaTelIbCKOH ayHHTOpHu oO HOBeHIIMX AOCTWOKCHHAX HM HepciieKTHBaXx B 
oOmacTH MeCXaHHKH, MallMHOCTpOecHHA, MHPOpMaTuKu VM BBIYMCIUTeIbHOM TeXHUKH. V3 manne ABJIACTCA cbopyMoM 
ASIA COTPYTHHAeCTBAa poccniickux MW HMHOCTpaHHbIxX YUCHBHIX, cioco6cTByeT COUIM2KeHHIO pocculickoro MW MHpoBoro 
Hay4HO-HHOpMallHOHHOrO TIPOCTpaHcTBa. 


2KypHal BKOUeH B Tepeyenb pelleH3sipyeMbIX Hay4HbIX 34aHHi, B KOTOPOM OJDKHBbI ObITb 
onyOJIMKOBaHbI OCHOBHbIe Hay4Hble pe3yIbTaTbI AMcceprayHi Ha concKaHHe y4ueHoH cTeneHH KaHAMAaTa 
HayK, Ha concKaHHe y4eHol creneHH fOKTOpa HayK (Ilepeyent BAK) no cueqyromuM Hay4dHbiM 
cienMaJIbHOCTAM: 


1.1.7 — Teoperuueckas MexaHuka, JMHaMMkKa MallMH (TexHHueckHe HaykKH) 

1.1.8 — Mexanuxka tecbopmupyemoro TBepyoro Tesa (TexXHHYecKHe, (PH3HKO-MaTeMaTHYecKHe HayKi) 

1.1.9 — Mexannka 2KHJIKOCTH, ra3a HW TW1a3MbI (TeXHH4eCKHe HayKH) 

1.2.2 — Matematw4eckoe MoyesIMpoBaHie, YMCJICHHbIC METOJbI HW KOMIUICKCbI IporpaMM (TexHv4eckHe HayKH) 

2.3.1 — CucTemHbIii anamu3, ypaBienve 41 oOpaboTka HHopMallun, cTaTucTuka (TexHHueckHe HaykKH) 

2.3.3 — ABTOMaTH3al[HA HM yIpaBleHHe TEXHOJIOTHYeCKHMH MpoleccaMy H TpOu3BOACTBaMH (TeXHH4eCcKHe HayK) 

2.3.5 — Maremarvueckoe 4 TIporpaMMHoe oOecreveHHe BLMHCJIMTCJIBHBIX CHCTeM, KOMIJICKCOB HM KOMIIBIOTepHBIX CeTeli (TeXHHyecKHe 
HayKH) 

2.3.7 — KommpriotepHoe MojesIMpoBaHle H aBTOMAaTH3alHA MIpOeKTHPOBaHHA (TexHHYecKHe, (PH3HKO-MaTeMaTHYecKHe HaykKH) 
2.3.8 — Mndopmatuka 4 HHC*OpMallMOHHEIe Mpoleccs! (TexHHueckHe HaykKH) 

2.5.2 — Mammmyopeyenne (TexHu4ecKHe HayKi) 

2.5.3 — Tpenue 1 W3HOC B MalliMHax (TexHM4ecKHe HayKH) 

2.5.5 — Texuosorna 1 ob6opyqoBaHve MexaHHueckoli M PU3MKO-TexHHYecKoH OOpaboTKH (TexHH4eCKHe HaykKH) 

2.5.6 — TexHosIOrMa MalMHOCTpoeHHaA (TeXHHYeCKHe HaykKM) 

2.5.8 — CBapka, poCTBeHHEIe IIpoleccbi H TexHOJIOrHH (TexHM4eCKHe HayKH) 

2.5.9 — Metogsr 4 mpHOopbl KOHTpOIA MH WHarHOCTHKM MaTepHasIOB, H3]eIM, BeLeCTB H NpHpoqHO cpegbI (TexHH4ecKkHe 
HayKu) 

2.5.10 — [uapapmyeckne MallHubl, BAKyyMHas, KOMMpeccopHasd TeXHHKa, THAPO- Mf THeEBMOCHCTeMBI (TeXHHYeCKHe HayKM) 
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Capyxansn Apecrak Apamancosn, JOKTOp TexHHYecKHX HayK, Tpodeccop, HalvouasibHbit yHHBepCHTeT apXHTeKTypbI H CTpoHTesiEctBa ApMerun (ApMeHHa); 
Cuyopos Baagumup Hukosaesuy, JoKTop TexHuueckux HayK, Poccuiickuii yHuBepcutet Tpaucnopta (Mocxsa, Poccuiicxaa Deyepaus); 
Conopbés Apkaguii Hukosaepny, JoKTOp Pv3HKO-MaTeMaTHYeCcKHX Hayk, Upodeccop, JloucKoli rocyqapcTBeHHBIt TexHHueCKHM YHHBepCHTeT 
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Poccniicxaa Dezepans); 
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Asm Maypoxny Xacan Asan, JOKTOp HayK (KOMIIbIOTepHasd HHDKeHepuaA), FOWeHT, YHuBepcutet Amb-Haxpeitn (Mpax); 
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Abstract 
Introduction. The modern development of stamping aircraft manufacturing is inextricably linked with the assessment of 


the limiting capabilities of sheet blanks. However, the issue of defect-free forming of blanks made of aviation aluminum 
alloys is understudied. The importance of this issue is due to the fact that aluminum alloys are often used in the 
manufacture of thin-walled products for aviation purposes. During the implementation of shaping processes, various 
defects may appear, specifically, corrugation or unacceptable thinning. In this regard, the objective of the work was to 
construct a diagram of the limit deformations of the base aviation alloys and to conduct a comparative analysis of the 
limit deformation curves for these materials. 

Materials and Methods. Logarithmic deformations with the property of additivity were used to account for large 
deformations. The construction of the diagram of the limit deformations was carried out in the formulation of the 
deformation theory of plasticity. The issue of constructing a diagram of limit deformations was considered on the basis 
of the positivity criterion of the loading force derivative. In the area of negative values of the smallest major 
deformations, the Hill criterion was used to construct the limit deformation curve, and in the area of positive values of 
the smallest major logarithmic deformations, the Swift criterion was used. When constructing the limit deformation 
diagram, a power approximation of the hardening rule was used. 

Results. The curves of limiting deformations for the following aviation alloys were obtained: AMg6, DI6AT, AMg2M, 
1201-T, AMcM. According to the comparative analysis of the areas of safe forming, the values of deformations of the 
beginning of necking and their influence on the change in the position of the curve of the limiting deformation of blanks 
were compared: the greater the deformation of the neck formation, the higher the position of the curve of the limiting 
deformations. The concept of the Keeler's limit deformation diagram was described. Approaches to the construction of 
the Hill-Swift criteria used on the basis of the results of tensile testing of sheet specimens were presented. 

Discussion and Conclusions. Based on the constructed curves of limiting deformations for aviation alloys, AMg-6, 
DI16AT, AMg2M, 1201-T, AMcM, the following has been found. AMg2M alloy has the largest area of safe forming, 
1201-T alloy has the smallest one. That is explained by the difference in relative deformations of the beginning of neck 
formation. The conducted research made it possible to evaluate the possibilities of defect-free forming of thin-walled 
blanks made of basic aviation aluminum alloys. The use of the constructed diagrams of limiting deformation will 


provide predicting the appearance of breaks in the process of forming sheet blanks. 
Keywords: sheet stamping, forming limit diagram, logarithmic strains, Hill-Swift diagram. 
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AHHOTalna 
Beedenue. CoppemMeHHoe pa3BHTHe WITAaMMOBOYHOrO aBHaCTPOHTebHOTO TpOW3BOACTBa Hepa3pbIBHO CBA3AaHO C 


OWeHKOH TpeebHbIX BO3MO2KHOCTeH JIMCTOBBIX 3arOTOBOK. OHaKO MaJIOM3yYeHHbIM ABJIAeCTCA BOTIPOC 
OesqedeKTHOrTO (OPMOH3MeHCHHA 3arOTOBOK V3 aBHallMOHHBIX aJIOMMHMeBBIX CIWIaBoB. BaxKHOCTb JaHHoro Borpoca 
CBa3aHa C TCM, 4TO AJIIOMHHVeBbIC CIIaBbl JOCTATOUHO YaCTO HCMOb3YIOTCA MPH W3rOTOBICHHH TOHKOCTeCHHBIX 
v3eIHH aBMalMoHHoro Ha3HayeHua. IIpu peann3ayuu mpoileccos dopmMoobpa30BaHHaA BO3MO2%KHO MOABIIeHHe 
pa3IM4HBIX edekToR — rodpoobpa30BaHuA WIM HeAOMycTHMoro yToHeHHA. B cBA3H C 3THM LWesIbIO paOoTHI ABIAIOCh 
MOCTpoeHve AHarpaMMBbI IIpeAebHBIX JecbopMallMit OCHOBHBIX aBHalMOHHBIX CIJIaBOB HM MIpoBejeHHe CpaBHUTeIbHOTO 
aHaJIM3a KPMBbIX peesbHOro AeopMupoBaHHA JIA TAaHHbIX MaTepHasIoB. 

Mamepuaavi u memoooi. Ina yaeta Oonpurux ZedpopMaluii Opi MCHOM30BaHbI JOrapudMuyeckHe Aedopmann, 
oOnaqaioujve CBOMCTBOM ayyuTuBHOcTH. IlocrpoeHue uarpaMMbI MIpeyebHBIX edopMalHi mopMou3MeHeHHA 
IpOBOJHIOCh B MOCTaHOBKe AecdopMallMOHHON Teopuu WWiacTH4HOcTH. Bonpoc nocTpoeHua WuarpaMMBI UpeesIbHBIX 
TedopMallHit paccMoTpeH Ha OCHOBaHHH KPHTepHA MOMORKUTEMBHOCTH MpOu3BOAHON cHIbI HarpyxeHus. B oOnacTu 
OTPHUaTesIbHBIX 3HAYCHHii HaMMCHbIMX YaBHBIX epopMalyvi Wit MOCTpoeHHaA KPHBOM MpeyembHOrO 
eopMupoBaHHA MCHONb30BaJICA KPUTepHu Xia, a B 30H MOMOKMTCIbHBIX 3HAaYeCHMM TIaBHbIX HaMMeCHbILIHX 
jlorapupMuyeckHx JedopMauHu — xkputepuit Caudta. [pu noctpoenuu AMarpaMMbl UpeyesbHoro ZedopMupoBaHHna 
HCHOb3OBAIACh CTEMCHHAA ANMpPOKCHMallHA 3aKOHA YIPOUHeHHA. 

Pezyismamoi ucciedoeanua. Wonyuenbl KpuBble MpeAebHEIX WedopMalui AIA aBHalMOHHBIX criaBos: AMr-6, 
JI6AT, AMr2M, 1201-T, AMuM. CormacHo mpoBeqeHHoMy cpaBHUTeIbHOMy aHau3y oOsacTel Oe3zonacHoro 
(OpMOv3MeHeHHA, COMOCTAaBJICHbI 3HaYeHHA WecsopMallui Havana WelikooOpa3z0BaHHA HM UX BIIMAHHe Ha H3MeHeHHe 
MOJOXKCHHA KPHBOM UpezesbHOro WesopMupoBaHHA 3arOTOBOK: 4eM OosbUIe DedopMaluA wWieitKooOpa3s0BaHHA, TEM 
BBIIMe MOOKeCHHE KPHBOM MpeAebHEIX AetopMauMi. OnucaHa KOHIeMUMA WHarpaMMbl MpesesbHBIX TePpopMalnit 
Kunepa. II[peactapnenbt NoxxXogbI K MocTpoeHulo KpuTepHeB Xuma u Candta, Ucnob3yeMbIX HO pe3yIbTaTaM 
MCIIbITAHHA JIMCTOBBIX OOpa3I{OB Ha pa3pbIB. 

O6cystcdenue u 3akniouenua. Ha ocHOBaHHM NOCTpOeHHBIX KPHBBIX IIPeCIbHBIX TedopMalui Did aBMallMOHHbIXx 
cumapop AMr-6, JII6AT, AMr2M, 1201-T, AMuM spiacHuau, uro HanOobity1o oONacTh Oe30MacHoro 
cbopMou3MeHeHHA UMeeT crap AMr2M, HauMeHbutyro — criaB 1201-T, aro oObaACHAeTCA OTIM4MeEM OTHOCHTECIIBHBIX 
epopMaunwi Havana welikooOpa3sopanua. IIpopeqenHoe MccieqOBaHve MO3BOMJIO OLCHHTb BO3MO%KHOCTH 
Oe3qecdeKTHOrO (POpMOM3MeHeEHHA TOHKOCTCHHBIX 3afOTOBOK M3 OCHOBHBIX aBHall[MOHHBIX aJIIOMMHMeBbIX CIIVIABOB. 
IIpuMeHeHuve MOCTpoeHHBIX MarpaMM [IpeyebHOro etopMupoBaHHA MO3BOJIMT MpOrHo3HpoBaTb MOABIIeHHe 


Pa3pBIBOB B Ipolecce (bopMooOpa30BaHHA JIMCTOBBIX 3AfOTOBOK. 
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Introduction. Modern development of the aviation production is inextricably linked with the study of sheet 
stamping processes. One of the key problems in the task of shaping is the defect prediction, in particular, thinning, 
ruptures, corrugation. These issues are related to the assessment of the limits of the blank. By forming limit of a sheet 
blank, we will mean the ability of the material to deform to the required geometry without necking or destruction. 

To date, the problem of predicting defects of sheet blanks in the stamping process is solved using the following 
methods: 

— empirical, based on mechanical tests for simple stretching of metal samples, thin sheets and tapes, bending tests, 
as well as pipe testing methods for flaring and broaching; 

— theoretical and empirical, which are based on the use and dissemination of the test results of samples for uniaxial 
tension to other schemes of deformation of blanks; 

— theoretical, which are based on the use of criteria for limit deformation, specifically, in the manufacture of thin- 
walled products. The founders of these methods were J. Sachs, R. Hill, A. D. Tomlenov, V. D. Golovlev, G. D. Del, 
Z. Marciniak, A. D. Matveev, J. D. Lubahn [1-20]. 

It should be noted that the disadvantage of empirical and theoretical-empirical methods is the limited use of the 
results. 

The most important step in solving the problem of predicting defects of thin-walled products was the development 
of the concept of the forming limit diagram (FLD) proposed by S. P. Keeler [16—20], which today is generally accepted 
in solving sheet stamping problems. FLD diagrams are widely used in AUTOFORM and PAM-STAMP 2G CAE 
software systems. 

Experimental methods for constructing deformation diagrams are based on the test methodology presented in the 
works of Marciniak and Nakazima. It should also be noted that the issues of plastic destruction of sheet blanks were 
considered in [17—19]. In the last decade, interest in the construction of FLD diagrams has grown significantly. Most 
foreign studies are aimed at experimental construction of diagrams for specific materials, as well as numerical modeling 
of shape-changing processes using finite element methods [20-30]. Theoretical aspects of the construction of DFD 
diagrams and deformation diagrams of the third kind are presented in [31, 32]. 

It is important to note that the diagram of limit deformations provides estimating the beginning of neck formation, 
which ends with the destruction of the sample in the process of deformation. The FLD diagram binds the values of the 
main logarithmic deformations acting in the plane of the sheet. Forming limit diagrams enable not only to predict the 
destruction of the blank, but also to assess the presence of other defects, in particular, wrinkling, thinning, which, in 
turn, reduce the quality of the stamped part. The main zones of the FLD diagram are the zones of destruction, possible 


ruptures, safe forming, probable formation of folds, and wrinkling (Fig. 1) [31]. 
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Fig. 1. Forming limit diagram: 1 — zone of destruction; 2 — zone of possible ruptures; 3 — zone of safe forming; 


4 — zone of probable formation of folds; 5 — zone of wrinkling [31] 


Materials and Methods. When describing the processes of forming thin-walled blanks, the moment of transition to 


the plasticity stage is determined in accordance with the Huber-von Mises criterion [33]: 


6; = 07 + 02 — 0,02 = or, 


where o,, 02 — primary true stresses, at 0; = 0 due to the light gauge of the blank; o; — intensity of true stresses; 


or — yield strength. 
According to the deformation theory of plasticity, the relationship between the intensity of stresses and the intensity 


of logarithmic deformations is defined as: 


osH{or-tos) : 


where e;, €2, — principal deformations. 

Since in sheet stamping tasks, the shaping processes can occur in several transitions, therefore, large deformations 
are considered. The use of relative deformations is unacceptable. In this regard, the deformed state in (1) is presented in 
true logarithmic deformations. 


Intensity of the principal deformations: 


2 
e; = i je? +e? + e1e2. 


The ratio of the principal deformations and the primary true stresses: 


=2 pa 2 
a= 2 B=. (2) 


Based on ratio (1) and (2), the relationship between a@ and P is determined by: 


B = 2at1 (3) 


24a’ 


According to expression (2), the Huber-von Mises criterion has the form: 


= a/1-B+ PR = ap, 


and the intensity of deformations: 
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e = eV tata, (4) 

To assess the onset of necking, the criterion of positivity of the derivative of the loading force is currently used, the 
founders of which were G. Sachs and J. D. Lubahn [17]. According to the described criterion, the deformation of the 
sample is stable with a positive increment of the tensile force. The moment of unstable deformation with subsequent 


stretching starts at AP = 0 and continues at AP < 0 (Fig. 2). 


Fig. 2. Indicator diagram of uniaxial tension of the sample [33] 

The deformation diagram of aluminum alloys obtained by the results of the uniaxial tensile test is approximated by a 
power function in the theoretical analysis of the diagrams of the limit deformations of the shape change [18, 19, 31, 32]: 
ds = Ae" umn o; = Ae”, 
where os = P/F — true stress; P — tensile force; F— current cross-sectional area of the sample; e = In(1+ 

Al/l)), while A and n — coefficients of the power approximation. 

Using the criterion of positivity of the derivative of the loading force and the power approximation of the 
deformation diagram of the third kind, it is possible to obtain a relationship between the limit deformation of the sample 
at the time of the occurrence of the diffuse neck and the power approximation coefficient n under condition 
AP = 0 (Fig. 3) [32]: 

Cn = Cin = 1. (5) 

It is significant that relation (1) is performed using logarithmic deformations. In case of using relative deformations, 


only approximate equality is possible. 


Pmax 


O71 


a) b) c) 
Fig. 3. Occurrence of a diffuse neck under uniaxial tension of a flat sample: 


a — sample; b — stress state; c— deformed state [33] 


Consider the stress state of a plate to which two tensile forces are applied at the edges, i.e., the blank experiences 
biaxial tension (Fig. 4). 
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P, F, G2 
a) b) c) 


Fig. 4. Stress-strain state of the plate under biaxial tension: a — plate; b — stress state; c— deformed state [33] 


Using the criterion of positivity of the derivative of the loading force, we determine the deformation limit value e;,, 
at the moment when forces P, or Pz are maximum. For the case when F,/F, = const, P;/P, = const, at the moment of 
maximum tensile forces dP, = dP, = 0. H. W. Swift [19] proposed the relation for the limit deformation: 

(l= Bp 467)" 


4 — 3B — 3B? + 4p? 
According to expressions (2)—(4), the relation for describing the curve of limit deformations at various exponents of 


Cin = 4n 


power approximation in the hardening zone has the form: 
A(e, —n)(e2 + 2e,)3 — 3(e, — 2n)(e, + 2e2)(e2 + 2e,)? — 
—3(e, + 2n)(e, + 2e,)(e, + 2e,)? + 2(2e, + n)(e, + 2e,)? = 0. (6) 

It is known from the experimental and theoretical studies that after the occurrence of a diffuse neck, plastic 
deformation of the sample continues. Thereafter, a localized neck may occur, which differs from the diffuse one not 
only in size, but also in that its occurrence and development are carried out under conditions of flat deformation with 
intensive thinning of the sample. 

According to R. Hill [20], the limit deformation criterion is determined by the moment of formation of the local 
neck, at which the increment of the total force is zero. In this case, the relation for the limit deformation is determined 
by the expression: 
ayt/2 


_ 95, A-B+B 
= 2n 1+B 


(7) 


Cin 


Taking into account expressions (2), (4), expression (7) for constructing a limit deformation diagram according to 
the Hill criterion has the form: 

e, tez —n=0. (8) 

It should also be noted that the described approach is consistent with the finite element method, which has been 


widely used to construct FLD diagrams of various materials in recent years [21—29]. 

Research Results. In practice, we will construct the FLD diagram using two criteria [31], namely: the Hill criterion, 
which is used for e2 < 0 according to (8); the Swift criterion, which is used for e, => 0 according to (6). Let us consider 
the application of these relations to construct the curve of the limit deformations of aluminum alloys widely used in the 
aviation industry at known values of deformation of necking: ¢,, = 0.18 (AMg6), €, = 0.16 (DI6AT), &€,, = 0.2 
(AMg2M), €,, = 0.06 (1201-T), €,, = 0.1 (AMcM) [34]. 

Determining e,, from formula e,, = In(1 + €,,) and using expression (5), we obtain the value of the strain-hardening 
indices: n = 0.17 (AMg6), n = 0.15 (DI6AT), n = 0.18 (AMg2M), n = 0.06 (1201-T), n = 0.09 (AMcM). Then, the 


Hill-Swift diagram, according to (6), (8), will take the form shown in Figure 5. 
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Fig. 5. FLD-Hill-Swift diagrams for various aviation alloys: 1 — AMg6 alloy; 2 — D16AT alloy; 
3 — AMg?2M alloy; 4 — 1201-T alloy; 5 — AMcM alloy 


Discussion and Conclusions. According to the constructed curves of limit deformations, AMg2M alloy has the 
largest area of safe forming of the five alloys studied, 1201-T alloy has the smallest one. That is due to differences in 
the deformation properties of materials, specifically, the difference in the deformations of the necking onset. In AMg2M 
alloy, the relative deformation of the beginning of neck formation is 20 %, and in 1201-T alloy — 6 %. 

Thus, on the basis of data on the hardening curves, power approximation and deformation of necking for aviation 
aluminum alloys AMg6, DI6AT, AMg2M, 1201-T, AMcM, curves of limit deformations of forming have been 
constructed, providing for the determination of the zone of safe deformation of sheet blanks. The study results are of 
practical importance in solving sheet stamping problems for these materials to predict unacceptable thinning, raptures, 


and wrinkling of thin-walled blanks. 
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Abstract 
Introduction. The paper considers the issues of constructing mathematical models of the momentless equilibrium stress 


state of elastic convex shells using methods of the complex analysis. At the same time, shells with a piecewise 
smooth (ribbed) lateral surface were considered for the first time. The work objective was to find classes of shells for 
which it is possible to build meaningful mathematical models. 

Materials and Methods. Using the methods of the theory of the discontinuous Riemann-Hilbert problem for generalized 
analytic functions, a criterion for the unconditional solvability of the corresponding static problem for the equilibrium 
equation of a convex shell with a ribbed lateral surface has been obtained. This criterion, combined with the methods of 
the theory of generalized analytical functions, is a tool for constructing mathematical models of the state of momentless 
stress equilibrium of elastic convex shells. 

Results. A method has been developed for constructing mathematical models of the momentless equilibrium stress state 
of a convex shell under the action of a variable external load and the condition of stress concentration at the corner 
points of the median surface. The introduction of a vector parameter, as well as the concepts of “order of quasi- 
correctness” and “quasi-stability”, into the boundary condition provided both quantitative and qualitative comparison of 
mathematical models. Classes of shells have been found for which the description of mathematical models is given in 
terms of the geometry of the boundary in the vicinity of the corner points of the median surface. The obtained result, 
when applied to shallow convex shells, provides a geometric criterion of quasi-stability. It is established that for a 
shallow shell, which is not quasi-stable, the only adequate mathematical model is a probabilistic one. 

Discussion and Conclusions. The proposed method for constructing a two-parameter family of problems with a 
modified boundary condition makes it possible to simulate the momentless equilibrium stress state for fairly wide 
classes of convex shells with a piecewise-smooth lateral surface under a sleeve connection. At the same time, the 
developed algorithm for calculating the boundary condition index allowed us to answer the question of the existence of 
an adequate mathematical model for a shell with a side surface of an arbitrary configuration, and for shells of a special 
type (specifically, shallow or shells of revolution), to formulate a geometric criterion for the existence of a mathematical 


model. 


Keywords: thin elastic shell, generalized analytic function, Riemann-Hilbert problem, index of boundary value 


condition, mathematical model. 
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AHHOTalna 
Beedenue. B padote paccMoTpeHbI BOMpOcbl MOCTpOeHHuA MaTeMaTHYeCKHX MOJ[elIeli Oe3MOMeCHTHOrFO COCTOAHHA 


Hallps9KCHHOTO paBHOBeCHA YIIPyIHX BBbIMYKIbIX OOOOUEK C HCIOJb30BAHHEM MeTOOB KOMIVIeKCHOrO aHasn3a. IIpu 
3TOM BIIepBble PpaCCMOTPeHbI OOOJOUKH C KyCOUHO-riaqKou (peOpucToH) OoKoBoi moBepxHoctTsio. Lembo padoTsr 
ABJIANOCb OTBICKAHHEe KIaCCOB OOOOUeK, DJIA KOTOPBIX BO3MOXKHO HOCTpoeHHe COepxKaTeIbHBIX MaTeMaTHYeCKHX 
Mogesen. 

Mamepuaavi u memoovi. C tlomoub1o MeTOAOB TeOpHu pa3pHIBHOH 3afauH Pumana-[umpoepta Wa OOOOWIEHHBIX 
aHaJIMTHYeCKHX PYHKUMM WOyIeH KpHTepHi OesycNOBHOM pa3peLIMMOCTH COOTBETCTBYIOIMeH CTaTHYeCKOl 3aa4n 
JUIA YpaBHeHHA paBHOBECHA BBITyKIOH OOONOUKH C peOpucToi OoKOBOM MOBepXHOCTHIO. ITOT KpUTepHi B COUeTAHHH C 
MeTOqaMH TeopHH oOoOnméHHBIX aHarMTM4YecKHX CyHKUMH WpelcTaBiaeT coOol MHCTpyMeHT mMocTpoeHua 
MaTeMaTHYeCKHX MOJeIelH COCTOAHHA Oe3MOMEHTHOFO Hallps#KEHHOTO PAaBHOBECHA YIPYTHX BBIIYKIbIX OOOOUeK. 
Pezynomamot ucciedoeanua. Pa3padoTaH MeTOA MOCTpoeHuaA MaTeMaTHYeCKHX MOseseH Oe3MOMeCHTHOrO COCTOAHUA 
HallpswKEHHOrO paBHOBeCHA BBbIMYKIOM OOONOYKH Mpu elicTBHH NepeMeHHOM BHeWIHeM Harpy3KH uM ycIOBHu 
KOHI[CHTpalH Hallps2xKeHHM B YIIOBbIX TOUKAX CpeAMHHOM WoOBepxHocTH. BBeyeHue B TpaHW4Hoe yCIOBHe BEKTOPHOTO 
TlapaMetpa, a TakKe MOHATHH «MOpANOK KBa3HKOPPCKTHOCTH» HM «KBa3HYCTOMYMBOCTb»> MO3BOJIAIOT MPOBeECTH Kak 
KOJIMYCCTBEHHOe, TaK H Ka4YeCTBeHHOe cpaBHeHHe MaTeMaTM4eCKHX MOjeNeH. Hatiqenbr Kiaccbl OOoT0YeK, WA 
KOTOPBbIX OMMcaHHe MaTeMaTHYeCKHX MOJeNei WaeTCA B TepPMHHaXx TeEOMeTPHH paHUI{bI B OKPECTHOCTH YIJIOBbIX TOUeK 
cpeqMHHOK noBepxHocTn. Ilomy4eHHbI pe3syibTaT B IPHMeHeCHHH K MOJIOFMM BBIITYKIbIM OOO0UKaM NO3BONIAeT JaTb 
reOMeTpHuecKHH KpuTepHii KBa3HycTOW4MBOcTH. YcTaHOBeHO, 4TO WIA MoMorow oOOoNOUKH, He ABIAIOMelicA 
KBa3HYCTOMYMBOM, CEAHHCTBCHHOM ajeKBaTHOM MaTeMaTHYECKOM MOJeJIbIO ABJIACTCA BEPOATHOCTHAA. 

O6cystcdenue u 3akniouenua. VIpeqnaraembiii MeTO MocTpoeHua AByxMapaMeTpHyeckoro ceMelicTBa 3aqja4 Cc 
MOJM@HUMPOBaHHbIM TpaHH4HbIM YCJIOBHEM MO3BOJIACT MOJeMPOBaTb COCTOAHMe Oe3MOMEHTHOTO Halips»KeHHOTO 
PaBHOBeCHA JIA JOCTATOUYHO WHPOKHX KaCCOB BbIIYKIIbIX OOOOUCK C KYCOUHO-raqKOH OOKOBOM MOBepXHOCTBIO TIpu 
ycIOBHM BTyIOYHON cBa3u. [pH 3ToM pa3pa0oTaHHbI aIrOopHTM BBIYHCeHHA WHAeKca rpaHw4yHoro ycuIOBHUA 
MO3BOJIA€T OTBETHTb Ha BOMPOC O CYIJeCTBOBaHHH ajleKBaTHOM MaTeMaTHYeCKOM MOeIM WIA OOOMOUKH C OOKOBOH 
HOBCPXHOCTbIO MpOM3BONbHOM KOHurypauHu, a AIA OOoNOueK crewMabHoro Bua (HalpuMep, WosorMx WIM 


o0om04eK BpalljeHua) = ccbopMysIMpoBaTb reoMetpuueckuii kKputepnit CylecTBOBAaHHA MaTeMaTH4eCcKOH MOJeIH. 


Kioueebie cioea: TouKas yiipyrat oOon0uKa, oOoOMjeHHad aHamMTuyeckad PyHKuuA, 3anqa4a Pumana-luupbeprta, 


WMAACKC PpaHH4HOrO YCJIOBHA, MaTeCMaTH4eCKaA MOJCJIb. 


Baarogapuoctu. ABtop BbipaxkaeT MPH3HATCIbHOCTb KOJJICKTHBY Kadezppl «Teopetuyeckad u mpukiaquasd 
MexaHuka»> JIoHCKOrO rocy4apcTBeHHOTO TeXHHYeCKOTO YHHBeEpCHTeta 3a OOCYKeHHe pe3yiIbTATOB paOoOTHI, TakxKe 
ONaroqapHOCTb AOKTOpy TexHHYecKux Hayk, upodeccopy A. H. Conopbésy 3a mpodeccuoHasbuyro NoAWepxKKy IpH 


HanwcaHun padortst. 
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Introduction. The issues considered in the paper were studied and described in the works of I. N. Vekua [1, 2] and 
A. L. Goldenveizer [3, 4], who initiated the application of methods of the theory of generalized analytical functions to 
the momentless (membrane) theory of thin elastic shells and the theory of surface bending. To date, the final results in 
this direction have been obtained for convex shells with a smooth edge (i.e., with a smooth boundary of its median 
surface). The most significant of them — the correctness and quasi-correctness of the key problem with a static 
boundary condition with simply connected and multiply connected median surfaces — are consequences of the fact that 
the index of the corresponding Riemann-Hilbert problem is an invariant of the connectivity of the surface. The author's 
application of I. N. Vekua's methods to the problems of the theory of convex shells with a piecewise smooth edge'[5, 6] 
recognized a connection between the “geometry” of the median surface in the vicinity of its angular point and the 
picture of the solvability of the corresponding Riemann-Hilbert problems with a discontinuous coefficient of the 
boundary condition. The use of these methods in [7—9] allowed us to obtain an effective formula for the index under 
some additional geometric conditions on the angular points of the surface, and, as a consequence, the geometric 
criterion of quasi-correctness of the main boundary problem. 

The purpose of this work is to construct mathematical models of the momentless equilibrium stress state of convex 
shells with ribbed side surfaces based on the geometric criterion of quasi-correctness of the main boundary problem. 

Materials and Methods. Let S be a simply connected surface of a given class of regularity [7] with a piecewise 
smooth edge L = ie L, and angulat points p, (i = Tei, n) . Let us define on S along L, piecewise continuous vector 
field r= {a(s),B(s)} , which admits discontinuities of the first kind at points Piss where a(s) : B(s) (a +67 =f 
B => 0) — tangential and normal components, s — a natural parameter, and we denote by J the homeomorphism of 
surface S, to complex plane z=x+iy given in [9]. Let area D=J (S ) be the image of the surface when mapped to 


plane z with boundary F and angular points q, =J ( Ae Consider the following problem (task R ): to find in domain 


D , complex-valued solution w(z) of the equation: 

w(2)-B(z)w(2)=F (2), <eD. (1) 
satisfying the Riemann—Hilbert condition 

Re{A(Z),w()}=1(S). Cer, (2) 

where 

(6) =s(6)[B(E)1(6)-4(5)s(6)], oS) 
s(C) = 5S, (C)+is, (C) , t(C) =, (C)+it, (C) ¢ a Sale Bee £ (i = 1,2) — real-valued functions, complex- 
valued functions ¥(C) ; r(C) Hélder functions on each of the ares I’, = iG Wie s(w, +iw,), B(z) ; F(z) a 
functions of class L, (D) , r>2. Here, it is assumed that the solutions of class W’” at the points of discontinuity q; 


admit “integrable infinity”, i.e., they admit an estimate |w(z)| < const -|z -q,| , 0<a, <1. Following [10], the class 


of such solutions is denoted by H’. 
As is known [11], the static boundary value problem of the momentless theory for an elastic convex shell with a 


ribbed side surface in the mathematical formulation is problem R , where w(z) —acomplex stress function, F (z) = 
a complex-valued function of the external load. In this case, the condition weH is equivalent to the stress 


concentration condition at the corner points of the median surface. We will construct a mathematical model of the 
equilibrium state of the shell based on the results on the solvability of problem R for a surface of a special type 


' Tyurikov EV. Obobshchennaya granichnaya zadacha Goldenveizera dlya bezmomentnykh sfericheskikh kupolov. In: Proc. XIV Int. Conf. 
“Sovremennye problemy mekhaniki sploshnoi sredy. Rostov-on-Don. P. 290-293. (In Russ.) 
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(canonical dome [9]). To simplify the presentation, we assume that at each corner point, the direction of one of the arcs 
coincides with one of the main directions k, (k,) and call the dome 2-canonical (1-canonical). Problem R_ for 


canonical dome K is called canonical if the direction of field r at each corner point p is the direction of the 
generalized tangent [7, p. 46]. Let us introduce the notation: 8 — ratio of the corresponding principal curvatures k, , 
k, at point p (0 <6, < 1) : P(v,) — corner point p with internal angle v,, T(v,) — the set (sector) of directions 
of the generalized tangent at this point, T — the set of continuous on L vector fields r , defining the direction of the 


generalized tangent at each corner point p(y ae As established in [7], canonical problem R is quasi-correct for any 


field reT ,if n=2.Problem R is family R' of problems (1)-(3), each of which is given by the selection of vector 
field r. Following I. N. Vekua [1], we will call problem R's -quasi-correct in class H’, if it is unconditionally 
solvable in this class, and its solution depends on s real arbitrary constants (.s — the order of quasi-correctness). 
Definition 1. Canonical problem R is called quasi-stable with respect to the field of directions of a generalized 
tangent if problem R'is s -quasi-correct for any field r ET . 
Remark I. Owing to the theorem on the solvability of the Riemann—Hilbert problem for generalized analytic 
functions [11], problem R is quasi-stable if and only if index « of problem R is an invariant of field re7 . Here, the 


index of the corresponding conjugation problem with coefficient A(G)=A(G)At (C) is called the index of the 


problem. In the case of k = —1, the order of quasi-correctness is s=K+1. 
The technique [6, 8] of calculating the index of a boundary condition of form (2) is used below, as well as the 


concept [10] of special node p_ of problem (1), (2) or special point q =J ( P;) of discontinuity of the boundary 


condition (2). Following [6], the direction of the generalized tangent at a corner point is called special if the 
corresponding point of the discontinuity of boundary condition (2), (3) is a special node of problem R. 


Definition 2. Corner point P(v) is called an instability point of problem R, if sector T(v) contains a special 


direction. 


We introduce the notation: v , o — one-sided limits at corner point p(v) of the unit vector tangent to L, where 
vector o sets main direction k, on the surface at point p, and inside corner v is given by pair (-v,0). 
Statement 1. If the direction of vector 1 at point p(v) coincides with the direction of vector v, then point 


q=J (p) is a special node of boundary condition (2) if and only if: 


1 
Vv = arccos ——. (4) 
1+6 
If vector v is replaced by vector ©, then: 
v= arcetg/t, (5) 
Where t — the only positive root of the equation: 
2 2 
> eee I a E i. (6) 
S+t S +t K(1+?)+4E1 t 
a a 1 
Here E=|"5~ |, K =kk, (12, p.164]. Ar thar, arceos—— < arcctgV/t . 
+ 
We denote 
1 
0 = arccos —, = arccte/t . (7) 
1+6 
The consequence of statement | is statement 2. 
Statement 2. Corner point p(v) is a point of instability of problem R if and only if: 
O<vsSu. (8) 
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Corner point p(v) is a 1-type (2-type) point if condition 0<v<0@ C <vs 4 is met. As established in [8], k -type 


(k = 12) point is a point of stability. If p(v) — is the point of instability, then the only special direction 1, of the 
generalized tangent divides sector T(v) into two connected sets 7" (v) and T* (v) : 
Statement 3. The index of problem R in class H” is calculated from the formula: 
K=4+)"(4-«,), (9) 
i=l 
where n — number of corner points of the boundary, «,=m for point of stability p(v) of m-type (m = 1,2) , and 
K, =S incase re T°) (v) (s = (2) for the point of instability p(v) : 
Remark 2. Formula (9) for the index of problem R in class H, of bounded solutions, according to [10], takes the form: 
Kk=-4+)°(3-K,). (10) 
i=l 
The condition that the solution to problem R belongs to class H’ is the boundedness condition of the integral of 


the shell stretching energy [13, p. 83] in the vicinity of the corner point. 
From formulas (9), (10), statement 4 follows. 


Statement 4. If boundary L contains a single 1-type corner point, then problem R is for sure solvable in class H* 
and has a unique solution Vr ET ; if n2=2, n=n,+n,, where n, — the number of corner points of k -type (k = 1,2) , 
then, problem R is quasi-stable with respect to r ET with the order of quasi-correctness 5 =2n,+n,—3. 


If boundary L contains instability points, then according to (9), problem R in class H” is not quasi-stable with 
respect to field r eT . However, it is possible to distinguish such classes of fields with respect to which problem R is 
quasi-stable with different orders of quasi-correctness. Classes of such fields can be set by choosing the direction of the 


generalized tangent in only one of the sectors T) (v) ; T°) (v) at each point of instability p(v) ‘ 


Remark 3. As it is easy to see, formulas (9), (10) together with statement 4 are valid in the case: 
H<vst+o", (11) 
where w — a sufficiently small value given by surface S$. At this, the condition of smallness @ is not mandatory. For 
example, for umbilical (k, = k,) point p(v) , it is enough to put @ = - : 


Let us now consider 1-canonical dome K. In this case, the description of the special nodes of boundary 
condition (2) is also given by statement 2, with the only difference that in equality (4) and equation (6), value 


k, k a ; 
8° =— <1 must be replaced by 8° = = >1. Here, for corner point p(v) of instability of problem R , inequality (8) 
1 2 


takes the form: 


Sv <arcoos——. (12) 


Based on the obvious graphical analysis of equation (6), we conclude: 


k, 
1° value p= arcetg Vt is a function of two parameters, namely: = (5, k,) , ke (0; +00), if 8 =— <1, and 
1 


p=p(8,k,), k, €(0;+00) , if &° =>; 
2 


a ; . k 
2° Jim. w(6.k)=%, limu(d.k)=% vk, (0:4). 8 =, 


1 


Oa: To. : k 
3 jim (5.42 )=3 dim 1(5,k,)=0 Vk, €(0;+400), if Oa) 


y) 


Mechanics 


21 


http://vestnik-donstu.ru 


Advanced Engineering Research (Rostov-on-Don). 2023;23(1):17—25. eISSN 2687-1653 


4° for any fixed 6 from the right (left) semineighborhood of the unit, function w as a function of argument k, (k,) 


is a slowly changing function [14] in the sense that u(5,. k,) € (0 4 for any fixed 6=6,, k; € (0; +00) (j = 1,2) : 


Research Results. We submit for consideration a family of surfaces S_, t€ [0,¢) , where t — a small parameter, 
each of which is the median surface of a thin elastic shell V, from some family {v,} , where V, and S, coincide with 
shell V and surface S , respectively. We assume that for Vt € [0,¢) , the following conditions are met: 

1) surface S_. is a canonical dome of the regularity class W*”, r>2, with internal angles of magnitude v j at 
corner points p, (i — Last) of boundary L,, where surface §, is a 2-canonical (l-canonical) dome S with 
boundary L; 

2) ratios of principal curvatures KO ; Kf) of surface S$, at corner points, coincide with values & at the corner 
points of surface S ; 

3) KI) =k, +e, (t), where lime, (t)=0, k =kj” (7=1,2). 


Let J be the mapping of surface S$, onto complex plane z=x+iy, given above, D, =J (S.) — a family of 
simply connected regions with corresponding [., = J (Z,) and corner points g, =J ( D; ) , bounded in plane z. Consider 
a family of problems R,, TE [0,¢) ; 

w(z)-B,(z)w(z)=F.(z), zeD,, (13) 
Re{A,(Z)w(S)}=7. (6). Cer. (14) 


where functions B. (z) LA, (C) rae (C) are defined by the middle surface S_ of shell V, according to [8], and 1, (C) 


has form (3). 

It is assumed that the lateral surfaces of shells V, have common edges containing points p, (i = aeglhhs Note that 
problem (13), (14) Vte [0,) can be considered as a family of problems R® , each of which is given by selecting a 
vector field on S, along L,. 

Consider 1-canonical dome S. Statement 4 is true. 


Statement 4. If in each of the corner points p(v) of dome S, the following condition is met: 


v< siebndie— ee (15) 
1+6 


then VTE [0,¢) problem R. is quasi-stable in class H” with respect to the field of directions of a generalized tangent 


with the order of quasi-correctness 5 =3n—3. If at each point p(v) , the following condition is met: 


arcetg Vi +o) <v< +00, (16) 


where @° is determined by remark 3, value My) = Wy (e) is given by family S,, TE [0,¢), then problem R, is quasi- 


stable in class H’ with the order of quasi-correctness s=2n-3. 
For proof, it is enough to use statement 3, conditions 1-3 and the known properties [1, p. 97] of mapping J. 


Statement 4 remains valid for 2-canonical dome S , if conditions (15), (16) are replaced by the conditions: 
v <arcctg Jt —@>, (17) 
1 
arecos —— < v <= +0" (18) 
1+6 2 


respectively. 
Now let us refine the concept of a corner point of k -type (k = 1,2) of surface S. We say that point p(v) of 


boundary T’ is a 1-type (2-type) point relative to family S_, te[0,8), if conditions (15) and (17) are met for 


1-canonical and 2-canonical points, respectively, (conditions (16) and (18) for 1-canonical and 2-canonical points, 
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respectively). Consider canonical dome S, each corner point of which is a point of 1-type or 2-type relative to the 


specified family S. Then, the consequence of statement 4 is: 

Statement 5. If p and q — the number of points of 1-type and 2-type, respectively ( ptq= n) , then VTE [0,¢) 
problem R, is for sure solvable in class H™ and quasi-stable with respect to directions r of the generalized tangent 
with the order of quasi-correctness s =3p+2q-—3. Specifically, if p=0, then q=2. 

The statement remains valid if class H” is replaced by class H, of bounded in D, solutions, and order 
S=3p+2q-3 by order s=2p+q-3 provided 2p+q 23. Thus, problem R, is not unconditionally solvable in the 
following cases: p=0, q<2 and p=1, g=1. 


The case of an umbilical (5=1) corner point pP(v) , in which any direction is the main one, requires special 
consideration. In this case, according to” == 3° any direction of the generalized tangent is a special direction of 
problem R . However, even in this case, statement 4 remains valid if, for family {s,} , we consider t #0, and values v, 


tt are replaced by value - F 


Discussion and Conclusions. The results obtained can be used to construct mathematical models of thin and 
shallow shells of the positive Gaussian curvature with ribbed side surfaces. The most complete and advanced results of 
both linear and nonlinear elastic shell theory are obtained for thin shallow shells. A detailed discussion of the concept of 
“shallow shell”, as well as a description of various versions of the theory is given in [15, p. 29]. The linear theory of flat 
convex curves was developed by I.N. Vekua [2, 16]. Within the framework of this theory, the issue of the realization of 
the equilibrium stress state of a shallow shell with a ribbed side surface, when a static boundary condition of a general 
form is met, is reduced to problem R discussed above. 


Let P be a shallow shell [2, p. 164] with a ribbed side surface, S — its median surface with a piecewise 
smooth edge. We assume that at each corner point of surface S, the condition of strong shallowness k, ~k, is 
satisfied, which is equivalent to the following: 


5~1, (19) 


k, k 
where 5 — any of the ratios -, —.. Condition (19) means that any corner point p(v) should be considered both a I- 


2 1 


canonical and a 2-canonical point. Consider corner point p(v) , at which ver and a family of surfaces S_, 
te[0,€), given by conditions (1)-(3). According to properties 1°-4° of values 8, 1 and statement 2, the 
corresponding problem R, Vt E(0,) is not quasi-stable in class H’, i.e., point P(v) at v a — the point of 


instability of problem R,. In this case, the natural presumption is that in the right part of formula (10), value k, 


(1 <i< n) corresponding to point p(v) is a discrete random variable with possible values 1 and 2. Thus, if at each 


corner point p(v) of the median surface, the condition v = a is met, then, owing to formula (9), the order of quasi- 


correctness of problem R is a discrete random variable taking values m, m+1,..., m+n, where n — the number of 


corner points, m= 2n—3 for class of solutions H *,and m=n-3 for class HA. 
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The method proposed above can be used to construct mathematical models of the theory of thin shallow shells with 
ribbed side surfaces of any configuration. To do this, it is enough to use the results’ on the solvability of problem R for 
spherical domes with a piecewise smooth edge. Let us consider for definiteness, the median surface S under the 


assumption that all the corner points P(v) of the boundary are “outgoing”, that is, y < 7. In this case, a corner point on 


k 
a spherical surface is a special node of the boundary condition if and only if y= = (k = sarees) It follows that the 


“outgoing” corner point of the median surface of the shallow shell is an instability point if one of the following 
k 
conditions is met: y = (k =1,2). Thus, formula (9) for the index can serve as validation for the following 


hypothesis: 

if problem R for a shallow convex shell is unconditionally solvable in a given class of solutions, then its quasi- 
correctness order is a discrete random variable taking integer values K, K +1,.... K+N, where N — the number of 
instability points, K — the number given by a set of corner points and a selection of continuous vector parameter f . 

In conclusion, we note that the same reasoning can serve as validation for this hypothesis, but carried out for regular 
convex surfaces satisfying the condition of local symmetry [17] at corner points. 
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Abstract 

Introduction. The article studies the problem of validating the specified levels of reliability during experimental 
development of a complex technical series system. Such tasks arise when it is required to make a decision on testing the 
system as part of a larger one or on the completion of experimental development and the start of series production. The 
study is aimed at validating the reduction of the experimental development time. The task is to determine whether the 
hypothesis H, is accepted or rejected. 

Materials and Methods. To implement the research objective and task, a critical area described by the inequality was 
constructed based on the test results. The formulation of the requirements validation task was based on well-known 
approaches to testing statistical hypotheses. The conceptual apparatus of information theory, probability, and statistics 
was involved. The theoretical and applied literature on mathematical methods in reliability theory was studied. The 
particular tasks of the work were solved by known ways. Thus, the probability of obtaining the exact number of 
successful outcomes in a certain number of experiments was determined by the Bernoulli scheme. The exact confidence 
interval based on the binomial distribution was derived from the Clopper-Pearson relation. The theorem of 
A.D. Solovyov and R. A. Mirny made it possible to assess the system reliability based on the test results of its 
components. 

Results. Control rules adequate to the stage of experimental development (with insufficient data on the technical 
system) and the stage of series production were mathematically defined. The probability of a successful outcome when 
testing technical systems was represented by: 

— the probability of event for a system element; 

— confidence value; 

— required scope of tests 

In these terms, the null and alternative hypotheses and the corresponding reliability control procedures were 


investigated. Two provisions were considered. The first one provided using the null confidence 


hypothesis H, = {P 2 ae and an alternative H = ip < Pt to confirm the requirements (P,, y) for the reliability 
indicator of one parameter for any (P,, y). In this case, one trouble-free test was enough. The second provision 


considered a sequential technical system with independent elements that were tested separately from the system 
according to the Bernoulli scheme for one parameter. We considered the requirements for the system in the form of a 


set of values (P,, y) and the requirements for any of its elements (P,,;, y) . They coincided when the planned outcome 
of the tests corresponded to the cases when the ratio P= lim:P,=P,, was fulfilled, and the null alternative 


lsis<N 


hypothesis was selected from the theory of statistical hypothesis testing. 
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Discussion and Conclusions. The experimental development strategy should be implemented in two stages: the search 
and validation of the reliability of the elements through a series of fail-safe tests. In this case, the planned scope of tests 
of each element is determined taking into account the confidence probability, the lower limit of the confidence interval, 
and the requirements for reliability indices of one parameter of the technical system. If the use of the null confidence 


hypothesis is acceptable, one fail-safe test is sufficient to confirm the requirements for the reliability index. 


Keywords: experimental development, testing of statistical hypotheses, reliability of a technical system, null 
hypothesis, alternative hypothesis, hypothesis of distrust, confidence hypothesis, confidence probability, scope of fail- 
safe tests, binomial type test model, Bernoulli scheme, Clopper-Pearson equation, theorem of A. D. Solovyov and 
R. A. Mirny. 
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AHHOTalna 

Beedenue. Cratbs llocBaljena npoOseme NOATBepxKAeCHHA 3a aHHbIX ypoBHel HajexKHOCTH IIpH 9KCIepHMeHTaIbHO 
oTpaOoTKe cyIO%KHOW TeXHWYeCKOH CHCTeMbI C MOCIIeCAOBaTeIbHbIM COeMHeHHeM 9JIeMeHToB. Takue 3aauH 
BO3HHKAIOT, KOI a TpeOyeTca IPHHATb pellieHve OO MCIIbITAHHM CHCTeMBI B COCTaBe Oolee KpyMHO HIM OO OKOHYAHHU 
3KCHepHMeHTaIbHOH OTpaOoTKH U 3allycKe cepHiiHoro mpon3BoycTBa. Lenb uccneqoBaHua — OO0CHOBaTb 
cokpallleHHe CpOKOB 93KCIIepHMecHTaIbHOM oTpaboTKH. 3aqaya — ompeyeuIMTb, IPWHUMaeTCA WIM OTKIOHAeTCA 
runote3a H,. 

Mamepuaaoi u memoooi. Jina peamm3aiuu Wen u 3ayauu padoTsl 10 pesyibTaTaM MCIbITAaHHM CTpOMTCA KpuTHYecKaA 
oOslacTb, ONMChIBaeMad HepaBeHCTBOM. WopMyMpoBKa 3aa4uu MOATBepKTeHHA TpeOoBaHui Oa3upyeTca Ha 
V3BeECTHBIX TOAxoqaxX K MpoBepke cTaTHcTWY4ecKHxX YTuMoTe3. 3afeiicTByeTca MOHATHMHEIM allapat Teopun 
WHopMallMu, BepoaxTHOCTU u cTaTucTuKH. M3y4eHa TeopeTuyeckad MU MIpHKaqHad WMTepaTypa O MaTeMaTH4eCKHX 
MeToyax B Teopuu HaylexHocTu. UacTHple 3afa4u paOoThI pellieHbI W3BeCTHBIMH ciiocoOamu. Tak, BepOATHOCTb 
MOJYYeHHA TOUHOTO YMCA YCIeLWIHbIX HCXOJOB B ONpeesIeCHHOM KOJIM4eCTBE IKCIeEPHMeHTOB OlipeyeseHa MO CxeMe 
Bepuymu. TouHbiii OBepHTeIbHbIM WHTepBasl, OCHOBaHHbIM Ha OMHOMMaJIbBHOM pacripeyeeHuu, WomryyeH v3 
cooTHomeHua Knonmepa — IIupcona. Teopema A. JI. Conospesa u P. A. MupHoro no3BouIMsIa OLCHHTb HajlexKHOCTb 
CHCTeMBI 10 pe3yJIbTaTaM MCIIbITAaHHN ee KOMIMOHEHT. 

Pezyivmamvl ucciedoeanua. MaremaTuyeckw ofipesesieHbl mpaBvsia KOHTpOJIA, afleKBaTHbIe 9TAarly 
3KCHePHMeHTAaIbHOM OTpaOOTKH (pH HeOcTaTOUHOCTH MaHHbIX O TeXHW4eCKOM cCucTeMe) HM 93Tally cepHiiHoro 
TIpov“3BoAcTBa. BeposTHOcTs ycreliHoro Ucxoya Ip MCIIbITaHHM TeXHMYeCKNX CHCTeM IipeycTaBsIeHa depes3: 

— BCPOATHOCTH COOBITHA JIA JIIEMCHTAa CHCTEMBI; 

— 3HadeHHe JOBepUTeIBHOM BeposATHOCTH; 

— TpeOyeMbIii OObeM UcHbITAaHHi. 

C 9oTHx No3HyHi UccHeqOBaHbI HyeBad U asbTepHAaTHBHAad THMOTe3bI H COOTBETCTBYIOMNIHe HM TIpollexypbl KOHTPOIA 


HajJl@2KHOCTH. PaccmMoTpeHbi ABa TWOJIO#KCHHA. Tleppoe AONYCKaeT HCHOJIb3OBaHHe HyJ1eBort THHOTe3bI WOBepuA 


Mechanics 


27 


http://vestnik-donstu.ru 


Advanced Engineering Research (Rostov-on-Don). 2023;23(1):26—33. eISSN 2687-1653 


H, = P2 Pe c ambTepHaTuBoit H =|P<P,} WIA NMOATBepxKyeHHA TpeOoBaHuit (P,, y) K MoKa3aTesO 
HaJje@2*KHOCTH OJHOrO Mapametpa pH soOnIX (P,, y). [pH stom AoctTarouHo oOfHOrO Oe30TKa3HOrO HCIIbITAHHA. 


Bropoe TIOJO7#KCHHE PaCCMaTPHBaeT NOCIICTOBATCJIBHYFO TEXHHYCCKYIO CHCTeMy C N He3aBHCHMbIMU SJIEMCHTaMH, 
KOTOPbIe€ HCHBITBIBAIOTCA OTICJIBHO OT CHCTCMBI IlI0 CXCMC Bepryms1u ASIA OTHOTO HapaMetpa. PaccMoTpuM TpeOoBaHHa 


K CHCTeME B BHC COBOKYNHOCTH BeM4HH (P,, Y) HW TpeOoBaHHa K jOOOMy ee sIeMeHTy (P,,, y). OHH coBnayator, 


eCJI TWIaHupyeMbiii MCXOJ, MCIbITaHHM COOTBCTCTBYe€T CJIYdaiM BbINOJIHCHHA COOTHOMICHHA P= lim iP, =P 
> 1sis<N a 


ae 
a HyJIeBad alIbTepHaTMBHad THMOTe3a BbIOUpaeTcA U3 TEOPHH MpOBepKH CTATHCTH4eCKHX FuMOTEeS. 

O6cystcoenue u 3aknio“enua. CTpatTeruto IKCIePHMCHTaIbHOM OTpaOoTKH CleAyeT peau30BaTb B Ba STAI: MOMCK H 
MOATBepKAeHHe HaexKHOCTH 3IEMeCHTOB CepHei Oe30TKa3HbIX HCIbITaHHH. B 9TOM cilyyae TaHHpyembili oObeM 
WCIbITaHHH + KaKAOrO 9ICMeCHTa OlpeyeteTcA C Yy4eTOM JOBeEPHTeIbHOM BePOATHOCTH, HWKHeM TpaHHUbl 
J{OBCPHTeIbHOTO HHTepBasla UH TpeOOBaHHit K NOKa3aTeAIM HajlexKHOCTH OJHOTO MapaMeTpa TeXHM4YeCKOM CHCTeMBI. 
Ecam JONycTHMO VCHONb30BaHHe HyeBOH THMOTe3bI AOBePHA, WIA MOATBepKCHHA TpeOoBaHHi K MOKa3aTesO 


Haq@2#KHOCTH JOCTATOYHO OJHOTO Oe30TKa3HOTO HCIIbITaHHA. 


Kiuoueevie cioea: SKCHepHUMeCHTasIbHaAa oTpadotTka, WpoBepKa CTAaTHCTHYCCKHX THHOTC3, HaeoxKHOCTb TexXHH4ecKol 
CHCTeCMbI, HyJIeBadt THMOTe3a, ajIbTepHaTHBHaA THMOTe3a, THHOTe3a HeAOBepHA, THMOTe3a AOBCpHA, AOBCpHTeCJIbHaA 
BepOATHOCTh, o0beM Oe30TKAa3HBIX MCHbITaHHH, MOJICJIb MCHbITaHHH OMHOMHAaJIbHOrO TUIa, CXeMa Bepuysuin, 


ypaBHenue Knonnepa — IIupcona, teopema A. J. Conosbezsa u P. A. Mupxoro. 


BaarogxapHocrn. ABTOPBI BbIparKaroT OnaroyapHocTh pelleH3eHTaM, UbA KPHTH4CCKadA OUCHKa pe ACTaBJICHHbIX MaTe- 
pHasIOB HW BbICKAa3aHHble UpeWIOKCHHA TO UX YCOBepMIeCHCTBOBaHHtO cmocoOcTBOBaIH 3HAYHTCJIBHOMY ITIOBbIMICHHIO 


KaudecTBa HacTosilei CTaTbH. 


Asia wWaTupospannsa: Lapes O.1O., apes 10.A. Mogtsepxazenue noxka3atemei HayexKHOCTH Ip 9KCIepHMeHTasIbHOH 
OTpaOoTKe CIO%*KHOM TeXHHYeCKOM CHCTeMbI C MOCJICIOBATeIbHBIM COeMHeEHHeM 93sIeMeHTOB. Advanced Engineering 
Research (Rostov-on-Don). 2023;23(1):26-33. https://doi.org/10.23947/2687-1653-2023-23- 1-26-33 


Introduction. Rational methods of validating the specified reliability levels are of current concern for experimental 
testing of a complex technical system when a decision is made on the possibility of testing it as part of a larger structure 
or on the completion of experimental development and the start of series production. The same tasks arise during series 
production, if it is required: 

— to assess the readiness of the enterprise to produce series products based on the test data of the pilot batch; 

— to make a conclusion about the compliance of the products with the requirements of technical documentation, 
taking into account the operating data. 

The study objective was to obtain an acceptable solution for planning and reducing the scope of tests using methods 
of interval estimation of reliability indices of sequential technical systems. To achieve the stated goal, it was required to 
determine whether hypothesis H, was accepted or rejected. 

Materials and Methods. It is reasonable to formulate the task of validating the requirements in terms of the theory 
of statistical hypothesis testing [1—4]. Let P be the reliability of the technical system, P7 — some fixed (required) level 


for P. Prior to testing about P, three initial assumptions can be made: P = P,, P < P,, P > Pr. 
Each of them is called a null hypothesis if it is written as: 


Het =P). B= iP oP.) =P Si 


Set H,= P = Pp} contains only one element, therefore, hypothesis H, = iP = Pt is called simple. 
Hypotheses of the form H, = iP < Pp} and H, = {P > Pr are called complex. Along with the null hypothesis 


expressing a pre-formulated point of view, an alternative hypothesis H is specified expressing the opposite statement 
A, (A, OH #1). We use the conceptual apparatus of the theory of information [1], probability and statistics, 
applicable to solving such problems. Consider two aggregates of sets H, and H: 

ASAP SP SP Se Pa (1) 
B= {P 22;) f= iP <P (2) 
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Hypothesis H, in (1) will be called rigid, or the distrust hypothesis. Indeed, in case (1), initially (before the test), we 
proceed from a position of distrust of the quality level of the system. Reliability index P is assumed to be no higher than 


a certain fixed level P’,,. Hypothesis H, in (2) will be called the confidence hypothesis, since in this case, it is initially 
assumed that reliability index P is not less than some fixed value Pr. 

The meaning of values P’, and Pr is different. In (1), P’, — such a rejected value that at P< P’,, the system is 
considered unacceptable. In (2), P; — such a value that at P = P,, the system is considered acceptable for use. 
Obviously, P, > P’,. 

Research Results. Thus, it is required to determine whether hypothesis H, is accepted or rejected. In the theory of 
statistical hypotheses, a critical area is constructed for this purpose based on the test results. It is described by some 
inequality. Moreover, the null hypothesis (due to the initial confidence in it) is adhered to as long as it is reasonable 
from the point of view of the accepted level of significance a. Therefore, as is already clear that the reliability control 


procedure in case (1) will be significantly different compared to case (2). 
Indeed, we will further make sure that to reject hypothesis H, in (1) and accept hypothesis H = {P > P.) of 
meeting the requirement for reliability indices, a critical area (or condition) should be used 
P>P*,, (3) 
where P — the lower bound of the confidence interval for P at the value of the confidence probability y = 1 — a; 
P’,, — the lower bound of the rejected interval for P at the value of confidence probability y = 1 — a. 


In case (2), the condition should be used to accept hypothesis H, about the compliance of the value of parameter P 


to the requirement 
P>P,, (4) 
where P — the upper bound of the confidence interval for P at the value of the confidence probability y = 1 — a. 
In (3) and (4), the requirement for the reliability index of one parameter P is understood as a set of values (P’,, Y) 
or (P,, y), given before testing. 


Let one successful test be carried out under the conditions of the Bernoulli scheme. Then, using the Clopper-Pearson 
relations, we find the lower and upper bounds of one parameter with the value, for example, y = 0.95: 


P=(1-y)'=1—y = 0.05; P=1. 


Here, even for very moderate values P’, € [0.05; 0.95] , condition (3) is not fulfilled, while (4) is fulfilled at any P;. 


Let us show the validity of the accepted position. 

First_position. If it is permissible to use the null confidence hypothesis H, = {P 2 Prt with alternative 
H = {P < Py , then one fail-safe test is sufficient to confirm the requirements (P,, y) for the reliability index of 
one parameter for any (P,, y). 

If the initial hypothesis H, is the hypothesis of distrust from (1), then a significantly larger number of tests are 
needed. Thus, for m = 0, we get n= log (1 = y)/logP, >>1. This is quite fair, because, when testing hypotheses, they 
initially proceed from the validity of the null hypothesis H,. 

At the stage of experimental development, there are no sufficiently complete data, therefore, it is reasonable to use 
control rule (3). At the stage of series production, one can proceed from the confidence hypothesis and use a 
significantly easier control rule (4). This is acceptable if, according to the experimental testing, condition (3) was 
fulfilled. 


Consider a system consisting of N independent elements connected in series, which can be tested separately. Then, 


the probability of a successful outcome when testing technical systems: 


P={IP. (5) 
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Here, P; — probability of the same event for the i-th element. The requirements for value P are specified in the form of 


a set of values (P,,y). It is required to plan a procedure for monitoring the reliability of one parameter for each 


element of the system, i.e., to specify Vi e[l,n] a pair (P;,, y). 
N 
Due to the multiplication of P; in formula (5), ratio [[4: =P must be fulfilled. Besides, y, = y . 


As a result, the required scope of tests n; of each element increases dramatically, and even with fail-safe results of all 
tests, it becomes unacceptable. At P; = 0.9 and y = 0.95, N= 100 and m=0, Vi= 1,N: 


P,, ~ P,'” = 0.999, n,, =log(1 — y)/logP,, ~ 3000. 
This method of planning is logically contradicted by inequality n,, >n,, following from P,; > P, when m; = 0. It 
is clear that with fail-safe outcomes, the required scope of tests n,; of i-th element, conducted separately from the 


system, should be equal to the required scope of tests of system n,. To avoid this contradiction, the theorems of 
A.D. Solovyov and R.A. Mirny should be used [5—7]. Thus, when m=0, Vi= LN: 


<i /n 
P=min:P, =f (n, 0,7) = (1-y)". (6) 
Here, P — the lower bound of the confidence interval for the reliability index of a technical system by one parameter, at 
value y of the confidence probability; P; — the value of the lower bound of the confidence interval for the reliability 
index of the i-th element of the system with the same confidence probability; » — minimum number of tests of system 
elements; f(n,0,y) — the root of the Clopper-Pearson equation: 


m 


1-y=)>) (",) P’*q' = B(n, P, m). (7) 
k=0 
According to [8-12]: 
f(nng.y) <P < f(a, [ng],y). (8) 


aS 


Here, q =1-P;P= (1 m, /n,); n= min nj 3 [ng] — integral part of the product nq: f(n.nq,y) — root 
of equation J,(n, P,ng +1) =1-y. 
From (8), it follows: 
P=min:P,=Pm, (9) 


lsi<N ~~ 
where P ,, — minimum of P; at the value of confidence probability y. 
P,=min:P,> P.S(P,2 Pp) OP, 2 P,)...A(Py 2 Pr). 


m 2 
1si<N 


This is true not only for case (1), when m=0, Vi= 1,N, but also for case (2), i.e., for the outcome 


n= (N,N, ....Ny), m= (mM, O, O, ..., O) of the tests, where n and m — the vector of tests and the vector of failures, 
if n =n. 
At this, only one element that has been tested a minimum number of times fails. Indeed, in this case, ng =m, —an 


integer. Calculations allowed us to establish that (9) is also approximately performed at the outcome of tests n, m, Le., 


in case (3), if for pair (n,,m,), min :P, =P.,, is possible. 


In all the cases mentioned (9), the lower bounds are not multiplied, and the system degenerates into one weakest 
element. This provides validating the following position. 

Second position. Consider a sequential system with N independent elements that are tested separately from the 
system according to the Bernoulli scheme for one parameter. The requirements specified for the system in the form of a 
set of values (P,, y), and the requirements for any of its elements (P,,, y;) coincide if the planned outcome of the 
tests corresponds to the mentioned cases of fulfillment of ratio (9), and the null alternative hypothesis is selected based 
on (1) and (2). 

Consequence. In case of (3), execution of (9) — the planned scope of fail-safe tests (V—1) of elements, it is 
determined from: 

n,; 2 N= log (1 —y) / logP,. (10) 
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The volume of failures of the conditionally first element at m, is from the ratio: 
P.=f(m,m,7) =P. (11) 
Proof (11) is based on the fact that the condition P = P,, = P, is fulfilled if P,, =P, = P, and P, <Vie[2,N]. 
The latter relation is satisfied if (10) is satisfied, because: 
120, = =p oP. 


In all the cases considered, it was assumed that there was no information about P before the tests, except for the 


obvious fact P e[0, il. However, it may be known that P = P,,, where Py=0. Hence, P = Pe [ Pi, 1] . Value Py 


can be found from test data or calculations at the time of planning reliability tests. There is no method for determining 
Py yet, and its development is the task of future research. But if value Py is known, according to the full probability 


formula, you can find: 


P =P+ qgP. 


From here: 
PaPor jit 2), P= PoP oP). (12) 


The latter relation is fulfilled due to monotonicity of dependence P = (P —P,,) inP: 


Taking into account (12), ratio (2) will take the form: 


P, +0 +P,)f(nnqy) <P <P, + 0 -P,)f(aingly - (13) 
Let the condition for making a decision on the compliance of the technical system with the requirements (Pr y) 


still be (3), where the lower bound of the confidence interval is determined from (13), taking into account 


P= P € [ Pi. 1| . Then, from the ratio: 


P2P, (14) 
we find the planned scope of trouble-free tests by one parameter for each of N elements: 
n, =n, =log(1 = y) log(P, —P,) / A -P,). C5) 
Value n’y decreases in Py. It means, no (10) at Py = 0 and n’y = O at Py = Pr. 


N 
Example. The requirements for the reliability indices P= []F of the system are specified for one parameter in the 


form of a set of values (P7= 0.90; y = 0.95). Number of elements of the technical system is N= 100. According to 
available data, P; > Py = 0.70. It is required to find the planned scope of tests for each of N elements, if the fulfillment 
of reliability requirements is checked by condition (14). From (15), we find: 

n; =Nn’o = log (0.05) / log (0.90 — 0.70) / (1 — 0.70) = 7. 

Note that at Py, = 0n;>n’y = 29. 

Ratios (10) and (15) make it possible to plan the required scope of testing of the i-th element of the technical system 
with a certain sequence of experimental development of the system. The whole process of experimental development is 
divided conditionally into two periods: the search and validation of reliability requirements for a decision on the 
transition to the next stage of testing or on the acceptance of a technical system for series production. In the first period, 
improvements are possible, and it is reasonable to use models with a variable probability P of a successful outcome of 
the system test. 

The data obtained in the first period can be used to calculate value Py. 

In the second period, we deal with an established version of the design of the technical system and technological 
process. This makes it possible to use the binomial type test models discussed above with constant probability P. 

Let the first period of developing N elements of the technical system be completed, then, the question is raised about 
validating the requirements for the reliability index of the system by one parameter. It is reasonable to validate 
reliability if a positive decision is made only in case of a fail-safe outcome of the last series of tests for each of N 
elements. This strategy is convenient because it is based on the minimum possible number of tests of elements of the 


technical system and provides simple analytical solutions (10) and (15). 
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In general, it makes sense to investigate a strategy that allows for failures of elements during testing and is based on 
optimization of some objective function. But here, we restrict ourselves to considering only the mentioned strategy with 
fail-safe final series. 

Discussion and Conclusions. The results of scientific research allowed us to formulate the following conclusions. 

1. Even with a large number of N elements of the system, it is possible to plan the scope of their tests. In this case, 
the methods of interval estimation of reliability indices of sequential technical systems provide obtaining an acceptable 
solution, but only for one parameter. 

2. The strategy of experimental development of technical systems is closely related to the method of validating the 
reliability of elements. The rational strategy of experimental development provides for confirmation of the reliability of 
the elements after the search period through a final series of fail-safe tests. In this case, the planned scope of tests of 
each of N elements does not depend on WN and is determined by ratio (15), which includes the requirements (P7, y) for 
the reliability indices of one parameter of the technical system as a whole and Py. The scope of tests obtained by (15) 
for each element of the technical system, for any number of them, is small if moderate requirements (P; = 0.80 ...0.95; 
y =0.90 ... 0.95) are specified for a system of N elements, but only for one parameter. At the same time, the scope 
decreases with increasing value Py. 

3. In series production, when testing upgraded technical systems, it is possible to use a control method with null and 
alternative hypothesis change. If the use of the null confidence hypothesis is acceptable, then one fail-safe test is 


sufficient to validate the requirements for the reliability index. 
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Abstract 
Introduction. The development of the polar regions of the World Ocean contributed to an increased interest in studying 


wave processes in water bodies with ice cover caused by the action of a mobile load. In most papers of domestic and 
foreign scientists, the ice sheet was considered as an elastic or viscoelastic plate loaded with a rectilinearly moving 
vertical force. However, when modeling the impact of vehicles on the ice cover, it is of interest to investigate problems 
in which the force moves along a more complex trajectory. Therefore, this study aims at developing a method for 
studying the behavior of the ice cover under the action of a force moving along a trajectory of a complex shape, obeying 
an arbitrary law of motion. 

Materials and Methods. A method for solving problems of the action of an arbitrarily moving force on the ice cover of 
a reservoir of finite depth filled with an inviscid incompressible fluid is proposed. The ice cover was considered as a 
viscoelastic plate lying on the surface of a liquid in a state of potential flow. A concentrated force moving along an 
arbitrary closed trajectory and being periodic in time was applied to the upper surface of the plate. Hydrodynamic 
pressure acted on the lower surface of the plate from the liquid side. Due to the periodicity of the load applied to the 
plate, an integral time transformation was used to solve differential equations describing the behavior of this system. 
Further, using traditional methods, formulas were obtained for calculating stresses and displacements in the plate and 
components of the velocity vector of liquid particles. These formulas were presented in the form of an iterated integral. 
Numerical methods were used to calculate integrals. 

Results. Calculations showed that the deflection of the ice cover increased markedly with the growth of speed and 
tangential acceleration of the load movement. An increase in the relaxation time of ice and a decrease in the radius of 
the trajectory of the load also caused an increase in deflection. The distribution of the fluid particle velocity vector over 
the depth of the reservoir was calculated. 

Discussion and Conclusions. The proposed method has shown its efficiency in solving problems about the impact of a 
moving load on the ice cover of a reservoir. With its help, the influence of the curvature of the trajectory of motion and 
the mechanical properties of ice, the kinematic characteristics of the movement of the load on the deflection of the ice 
cover was investigated. The method simulates the impact of vehicles on the ice cover of a reservoir. The results of its 


application can be used under the construction of ice roads or airdromes on ice. 
Keywords: infinite ice cover, moving load, arbitrary closed trajectory, variable speed. 
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AHHOTauHA 

Beedenue. Ocpoenve MonApHElx paltioHop Muposoro oKeaHa CMOcOOcTBOBAaIO HOBbIINeCHHIO MHTepeca K M3y4eHHIO 
BOJIHOBBIX IIPOIeCcCOB B BOJOeMAaX C JICXAHbIM MIOKPOBOM, OOYCJIOBJICHHBIX JelCTBHeM MOJBYWXKHOM Harpy3KH. B 
OOJIbIIMHCTBe paOOT OTeYECTBEHHBIX HW 3apyOexKHbIX YYCHBIX JIeqAHOM MOKPOB paccCMaTpHBaJIcA KaK ylipyrad WIM 
BA3KOYUpyrad WJiacTHHa, HarpyKeHHat MpAMONMHeMHO BWKyWelicd BepTHKaIbHOM cHoM. Ognako mp 
MOJCJIMpOBaHHH BO3eHCTBUA TpaHCIOpTHBIX CpeACTB Ha JIeqAHOM MOKpOB MpecTaBlaeT MHTepec paccmMoTpeHne 
3aa4, B KOTOPBIX CHJIa JBMKeTCA 10 Oosee CO%KHOM TpaekTopun. IlosToMy WebIO JaHHOrO MCCeAOBaHUA ABJIACTCA 
pa3pa0oTKa MeTOJa UccueqOBaHHA NOBeeHHA JeqAHOTO MOKpoOBa MOA WelMCTBHEM CHJIbI, JBMOKyIUelca 10 TpaeKTOpHu 
CIOXKHOM POPMBI, HOLYHHAACh MPOVM3BOJIbBHOMY 3aKOHY JBYDKeHHA. 

Mamepuanoit u memoooi. Upeqioxer MetToy pelieHua 3aa4 O TelCTBHM JBYKYLWIelicaA TIPOM3BOJIbHBIM OOpa30M CHIIBI 110 
JI€IAHOMY MOKpOBy BOJOeCMa KOHCYHOM TIyOMHBI, 3allOJIHeHHOTO HeBASKOM HecxKHMaeMOl 2K AKOCTHIO. JleqAHOM MOKpoB 
paccmMaTpHBalicd KaK BA3KOYIIpyrad@ IWiacTHHa, JlexKalllat Ha MOBEPXHOCTH 2 KHKOCTH, HaxOsAllelica B COCTOAHHM 
MOTeHIMaIbHoro TeyeHit. Ha BepXHIOIO NOBEPXHOCTb IWJIaCTHHbI MpHIO%KeHa COCpesOTOYeHHa CHa, TBWOKyWaica 10 
MIPOH3BOJIbHOM 3AMKHYTOM TpacKTOPHH VM ABIIAIOWMIAACA MepHOAMYeCcKON Ho BpemMeHu. Ha HYKHIOIO HMOBEPXHOCTb TIaCTHHbI 
CO CTOPOHbI %KHAKOCTH JelicTByeT THqpoqMHamuyecKoe JaBsieHHe. B cusly MepHOAM4YHOCTH NpHNOKeHHOM K TWIacTHHe 
Harpy3kKH TPH pelieHun udpdepeHuMabHbIX ypaBHeHHi, OMMCbIBaIOWIHX MOBeyeHHe WaHHOM CHCTeMbI, MCTIOJIb30BaJIN 
WHTerpasibHoe ImpeoOpa3o0BaHve m0 BpeMeHH. JJaee, NPHMeHAA TpawMUMOHHbIe MeTOJbI, NOMyYMIM COpMybl [yA 
BBIYHCIICHHA HallpsKeHHM HW MepeMeLeHHH B IWIaCTHHe H KOMIOHEHT BeEKTOpa CKOPOCTH YacTHL %KHAKOCTH. ITU PopMysBI 
TIpecTaBMIM B Be MOBTOpHOro MHTerpaya. Jit BEIYHCIeCHHA MHTerpasiOB MCHOMb3OBaJIN YMCICHHbIC METOJEL. 
Pezyismamot uccredosanua. Pacuetel MoKa3asIM, YTO MporHO WeqAHOrO MOKpOBa 3AMeTHO YBeIMYMBaeTCA C POCTOM 
CKOpOCTH UM KacaTeIbHOrO yCKOpeHHA JBHKeHHA Harpy3KH. K pocty nporuba TakoKe IPUBOAAT yBeMyeHHe BpeMeHH 
pewlakcallMM jIbya MW yMeHbIeHHe payuyca TpaeKTOpHu ABwKeHHA Harpy3KH. PaccuuTaHo paciipeyeseHve BeKTOpa 
CKOPOCTH 4acTHL *XHAKOCTH M0 rryOuHe BosOema. 

Oécystcoenue u 3akio4uenua. I[penioxKeHHbIM MeTOA MoKa3aI CBOIO 9*pPeKTMBHOCTh MPH pelwieHHM 3ayjay oO TelicTBHH 
NOABWKHOM Harpy3KH Ha JIeqaHOM MOKpoB BofOemMa. C ero MOMOLIbIO HCCIeqOBaHO BIMAHHe Ha UporHMO ezAHOrO 
MIOKPOBa KPHBH3HbI TpacKTOpHH ABHKeHHA M MeXaHHYeCKHX CBOMCTB JIba, KHHEMaTHYCCKHX XapakTepHCTHK JBWOKeHHA 
Harpy3KH. Metoy, MojeMpyeT Bo3yelcTBHe TpaHCHOpTHBIX CpeJICTB Ha JIeqAHOM MOKpoB BoyOema. Pe3ybTaTHI ero 


TIpHMeHeHHA MOryT ObITb HCHOJIb30BaHbI TIPH CTPOHTCJIBCTBE JICQOBBIX JOPOr WJM aSpPOAPOMOB Ha JIb]Ty. 


KoroueBbie c1I0Ba: OeCKOHeUHBIN Je AHOU TIOKPOB, ABYMKYUWaACA Harpy3Ka, MpOW3BOJIbHaA 3aMKHYTaA TpaecKTOpHA, 


TlepeMeHHawA CKOPOCTb. 


BuaarogapHocrn. ABTOp BbIpakKaeT OnarojyjapHOcTb peueH3eHTaM 3a YKa3aHHble 3aMCYaHHA, KOTOPbIe MO3BOJIMJIM 


TLIOBbICHTb KAYCCTBO CTaTbH. 


Aaa unTupopanua. Tanadypoun A.B. Metoy peulenua 3aqauu 0 JBWKeHHM Harpy3KH MO jeqtHOMy MOKpoBy 
BoqoOeMa 0 CiIOKHOM Tpaextopuu. Advanced Engineering Research (Rostov-on-Don). 2023;23(1):34-40. 
https://doi.org/10.23947/2687-1653-2023-23-1-34-40 


Introduction. Recently, in connection with the development of the northern territories and areas of the World 
Ocean, great attention has been paid to the studying the behavior of the ice cover under the action of external load. A 
large number of works by domestic and foreign scientists are devoted to these issues. The impact of a pulsed moving 
load on a viscoelastic floating plate was considered in [1]. The effect of a moving load on the ice cover frozen to the 


channel walls was presented in [2, 3]. The fluctuations of the ice cover caused by a load moving at a constant speed 
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were considered in [4]. The behavior of a semi-infinite ice sheet under the action of a uniformly moving load on it was 
studied in [5]. The propagation of waves excited along a channel with an ice cover was discussed in [6]. In [7], 
nonlinear models were used in simulation. The monograph [8] presented the results of studies of surface waves in a sea 
with a floating broken and solid ice cover. The ice failure under the action of a moving load was considered in [9]. In 
the above works, the rectilinear load movement was mainly considered. The objective of this study was to develop a 
method for solving the problem of the impact of a force moving across the ice cover in an arbitrary way. Achieving this 
goal makes it possible to more accurately simulate the impact on the ice cover of vehicles, whose movement often 
occurs along rather complex trajectories and according to a complex law. 

Problem Statement. The oscillations of an infinite ice cover lying on the surface of a reservoir of finite depth under 
the action of a force moving arbitrarily along a closed trajectory were considered. The ice cover was modeled by a thin 
viscoelastic plate, whose mechanical properties were described by the Kelvin-Voigt model. The reservoir was filled 
with an incompressible liquid. 

Materials and Methods. The ice cover bending was described by the differential equation [7]: 

(1+ ,0;)A?W + c?ZW + kW + bd,®|,-9 =”, (1) 
where W(x,y,t) — ice cover deflection; D=Eh'{ 12( 1-1”) ); E — Young's modulus; 1 — Poisson's ratio; h — thickness 
of the ice cover; t, — relaxation time of deformations; A?=0,’+20,7 a7 +0,*; P,— ice density; p, — water density; c= 
p,h/D, k=p,.g/D, b= p,/D; P(x,y,t) — load acting on the ice surface; ®(x,y,z,t) — potential of fluid movement. The fluid 
behavior was subject to the equation: 

A®=0. 

Boundary conditions at the ice-water interface at z=0 and at the bottom of the reservoir z=—H (H — reservoir depth) 

had the form: 
0,W=0,Q|.-o, 0-D| ._., =0. 

It was supposed that force P (x, y, t) moved along an arbitrary closed trajectory y in an arbitrary way. It was assumed 
that P=P(s(t)), where s — arc coordinate measured from some fixed point of the curve y. If function s(t) was periodic 
with period T, then it was quite obvious that value P(s(t)), would also be periodic, and also with period T. The 
parametric assignment of the trajectory was taken as: 

x = X(t) 
t = Yo(t) 
where ¢ — time. 

Considering the steady-state process and applying the integral transformation with respect to variables x and y as 
well as the final integral transformation with respect to t on the interval [0; 7], we obtained: 

(1— i@)t) p'Wo— °°,” Wot kKWotbiaj,Do|--c=Po/D, 
8/@Mp — p’ Do=0, 
where Wy, ®) — images of unknown functions W and ®, p*=’+ a’; 4; a — parameters of integral transformations, 
respectively, by variables x and y; o,=2an/T, n=0, 1, 2,3, ... . Solving the second equation allowing for the boundary 
conditions, we obtained: 
Do=-i0,,Woch(p(z+H))/psh(pH). 

Then, the first equation followed: 

Wo=PA((c°+beth(pH Vp) ©, -p*-k+i OnTop’). 
The moving concentrated force was approximated by the function: 

P(x,y,t)=€° exp(-€°((x-X,(t)) + (y-yolt))”))/a 
where ¢ — parameter. 

Further, by performing quite obvious transformations, we got: 


1 i 2 ( 7,,cos(@,(t—-T)) %,Sin(@, (t—T seid 
w= sofll esl ( ( ve 2 ( ( 2, (ome Jos 


n=l r, 


n n 


rin=(c+bcth (PHYp) Oy P'-k,, ry= top", R'=5 +B, 
5=x(T-x, B= yo(t)V. 
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Knowing W, it was possible to determine formulas for calculating the components of the displacement vector and 
the stress tensor at any point of the ice cover using available ratios. 

In the numerical implementation of the proposed method, a problem arose related to the need to calculate integrals. 
One of them as a subintegral function had a strongly oscillating function for sufficiently large n, and the second was 
improper with an infinite upper limit. 

To calculate integrals from strongly oscillating functions, a quadrature formula was used, obtained by the cubic 


spline method [11]: 


b ; 1 N-1 eon =p 

10x 
Jef (x)dx~- - (M,.-M,). 
a @ j=l Si 


where h; — lengths of elementary segments into which the interval was divided [a; b]; S(x) — approximation of f(x) by 
cubic spline, M; = 5 (Xp): 
When calculating the improper integral, approximate ratio ss f(p)dp = f f(p)dp was used for a sufficiently 


large A. Value A was chosen so that the estimate of the allowable error | if i. f (p)dp| was small enough. 

Using traditional methods, it was possible to obtain the required estimates for the calculated quantities. For example, 
for deflection value of the ice cover, this estimate had the form: 
(c? +bcth(AH)/A)o, +A*(1+0,1, )+k 


no 


loa) 


fw(p) dp 


A 


_A2s4p2 
<4e2e 4 


[(? +beth(AH)/ A)@, — A’ -k] +[o,1,A*} 


w=) w(p)dp + J,” w@)dp. 
When calculating the sum of a series, the Lanczos sigma multiplier method was used to accelerate its convergence. 
Research Results. The calculations were carried out for the case of the action of a single concentrated force that 
moved along a closed curve, shown in Figure 1: 


“si =f; 
cae 
2sin( Tp) 
It was assumed that the thickness of the ice cover was h=0.25 m, Young's modulus of the plate material 
E=500,000,000 N/m’, Poisson's ratio u=1/3, ice density p=900kg/m*, water density p=1,000 kg/m’, reservoir 


depth H=5 m, e=2.5. The radii determining the shape of the force trajectory were assumed to be R;=15 m, Ro=9 m, 
R3=3 m (Fig. 1). 


ans, te[0:T]. 


Y,m 


4 -3 -2 -1 O 1 2 3 4 X,m 
Fig. 1. Trajectory of the concentrated force 


Figure 2 shows the change in the deflection of the ice cover (a), as well as the maximum values of normal stresses 


(at z=th) Sx and Sy (b and c) and at the movement speed of single concentrated force v =7.4022 m/s and tangential 


acceleration w,=0. 
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Fig. 2. Change in deflection and stresses of the ice cover: a — change in deflection; 
b — change in stress S,; c — change in stress S, 


On these graphs, the red dot indicates the place of application of force, and the blue color indicates the trajectory of 
movement. 

Figure 3 shows the dependence of the change in the maximum value of the deflection of the ice cover on the force 
speed rate (a) at the moment t=T/2. The position of the force on the trajectory at the moment is marked by a red 
dot (Fig. 1), while the tangential acceleration w,=0. Two graphs are given, one of which corresponds to the radius of the 
trajectory R,=15 m (solid line), and the second — to radius R;=5 m (dotted line). Graph b in Figure 3 shows the 
dependence of the maximum value of the ice cover deflection on the tangential acceleration rate at moment t=T at the 
point of the trajectory (2; 0) in Figure | at speed v=0. The graph in Figure 3 shows the dependence of the maximum 


deflection value on the value of the strain relaxation time. 
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Fig. 3. Dependences of the ice cover deflection on: a — speed; 
b — tangential acceleration; c — relaxation time 


The liquid behavior is shown in Figure 4. It illustrates the distribution of the speed vector of liquid particles over the 


depth of the reservoir at t=T/2, force speed v=3.7011 m/s and acceleration w,=0. 
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Fig. 4. Distribution of the speed vector of liquid particles over the reservoir depth 


Discussion and Conclusions. The influence of the movement trajectory curvature and the mechanical properties of 
ice on the ice cover deflection, speed and acceleration of the load movement was investigated. Calculations have shown 
that the deflection of the ice cover depends significantly on the movement speed of the force and acceleration of 
movement. 

However, the nature of qualitative changes in displacements and stresses in the ice cover caused by the action of a 
moving force with a change in speed and acceleration of movement changed slightly. 

Mechanical properties of ice, specifically, the relaxation time, had a noticeable effect on the deflection of the ice 
cover. 

The results obtained and the proposed method for solving such problems can be used in the construction of ice roads 
or airfields on ice. Moreover, the proposed solution method has shown its effectiveness and can be used to solve other 


similar problems. 
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Abstract 
Introduction. Data on the occurrence of initial failures obtained through testing on standard samples cannot always be 


extrapolated to real welded joints and structures. This is due to the difference between the concentrators in the joints, 
because after welding there is a significant structural and mechanical heterogeneity of the heat-affected and stress 
concentrator zone. Extended, deep concentrators are considered as crack-like defects, at whose vertices a volumetric, 
multiaxial stress state is formed. The paper addresses the issue of constructing critical diagrams of the onset of the 
limiting state at the concentrator vertex, which depends on the level of external load and the theoretical concentration 
coefficient. 

Materials and Methods. Analytical methods were used to study the stress state. The literature on the topic was 
analyzed. The features of proven physical models and patterns of behavior of materials were taken into account. The 
characteristics of steel alloys were taken from open sources and summarized in a tabulated form. Nonlinear equations 
were solved in MATLAB applications. The diagrams constructed by the authors enable to track the correlation of the 
dangerous level of the theoretical stress concentration factor and the level of external load. Curve Fitting Toolbox 
MATLAB was used to design the graphic part of the work. 

Results. The characteristic of damage from stress concentrators in welded joints was given. The crack propagation in 
the fusion zone was shown. The conditions stimulating and inhibiting destruction were indicated. The theoretical stress 
concentration factor a; was specified. It was shown how this indicator depended on the width, the height of the seam 
and the thickness of the welded part. Acute stress concentrators with theoretical concentration factor a, = 5...14 and 
more were studied. For this case, an approximating formula was given that took into account the maximum stress in the 
concentrator in the first half cycle, the initial deformation, and the load ratio. Through those elements, an indicator of an 
increase in maximum stresses was set depending on the number of loading cycles. The flow condition, the stress state, 
and the overvoltage factor, which took into account the increase in the first principal voltage for a combined stress state, 
were analytically shown. A model of the critical state at the apex of an acute stress macro concentrator was described. It 
was presented as the dependence of the relative stresses of the initiation of destruction 62°/o)2 on the concentrator. 
Possible variations of this model were analyzed. The dependences of relative values o2°/o9 2 on the theoretical 


concentration factor a? = 


a; were presented. To check the physical adequacy of this model, graphs were constructed 
that reflected changes in the relative stress of the external load at a critical state at the stress concentrator apex. The 
inevitability of bifurcation as a result of the studied processes was validated. Two directions of further development of events 
were indicated: brittle destruction and loss of stability of the stressed state with the transition to an increase in plastic 


deformations. The moment of bifurcation was defined as a critical state in the focus of the concentrator. 
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Discussion and Conclusion. The analysis and calculations performed within the framework of the presented scientific work 
enabled, in particular, to draw conclusions about the role of key factors of the processes under study. It was established, for 
example, that the operation of a steel alloy at a high theoretical stress concentration factor depended on the characteristics of 
the stress state. In a rigid state, it was possible to inhibit shear deformation and the onset of the limiting state at a lower value 
of the theoretical stress concentration factor. With the usual strength of steel (in comparison to high), a greater impact of the 
volume of the stress state on the value of the theoretical stress concentration factor was recorded. The probability of failure 
depended on the resistance of the material to the growth of a macrocrack. In future research, it is possible to refine analytical 


models and results, evaluate effective stress concentration factors. 


Keywords: welded joint, theoretical concentration factor, defects of welded joints, volumetric stress state, stress 


concentration, yield strength, macrocrack. 
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AHHOTalna 

Beedenue. Jjanupie 0 BO3HHKHOBCHMH HayaJIbHbIX pa3pyUIeHHii, MOMyYeHHbIe MCUbITaHHAMH Ha CTaHJapTHBIX 
oOpa3lax, He BCera MOXKHO SKCTpalloMpOBaTb Ha peasiIbHble CBapHble COCAMHCHHA HW KOHCTpyKUMH. ITO obycOBIeHO 
OTJIMYHAMH KOHIWCHTPaTOPOB B COCAMHEHHAX, T. K. WOCIe CBaPKH BOSHHKaeT 3HadHTeIbHad CTpyKTypHO-MexaHHyeckar 
HEOJHOPOAHOCTh 30HbI TePMHMYeCKOrO BIMAHHA WM KOHICHTpaTopoB HalipmKeHuli. IIpotaxenupie, rryOoKue 
KOHI[CHTpaTOpbl pacCMaTpUBaIOTCA KaK TPeLIMHOMOAOOHbIe WedeKTbI, B BeEPWIMHaX KOTOPbIX OOpa3yeTca OOBeEMHOE, 
CIO%KHOe HallpwKeHHOe cocTosHHe. PemlaeTca BOIpOc NOCTpoeHHA KPpUTM4YeCKUX WMarpaMM Hayasia BOZ3HUKHOBCHHA 
TipefesbHOrO COCTOAHHA B BepIHHe KOHUeHTpaTopa, KOTOpoe 3aBHCHT OT YyPpOBHA BHeIHel Harpy3KH U 
TeopeTHYecKOrO KOIPPUUHeEHTa KOHMCHTpalHH. 

Mamepuanvi u memoooi. J\na uccheyqoBaHva HallpsKeHHOTO COCTOAHMA 3a{eMCTBOBaIM aHasIMTHYeCKHe MeTOJBI. 
IIpoananu3supoBaHa IMTepaTypa moO Teme. YuTeHbI OCOOeHHOCTH MpoOBepeHHbIX (v3H4eCKHX MOJeeH MU 
3aKOHOMEPHOCTH MOBeeCHHA MaTepHalloB. XapaKTepHCTHKM CIIaBOB CTaJIM B3ATbI M3 OTKPbITbIX HCTOUHHKOB HU 
oOoOmjeHbI B Bue TabsuUbI. HermHerHble ypaBHeHHA peach B IpHkKaqHbIxX WporpamMax Matlab. Hoctpoenupre 
aBTOpaMH JMarpaMMbl MO3BOJAIOT OTCIICIHUTb KOppeAMIO OMacHOro ypoOBHA TeopeTHYecKoro KoopPunneHta 
KOHI[CHTpalMH HallpsKeHHu WM YypOBHA BHelwHei Harpy3KH. [na odopmieHua rpadbuyeckoi actu padoTsl 
ucnomb30BaH Curve Fitting Toolbox Matlab. 

Pezyismamoi ucciedosanua. Jlana xapakTepucTHKa pa3pyWIeHHii OT KOHUeHTpaTopoB HalipmKeHuii B CBapHBIX 
coeqMHeHuax. HariaqHo 10Ka3aHO pa3BHTHe TpelHH B 30HE CIaBIeHHa. YKa3aHbl YCNOBMA, CTHMyIMpyroulve U 
TOpMO3AlHe paspyuieHve. OnpeyereH TeopeTuyeckHi KoIPPHUMeHT KOHUeHTpalun HanpwKeHHu ar. IloKa3aHo, 
KaKHM OOpa30M laHHbI WOKa3aTeIb 3aBHCHT OT WIMPHHbI, BbICOTBI WIBa HM OT TOJUIMHbI CBapHBaeMOM jeTasMU. 
PaccMOTpeHbI OCTpble KOHI[CHTpaTOphI Hallps.KeHUM C TeopeTHYeCKHM KOIPUUMeEHTOM KOHIeHTpayHH a, = 5...14 1 


Oosree. Jina 9TOTO Cirydat UpHBOAHTCA allipoOKCHMupyroular cbopmyzia, KOTOpadA YU4HTbIBACT MaKCHMasIbHOoe 
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HallpsKeHHe B KOHIeCHTpaTope B HepBOM NOYyUMKIe, HCXOAHY!IO WedopMalHio U KOIPPuUNeEHT acHMMeTPHH WMKIA 
HarpyxeHnua. Uepe3 3TH 3JIEMeHTHI 3a{aeTCA MOKa3aTeJIb MOBbIMICHHA MAKCHMAJIBHBIX Hallps.KeHHM B 3ABHCHMOCTH OT 
YMClla UMKIOB HarpyxKeHua. AHaIMTHYeCKH MOKa3aHbl YCIOBHe TeKYYeCTH, HallpmKeHHOe COcTOAHHe u KoapPuynHeHT 
TlepeHalpwKeHHA, YUYHTbIBAIOWIMH WOBbIMIeHHe WepBoro TaBHOrO HallpsKeHHA AIA CiOXKHOTO HaliIps»KeHHOTO 
cocTosHHa. OnvcaHa MOJelIb KPHTHYCCKOLO COCTOAHHA B BEPIUHHe OCTPOrO MaKpOKOHIeHTpaTopa HallpsKeHHH. Ona 
TIpeycTaByleHa Kak 3aBMCMMOCTb OTHOCHTEJIBHBIX HallpwKeHUi 3apoxKeHHA paspyuleHuA Of°/092 OT KOHLeHTpaTopa. 
II[poaHanu3upoBaHbI BO3MOXKHbIe BapHaluu sTOM MoyemH. IpeyctaBeHbl 3aBHCHMOCTH OTHOCHTCJIbHBIX 3HAa4eHH 
of /0o,2 OT TeopeTH4eckoro Koop puyHeHta KOHYeHTpayun af = ar. La mposepKu duswyeckot ayeKBaTHocTu 
aHHOM MOjleIH MOCTpoeHbI rpatbHKH, KOTOpbIe OTPaxKalOT W3MCHEHHA OTHOCHTeIbHOTO HalipsKeHHA BHeLIHEeL 
Harpy3KH IIpH KPHTHYeCCKOM COCTOAHHM B BepIHHe KOHI[eCHTpaTopa HalipsxKeHHi. OOocHoBaHa HeH30e%xKHOCTb 
Oudypkalluu Kak pe3yIbTaTa UCcueAyeMbIx MpoweccoB. Yka3aHb! Ba HallpaBeHHA WaIbHelWero pasBuTHA COOBITHH: 
xpylkoe pa3pyllieHve HW MOTepA yCTOM4YMBOCTH HalipsxKeHHOTO COCTOAHHA C THepexOOM K pocTy mlacTH4ecKHx 
epopmanni. Moment Oudypkalnu onmpeyesieH Kak KPHTHYeCKOe COCTOAHHE B OUare KOHICHTpaTopa. 

O6cystcdenue u 3aKiio4uenua. AHaiu3 M pac4eTbl, BBIMIOJHeHHbIe B paMKax IpejcTaBIeHHOM Hay4HoH padoTHEt, 
MO3BOJIMJIM, B YaCTHOCTH, CileaTb BbIBOALI O POM KIIOYEBLIX (PaKTOPOB MCceyeMbIxX poweccoB. YcTaHoBJIeHoO, 
HallpHMep, 4TO padoTa cilaBa CTaIM IPH BbICOKOM TeOpeTHYeCKOM KOSPPHUMEHTe KOHUeCHTpalMu Hanps9KeHHH 
3aBMCHT OT XapaKTepHCTHK HallpsKeHHOTO cocTosHHA. II[pu x%XeCTKOM COCTOAHHM BO3MOX%KHO CiepoxKMBaHHe CABHTOBOM 
esopMalwu MW HacTyMieHHe MpeyeubHOrO COCTOAHHA Mp MeHBbUIeM 3HayeHHH TeopeTH4ecKoro KoIdpPuuNeHTa 
KOHI[eHTpalwHu HanpsKeHH. [Ip oObraHolt MmpouHOcTH cTalIH (B CpaBHeHHH C BbICOKOM) PuKcupyeTca Oobiee 
BIIMAHHe OOBEMHOCTH HallpwKeHHOTO COCTOAHMA Ha 3Ha4eHHe TeopeTMYecKoro KospPuuMeHTa KOHWeHTpalHu 
HallpwKeHHH. BepoaTHOCTb pa3pyWIeHHA 3aBHCHT OT CONPOTHBIAeMOCTH MaTepHalla pocTy MakpoTpelHHBl. B 
OyHYUIMX U3bICKaHHAX BO3MO2%KHO YTOYHeHHe aHasIMTHUeCKHX MOJeNel MU pe3yIbTAaTOB, OWjeHKa 93pdeKTHBHBIX 


KoIPPUUMEHTOB KOHUCHTpalHn HalpsKeHHH. 


KsoyeBbie CaOBa: cBapHoe coeqMHeHHe, TeopeTHYecKHH KoopPPUUMeCHT KOHLCHTpaluu, edeKTbI CBapHbIxX 


coeMHeHHH, o0beMHOE Halips»KCHHOS COCTOAHHE, KOHUCHTpalnA Halips»KeHni, TipewqeI TeKyyecTH, MakpOTpel{HHa. 


Baarogapuoctn. ABToppl BbIpaxkaloT MIpH3HaTeIbHOCTb COTPYAHHKaM JellapTaMeHTa MOpPCKOM TEXHUKH MH TpaHciopTa 
JIB®Y 3a KOHCYIbTaTHBHYIO TOMO. B pa3spaOoTKe TeMBI, a TaKXKe PCL|CH3CHTAaM 3a LCHHbIe 3aMe4aHHA 110 CTpyKType 


pykonucn. 


Aisa uuTHpoBanua. Monoxos K.A., Hopuxos B.B., Jabane3 M. OneHka MosBIeHHuA HadasIbHbIX pa3pyIIeHHi OT 
KOHICHTpaTOpOB HalipsKeHH B CBAPHBIX COCAMHEHUAX HM 3JIEMeHTaX KOHCTpyKuMi. Advanced Engineering Research 
(Rostov-on-Don). 2023;23(1):41—54. https://doi.org/10.23947/2687-1653-2023-23-1-41-54 


Introduction. Static and fatigue strength is reduced due to defects in welded joints. These can be: 
— cracks formed during or after welding; 


— stress concentrators (undercuts, incomplete fusion, pores, welding craters, high rippling, abrupt change in the 
shape of the weld, rolls, etc.) 

In the latter case, the determining factors will be: 

— shape, size of the concentrator and the position in the welded joint; 

— stress state indicator at the top [1]. 

Under cyclic loads, cracks often occur and develop in welded joints, provoked by stress concentrators [2]. The risk 
also depends on how close the stress state of the welded joint with the concentrator is to the occurrence of a macrocrack 
and its spread at operating o,, and cyclic loads o_,. Obviously, the latter depends significantly on the size, shape of the 
macrodefect and the stress state at its top. For example, an oblong and narrow macro stress concentrator is more 
dangerous than a rounded one. Probably, a crack is immediately formed at its top, and together with it, the initial 


concentrator will represent one long macrocrack [3]. The size of such a “total” macrocrack may be critical, which will 
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reduce the static strength [4]. It depends on its length and mechanical characteristics of the material at the top of the 
initial concentrator. 

The study objective is to develop an analytical model for assessing the dangerous level of stress concentration. The 
case in hand is about the formation of a crack in the top of the concentrator, a sharp decrease in the bearing capacity of 
welded joints and structural parts. 

Materials and Methods. Scientific research within the framework of the stated topic was based on known physical 
models and patterns of behavior of materials. The theoretical and applied literature was analyzed. Illustrative and 
reference materials were extracted from the sources. 

The stress state was studied by analytical methods. The relationship between the rigidity of the stress state and the 
value of the theoretical stress concentration factor was established. The results were presented in the form of diagrams. 
This visualization method made it possible to track the correlation of the dangerous level of external load and the 
theoretical stress concentration factor. 

To check the results of calculations based on the physical adequacy of the authors’ model, dependency graphs for 
the critical state were constructed. The paper used information about widespread structural steels of ferrite-perlite class 
in the state of delivery (steel 10, 22K, 50, St3sp, 37KHN3A, 30KHGSA, etc.). Their mechanical characteristics were 


obtained from open sources! and summarized in Table 1. The data was visualized in the Curve Fitting Toolbox 


MATLAB. 
Table 1 
Characteristics of steels 

Steel grade 0;, MPa 0,, MPa m Ox 
10 320 190 0.17 0.73 
15G 410 245 0.148 0.55 
St3sp 450 270 0.16 0.71 
22K 540 310 0.16 0.69 
50 680 350 0.16 0.62 
1OKHSND 540 390 0.132 0.71 
37KHN3A 1014 743 0.12 0.6 
30KHGSA 1750 1360 0.09 0.44 


Mathematical apparatus was applied to derive formulas. When solving nonlinear equations, the MATLAB 
application software package was used. 


Research Results 


Characteristics of damage from stress concentrators in welded joints. Since early 2000s, stress concentrations in 


welded joints and structures have been studied in relation to industrial tasks [2—3, 5—9]. 


Figure 1 shows examples of destruction of welded joints due to stress concentration. 


Saat 


a) b) 
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Fig. 1. Destruction of cross and butt joints” from stress concentrators: 


a, f— from the incomplete fusion zone, square butt joint; b — from undercutting, V-weld; c — from a sharp transition in the fusion 


zone; d — undercut; e — from the concentrator formed by a step transition at the edge junction [10] 


The destruction of cross joints with large non-fusion of the edges (indicated by arrows in Fig. 1a) is possible under 
the following loads: 

— static (as a result of embrittlement of the fusion zone); 

— cyclic (cause cracking at the tops of the concentrator). 

In [5], the stress distribution for butt, corner, cross joints and intermittent connections of ship structures is analyzed. 
In butt and cross welded joints, the area of the base metal adjacent to the weld is of particular importance. This is the 
weakest cross-section that determines the strength of the joint at variable stresses. It should be noted that the welds are 
always quite long compared to the thickness of the metal and the stress concentrator (undercut in the butt or cross joint). 
This prevents a shift in the area of the concentrator and creates a biaxial, and more often a multiaxial stress 
state [11, 12]. When considering the destruction of welds with concentrators (see Fig. 1 a, e, f), we take into account the 
concept of the occurrence of plastic deformation in the top region [5, 7]. In addition, it should be noted that the 
concentrators are adjacent to the fusion zone, which, after welding, experiences longitudinal residual tensile stresses. 
They affect the degree of rigidity of the stress state in the concentrator focus [1, 13-14]. All this prevents the passage 
of shear in the concentrator focus and can contribute to the occurrence of a minimal brittle macrocrack with a high 
theoretical stress concentration factor a, even from static load. 

In [9], indicative failures from welded defects in large-diameter pipelines are described. It is established that with 
long-term operation (more than 20 years), fatigue cracks appear in the undercuts of longitudinal welded pipe joints. 
Cracks in undercuts in the fusion zone, as a rule, branch out and spread in two directions (fig. 2): one moves away from 


the fusion zone, the other develops along the fusion zone (along the weld). This indicates a strong impact of the 


* Molokov KA, Novikov VV, Turmov GP. Osnovy raschetnogo proektirovaniya svarnykh konstruktsii. Vladivostok, 2019. T. 1. 204 p. (In Russ.) 
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mechanical, structural characteristics of the material’and mechanochemical heterogeneity on the propagation of initial 


destruction [13, 14]. On the other hand, the crack branching can complicate the process of further destruction and the 


transition of the crack into the main one. 


Fig. 2. Destruction in the pipe welded joint from the stress concentration in the fusion zone. Steel 14HGS. Single arrow shows 
fatigue cracks; double arrow — pipe burst edge; strokes — weld outlines. Here, b — weld width; h — weld height; 


5 — weldment thickness; y — angle of transition from base metal to weld [9] 
The theoretical stress concentration factor in the fusion zone can be determined from the formula that takes into 
account the concentration effect of a non-melting transition in this zone [5]: 
2 
a, =14+11h CMY ** (h/in(90°/7), (1) 
where y = A[90 exp(—p/po) + B] — angle of transition from base metal to weld; p — radius of transition from base 
metal to weld; pp = 1mm; A = 0.94...0.17; B = 0.8; b — weld width; h — weld height; 6 — weldment thickness 
(fig. 2). 
However, this dependence does not take into account the changed mechanical characteristics of the material. In 
these zones, hardening structures are usually present, whose mechanical characteristics can differ significantly (up to 30 


%) from the characteristics of the weld and the base metal. 
It is possible to indirectly control crack formation and the stress state in the region of the top of the stress 
concentrator through the concentration factor of their intensity. Accordingly, when a sample fails due to fatigue, one 


macrocrack is formed at the concentrator top, and in a corrosive medium — a group (fig. 3). 


) b) 


a 


Fig. 3. Fatigue failure of steel 45 from a V-shaped stress concentrator: a — in air; 


b — in water — in a corrosive environment. An increase of 134 times* 
After the formation of a group of macrocracks (fig. 3 b), the role of the macro-concentrator decreases, and further 


development of destruction can be suspended or inhibited. This was validated by experiments that indicated an increase 


> Molokov KA. Otsenka povrezhdennosti ferrito-perlitnykh stalei v usloviyakh malotsiklovogo nagruzheniya. In: Proc. Conf. “Nauka. Innovatsii. 
Tekhnika i tekhnologii: problemy, dostizheniya i perspektivy”. URL: https://elibrary.ru/item.asp?id=2375224 1 &selid=46181945 (accessed: 


31.10.2022). (In Russ.) 
4 Molokov KA. Otsenka povrezhdennosti ferrito-perlitnykh stalei v usloviyakh malotsiklovogo nagruzheniya. 
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in fatigue strength in samples with stress concentration in a corrosive environment”. Such a group of cracks loosens the 
material, changes the stress state above and below the surface at the top of the macro-concentrator, which prevents the 
formation of a leading crack that can develop further at a lower level of cyclic loads. 

Consider the situation on the surface or in the thickness of metal in massive bodies under stresses 0, at infinity in a 
plane-stressed state. In this case, it is reasonable to use the Kolosov and Inglis’s solution for an elliptical shape 
concentrator: 

Op = ATES 142 a/p, (2) 
where P 4 Oma, — radius of curvature and maximum component of stresses on the surface, respectively, at the top of 
the notch; a — semi-major axis, or half of the longest area perpendicular to the direction of the external load field. 

Dependence (2) is well consistent with engineering practice. In accordance with the accepted concept, the 
theoretical stress concentration factor at the notch is determined by the depth of the notch and the contour curvature 
radius at its top, but does not depend on the shape of the contour. 

Due to the development of information technologies, the evaluation of the theoretical stress concentration factor 
presents no great difficulty. The results of calculations of stress and strain concentrations in elastic and elastoplastic 
formulation of the problem can be obtained through the CAE finite element analysis [16—18, 12]. 

The experiments showed that at a certain stress concentration, external static load, and elastic-plastic stage of the 
material, an initial fracture crack appeared at the concentrator top at a certain distance from the surface. It rapidly 
propagated to an axis perpendicular to the direction of tension of the notched plates. The propagation is visible in the 
pictures (fig. 4). The process was fixed in the plates under a plane stress state and under a plane deformation. In the 
latter case, the appearance of a crack in the thickness of the metal under the surface was detected using an ultrasonic 


flaw detector. 


a) b) c) d) 


Fig. 4. Initiation of a subsurface crack at the top of the concentrator under static tension: a — macrocrack initiation; b, c — size 


increase in depth and towards the surface; d — crack exit to the surface [19] 


Thus, under the given conditions and a sufficiently acute stress concentrator, a zone of a multiaxial stress state was 
formed, in which the continuity violation was possible. There would be a higher level of the onset of fluidity, and the 
tensile strength in a weakened place due to the concentrator could increase to 20-60 % (it depends on the material) [11]. 

Stress concentration and stress state conditions. Consider rather sharp stress concentrators with theoretical 
concentration factor a, = 5...14 and more. As a rule, the defects mentioned above act as such concentrators. At the 
top of the concentrator, the stress state can be considered the same as in plane deformation. In [18], the results of the 
analytical solution and experimental data for sharp concentrators under plane deformation were presented, the region of 
admissible values [a,] and the radius at the top [p] for steel 09G2S were established. However, the impact of its 
mechanical characteristics was not described. Moreover, it is not clear whether the analytical solutions are based on 
such initial parameters of 09G2S as the hardening coefficient m, yield strength 692, ultimate strength og, critical 
fracture narrowing @,, etc. It is known that ductility and strength affect significantly the level of allowable stress 
concentration and the occurrence of brittle fracture at the top. 


For hardening materials, cyclic stresses at the top of an acute concentrator aggravate the possibility of brittle 


> Sergeev NN, Sergeev AN. Mekhanicheskie svoistva i vnutrennee trenie vysokoprochnykh stalei v korrozionnykh sredakh. Vologda: Infra- 
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fracture, therefore, we will take such ay as a hazard criterion, at which brittle fracture is likely to occur under static load 
conditions. The increase in maximum stresses depending on the number of loading cycles N is given by the 
approximating formula: 


Omax = Omax + A ‘lg(N), (3) 
where Oinax —- Maximum stress in the concentrator in the first half cycle; A —coefficient depending on the level of 


initial deformation and the asymmetry coefficient of the loading cycle at the top of the concentrator. 

We accept the obvious assumption that the appearance of brittle fracture may be due to the following reasons: 

— an increase in the intensity of three-dimensional stretching in close proximity to the top of the concentrator; 

— an increase in the resistance to plastic deformation of the material near the top as a consequence of shear resistance 
in acute stress concentrators. 

The latter circumstance turns out to be the stronger the higher the loading speed and the sharper the stress 
concentrator. 

Consider the elastic-plastic plane stress at the top of the concentrator and the elastic state in the gross section: 
Omax = ApOy > 09.2, Oy © O92 and Oy < 0g. Here, 69. — conditional yield strength of the material under uniaxial 
tension. At a, = 7.3, stiffness of the stressed state is close to the state for a crack. Safety coefficient m;,4p for a crack- 
like defect with a plane deformation differs slightly (within 10 %) from 77f,p for a crack, i.e., for a state where there is 


an approximate equality Mpapn ~ Mpap: 


1-(60/or )?/ap-V1-(60.2/nr/oo)* 1-(69.2/nr/o0)? 
1-(6onp/or)2/a%l1—(69.2/00)2 ¥1-(60.2/60) * 


Here, n; — yield strength margin coefficient, og — local strength of the material at the top of the crack or sharp 


(4) 


concentrator. 

In [11], the term “local yield strength” is applied to acute stress concentrators when the area of the concentrator 
adjacent to the top, experiences a complex stress state (CSS). Obviously, this term denotes the material flow stress on 
the concentrator circuit. There is plane stress state on the surface at its top, and a triaxial CSS appears under the surface. 
There, stiffness of the stress state increases sharply, therefore, the initial destruction of continuity is formed at some 
distance from the surface (fig. 4). It is established that the onset of local fluidity in the stress concentrator zone does not 
coincide with the level determined by the calculations based on strength criteria, particularly, according to von Mises— 
Huber—Hencky. As the stress concentration increases, the difference increases between: 

— yield strength 69 2, theoretically calculated according to this criterion; 

— experimentally determined stress value of the local flow of material 05. in the region adjacent to the top of the 
concentrator. 

The results of experiments on plane samples of various steels and alloys [11, 16] show that ratio 06./09.2 is well 
approximated by linear dependence on ar. Then, the yield condition can be written as: 


0; = 092(0.9 + 0.1a,), (5) 
where o; — von Mises stress intensity. 


The authors limited their experimental studies to values a; = 10, and the tests were carried out only on plane 
samples with centrally located concentrators. In case of plane deformation, stiffness of the stressed state turns out to be 
slightly higher. Therefore, it is logical to assume that fluidity at some values a; would occur no earlier than at the 
plane stressed state. And it is likely that the slope of the straight line according to equation (5) could be somewhat 
different [11]. However, we focus on this dependence. 

For the moment of the beginning of yield, assume that with the transition from elastic deformations to elastoplastic, 


the ratio of the second and third main stresses to the first 06; = Omax does not change. This has been proven 
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experimentally. 

We equate to (5), the intensity of the stresses of the beginning of von Mises yield according to the fourth (energy) 
theory of strength, and introduce the coefficients of the relationship between the main stresses under plane deformation: 
0, = 0;/D; 02 = qo,; 03 = Ur (14+ q)o, [15]. Here, nu; = 0.5 — Poisson's ratio in the plastic region, and D — 
overvoltage coefficient, which takes into account the increase in the first main voltage in the case of CSS. The equation 
obtained with respect to q has two solutions. The first is for tensile 0. component, the second — for the compressive 
one. Tensile component o, increases stiffness of the stress state in the case of CSS. After transformations and 


compressions, we can express it like this: 


q=1-S2(14+2), (6) 

For the limiting case 0, = 69 2, it is possible to neglect the change in the radius at the top of the concentrator and 
the onset of some initial global fluidity for the net cross section. Then, the only value a, satisfying the stiffness of the 
stress state equivalent to the crack is 7.33. This validates approximate equality (4) obtained earlier. For other values, as 
can be seen from (6), proportionality q will vary depending on value ar and oy. 

For a crack, coefficient D is determined from the equilibrium condition in the elastoplastic region and is calculated 
for plane deformation from the formula: 

d= me (7) 
where m — power hardening coefficient; 1 — Poisson's ratio in the elastic region. 

For ferritic-pearlitic steels, this coefficient takes values from 0.22 to 0.26. 

Multiplier before D in formula q = 1 — 2D/V3 for a crack is 1.156, and q for ferrite-pearlite steels is 0.73. Indeed, 
when substituting a7 = 7.3 (6) for og = 092, we get q = 0.73. Next, we find a7, if 092/04 > 1, at which the same 
stress state is realized as at ar = 7.3 and 6092/0, —> 1. 

In [20], the impact of a7 on the effective stress concentration coefficient K, is investigated, and for a typical 
aluminum alloy, it is shown that, in the range a; from 7 to 13, global extremum of the value K, is observed. It is quite 
probable that such a significant maximum value may indicate brittle fracture at the initial stages of cyclic loading or 
discontinuity at the top of the concentrator. In this case, the resource will depend only on the further ability of the 
material to prevent the propagation of macrocracks. 

Model of critical state at the top of an acute macro-stress concentrator. Almost all defects in the welded 
connection of pipelines, hull ship structures, etc., create a stress concentration. An analytical model developed on the 
basis of the concept of “local stress concentration factor” is presented in [9]. In addition, a simple dependence was 
obtained based on the Neuber formula linking the theoretical generalized concentration factor and the concentration 
factors of the intensity of elastic stresses and deformations. As a result, for the critical value of plasticity in the 


concentrator focus, the following dependence was obtained: 


opar!/ tm) 


Eo = (8) 
where a, — critical value of the theoretical stress Seen factor; og — modulus of rupture; EF — modulus of 
elasticity; m — exponent of power hardening. 

Dependence (8) gives such high values that are rarely found in practice. For example, with critical logarithmic 
deformation for steel 50, the critical concentration factor is 26, and for high-strength steel 37KHN3A with 
Op = 1.014 MPa — more than 18. It can be concluded that such plasticity is unattainable until continuity is broken in 
the first loading cycles. This is not surprising, since it is known that the maximum plasticity value for samples with 


stress concentrators is significantly lower than for samples without stress concentration. In addition, €, should be 


determined not only by the properties of the material, but also by the conditions for the development of plastic 
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deformation before destruction. Similar results are given by the dependence used in [9]: 


Eve ‘oO 

pasp’Cpasp 9 
a 9) 
S.A. Kurkin used it for a comparative assessment of the material sensitivity to stress concentration. Here, ar 


ar = 


characterizes the concentrator in which the crack occurs at specific load o,,. Next, we will use the concept of “the most 
probable initiation of a macrocrack in the concentrator” (born crack) «8° for this case. The true deformation and stress 
at the moment of occurrence of the discontinuity at the top of the concentrator are denoted respectively by 
Epasp ANd Opasp- They can be determined from the generally accepted formulas Eas) = Exp = In{1/(1 — ~,)] and 
Opasp = Sorp = On(1 + 1.4¢,). Here, Son) — true stresses of destruction at uniaxial stress state. Note that the fluidity in 
the concentrator occurs according to (5). Replace oy = 692 and 092 with 092. As a result, we obtain for the critical 


states (a; = 8°) relative stresses of the nucleation of destruction 62°/o9.2 from the concentrator: 


ob — Vv E-€pasp’Opasp (10) 


Son a2©(0.9+0.102°)o9 2° 
Note that for a very small radius of the concentrator top (p < 10d,, where d, — average grain diameter), crack 


propagation is limited at constant a7. As shown in [6], it depends on d, of steel. Formulas (8) and (9) are attractive to 
an engineer because they provide comparing the material sensitivity to stress concentration and using the initial data of 
the mechanical characteristics of steel. 

We can get a slightly different model that takes into account the CSS in the area of the top of the macro- 
concentrator. The author [12] uses the strength criterion of brittle fracture in the structural material 0; > Sorp, where 
Sorp — Stress of normal fracture in the approximation to a uniaxial stress state. Significantly, this characteristic of the 
material does not depend on the temperature of its test. We use this criterion in a slightly different way — for the 
moment of occurrence of a discontinuity in the area of the top of the stress concentrator. Suppose that S,., is achievable 
in the case of constant energy equality at CSS and uniaxial stress state, i.e., with loss of plastic stability and transition 
from a volumetric stress state (when 0; is small) to a uniaxial stress state with an inevitable increase in deformations. 

We use power approximation of the deformation diagram within the framework of the deformation theory of 


plasticity. We equate true stresses 0; of the deformation diagram in the elastic solution function of the problem m 


Oo; = F(o”) to Sorp. AS a result, we get: 


matt 1 wy)? 
5,= fom (0) = Soup = op (1 + 1.49%), (11) 
where o” ) __ uniaxial stresses in the elastic solution to the problem. 
Writing (11) with respect to o” - we obtain an expression for the elastic solution to the stress concentration problem 
on the other hand: a” = Ar * Oy. Thus: 
a41 
m 
Op: Oy = eer (12) 
m 


Here, 06, — nominal workload on the welded joint. In hard points of welded structures, joints or under overloads, it 
can reach the value of the yield strength of the material. Let us imagine a special case under conditions of cyclic loading 
with hardening of the material for the moment of formation of the born crack discontinuity — failure initiation. After 


transformations (12), we write down a simplified formula: 
1+m 
abe = [eee re (13) 


90.2 
where «2° — value of the theoretical stress concentration factor at which the nucleation of destruction occurs in the 


region of the top of the macro-concentrator; 09.2 — yield strength of the material at the top of the concentrator, which 


can be increased according to (5) in the case of CSS. 
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At the moment of equality ap = a%°, the nominal stresses must be equivalent to the nucleation (destruction) 
stresses o2°. For this situation (13), we can write in relative values 07/09 and construct the dependences of the 


relative values of the stresses of fracture formation on the theoretical concentration factor a?° = ar: 


-2m 
2-(ap°)"*t .op-(70K+5) 


be 
ar +9 


on = (14) 

Unfortunately, it is not possible to write (14) with respect to a%° at of /o9> = 1. However, a?° can be found 
numerically for specific load 02°, or diagrams of the danger of discontinuities in the form of a macrocrack in the focus 
of the concentrator can be constructed. 

Dependences are analytically obtained for calculating critical theoretical stress concentration factors related to the 
mechanical characteristics of the material and a given external static load under the condition of a volumetric stress 
state for concentrators in welded joints. It is shown that the volume of the stress state in the concentrator focus affects 
significantly the value of the critical theoretical stress concentration factor. It is established that the complex stress state 
in the concentrator focus can be controlled by the geometric characteristics of the concentrator itself and its location 
relative to the external stress field. 

Experimental data of steels (see Table 1) are obtained from literary sources. These are ferrite-perlite materials for 
which formula (14) has shown good agreement. Dependence (14) is attractive because standard mechanical 
characteristics are used as initial data. It can also be used to assess the risk of defects in the fusion zone (fig. 1,fig. 2), 
where hardening structures are formed, whose properties differ significantly from the initial characteristics of the 
welded steel. 

Solutions for (9) and (14) are used more often for the most brittle steels at low nee However, in the area of high 
stress concentration, critical ratios o2°/o 92 for different steels will be very close in values. This can be explained by 
the fact that with the CSS, a stress state similar to a dangerous case is created for all steels in the top region. In the 
region of small stress concentrations, solutions to model (9) show higher critical stress values 02°. Note that the 
average value «3° at o2° © oo is ~8.5. For steels with different mechanical characteristics, it does not change as much 
as according to the calculated results of model (9). 

Figure 5 is constructed to validate the calculation results based on the physical adequacy of model (14). These are 


graphs of critical state dependences 02°/o92 from a°. 
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Fig. 5. Dependences of changes in the relative voltage of the external load 02/09 from aS for the critical state at the top of the 
stress concentrator. Solid curves — formula (14); dotted curves — formula (10) 


Model (9), which does not take into account the volume and increase in yield stresses at the top of the concentrator, 


gives a?° = 91 at stresses 02° = 0.509 for steel 10, which is unlikely. If we take into account the volume in this 
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model, then «8° = 33, which is significantly less. Taking into account the volume of the stressed state, the critical 
plastic deformation is reduced to the occurrence of destruction during static tensile of the samples. This was validated 
by experiments with various steels and alloys®. Moreover, the researchers obtained less ductility of the steels on 
corroded samples. Here, it was shown that a volumetric stress state was created at the tops, and the volume prevented 
the free flow of plasticity. The plasticity before the destruction of these samples was 40 %, and the yield strength 
increased to 27 % [21]. It can be concluded that deeper concentrators cause a greater increase in the yield strength. The 
results of indirect experiments allowed us to assert that under hydrostatic load, the fluidity is not fixed, in this case the 
material only embrittles. However, in the focus of the concentrator, with an increase in the external load, the volumetric 
stress state cannot remain stable enough for the development of a higher volume level. Therefore, bifurcation is 
inevitable. Further, the development can go in two directions: 

— fragile destruction or its beginnings in the focus; 

— loss of stability of the stress state and, as a consequence, the transition to a sharp increase in plastic deformations. 

Probably, the moment of bifurcation should be considered a critical state in the focus of the concentrator. Perhaps, 
this moment depends significantly on the depth of the concentrator, its relative length and depth to the thickness of the 
part, orientation to the external stress field, as well as on the mechanical characteristics and extent of the concentrator 
focus. All this requires additional research and more precise solutions. 

Discussion and Conclusions. Defects and design features of welded joints (weld groove geometry, incomplete 
penetration, undercuts, pores, hardening structures, etc.) reduce static and fatigue strength. Stress concentrators are 
typical for welded joints, having a long and shallow shape relative to the thickness of the parts being connected 
(undercuts) or long and deep (incomplete penetration, etc.). All of them reduce fatigue strength. The analytical models 
obtained in the presented work allowed us to draw a number of conclusions. 

1. The operation of the material at high values of the theoretical stress concentration factor depends on the stress 
state and its rigidity, as well as on mechanical and structural characteristics. The latter may differ from the parameters 
of the source material, since the tops of the concentrators may be located in the fusion zones of welded joints. 

2. A rigid stress state in the stress concentrator focus can cause an increase in the yield strength. In this case, the 
passage of shear deformation is restrained, and the onset of the limiting state is achieved at a lower value of the 
theoretical stress concentration factor and is characterized by discontinuity at a constant value of the external load. 

3. Comparison of high-strength steels with ductile and ordinary strength steels at the same level o7°/o 2 . In the 
second case, the formation of the volume of the stress state will have a more significant effect on the change in the 
critical value of the theoretical stress concentration factor, at which the nucleation of destruction occurs in the region of 
the top of the macro-concentrator «° (fig. 5). 


4. Comparison of models (14) and (10) allows us to draw a certain conclusion. Taking into account the rigidity of 
the stress state and the increase in the yield strength at o2°/o)> > 1 from below, (10) gives estimate a° with a 


margin of safety, if we are talking about 02/092 = 0.8 (30KHGSA). For less durable steel (14), the margin of strength 
at 02/69 2 = 0.5 (37KHN3A), etc. Thus, it is preferable to take into account the volume and use model (14). 


5. The dependences on the diagram of the theoretical stress factors are more densely grouped in (14), than in ( 10) 
due to very similar stress states in the region of the tops of acute stress concentrators. Therefore, the mechanical 
characteristics of steels have a secondary effect on a°. 

6. With an increase in the theoretical stress concentration factor and a decrease in «2° at the same level 02° /o9> of 
the external load, the influence of the structural factor in steel on the occurrence of initial destruction (formation of 
discontinuity) reduces. However, for further destruction, the resistance of the material to macrocrack growth is 
essential. After its occurrence, other criteria of mechanics and kinetics of destruction are triggered. 

The presented study can serve as a prerequisite for the development of analytical models to assess the residual life of 
welded joints and structures exposed to cyclic loads. Further research will presumably refine the analytical models. The 
authors will check and test the results through CAE modeling and proceed to the evaluation of effective stress 


concentration factors. 
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Abstract 
Introduction. Polyethylene is the most widely used material in various fields of the national economy, and products 


made of it have essential advantages, such as lightness, insolubility in organic solutions with quite satisfactory strength. 
However, the mechanism of its destruction is quite complex and depends on the working conditions and substances, 
which are in contact with it. The research purpose was to study the polyethylene resistance to oil and associated water 
under the static and dynamic conditions and at room temperature. 

Methods and Materials. The research was carried out on a laboratory device for passing various liquids (oil, associated 
and distilled water) in polyethylene pipes, assembled by the author of the article. While working, methods of statistical 
and dynamic data processing were applied. 

Results. Based on the results of previous experiments, graphs were plotted for the change in the weight and volume of 
immersed granules over time. They have shown that polyethylene has a quadratic dependence, and diffusion for the 
three liquids studied in this work (distilled water, accompanying water and oil) is described by Fick's law. This indicates 
the fact that the rate of liquid diffusion through polyethylene is the key factor. 

Discussion and Conclusions. The obtained results have shown that the rate of liquid diffusion through polyethylene is 
the key factor. Immersion in oil has a greater impact than immersion in associated or distilled water due to the presence 
of salts. It has been found out that the relative change in the weight and thickness of the polyethylene pipe walls through 
which oil passes is greater than those through which the associated water passes. Moreover, the microscopic cross 


sections images in the samples before and after the tests have confirmed the obtained results. 
Keywords: Polyethylene, associated water, solid materials, salts, distilled water. 
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AHHOTalna 
Beedenue. WlonuzTuNeH aBiaeTca HawOosee WIMPOKO IpHMeHACMBIM MaTepHasIOM B pa3IM4HbIX OOAacTAX HapowHOro 


XO3AHCTBA, a W3AesIMA U3 Hero oOsaqaror CYIeCCTBCHHbIMH JOCTOHMHCTBaAMH, TaKHMH Kak JIGKOCTb, HepacTBOPUMOCTh B 
OpraHHdecKHxX pacTBOpax IIpH BIIOJIHe YAOBJICETBOPUTeIbHOW TIPOUHOCTH. OHako MeCxaHH3M ero pa3pyilieHisA 


HOCTATOUHO CJIOHKCH HU 3aBUCHT OT ycIOBHit padortsi MW KOHTAKTHpyFOuIMx C HUM BelLIeCcTB. Lemp MccIeqOBaHui 
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3akI0Yuaacb B W3yYCHHH yCTOMYMBOCTH NMomMITHIeHa K HePTH U COMyTcTByIOWleH Boe MpH cTaTM4ecKOM U 
J{MHaMH4YeCKOM PpexKUMe HW KOMHATHOH TemMilepatype. 

Mamepuaavi u memoooi. UccneqoBaHua MpoBOAMIMCbh Ha JaOopaTOpHOM ycTpoiictTBe Aya MpOmycKaHuA pa3zIMYHBIX 
*UAKOCTeH (HeTb, COMYTCTBYIOWad HW JMCTHIWIMpOBaHHad BOJa) B MOJM3ITHIICHOBLIX TpyOax, CoOpaHHOM aBTopaMHi 
ctatbu. IIpu padote Obi IpHMeHeHbI MeTObI CTATHCTHYeCKON U WHHaMH4ecKON OOpaooTKH WaHHBIX. 

Pezyismamot uccnedosanua. papuku u3MeHeHHa Beca HM OOBEMa MOrpy2KeHHbIX TpaHyJI BO BPeEMeHH, MOCTPpOeHHEIe 
TO pe3yIbTaTaM TPOBeCHHbIX 9KCICPHMCHTOB, TOKa3aJIM, YTO AIA MONMSTMeHa XapakTepHa KBaflpaTHuHad 
3aBHCHMOCTb HU AMdpy3vA WA Tpex WccleqyeMbIxX B aHHOM padOoTe %xMAKOCTeH (AMCTHIMpoBaHHad BOA, 
CONYTCTByIOINat BOLa U HedTb) OMMCbIBaeTCA 3AKOHOM Muka. DTO yKa3bIBaeT Ha TOT (pakT, YTO CKOpocTS Aupdpy3uu 
XKUKOCTH Yepe3 NOMMITHICH ABJIACTCA HaMOOJIee BaXKHbIM (PaKTOPOM. 

O6cystcoenue u 3akmio“uenuA. DKCIICPUMCHTEI UM MOUYYeCHHbIe pe3yJIbTAaTbI WOKa3ayIH, 4TO ckKOpocTs Audpy3un 
XKUKOCTH Yepe3 MOMMITHICH ABJIAeTCA HaHOosee BaxKHbIM dakTopom. IlorpyxxeHHe B He*Tb OKa3bIBaeT OosbUIee 
BIMAHHe, YM TOrpy2KeHHe B COMYTCTBYIOMYy!IO WIM JMCTHIWIMpoBaHHy!0 BOy H3-3a MpvcytcTBuA comet. brio 
oOOHapyxKeHO, YTO OTHOCHTeIbHOe H3MCHEHHE BeCa H TOJILIMHbI CTCHOK HOJIMITHIICHOBBIX TpyO, HO KOTOPbIM IIpoOxOJuT 
HeTb, OobUe Tex, 10 KOTOPbIM IIpOXOAHT COMYyTCTBYIOWat BONa, a MMKPOCKONMYecKHe H300paxKeHHA Cpe30B B 


oOpa3lax TIO HW Tocue UciibITaHHi TIOATBCPAYIIN MOJTYACHHbIC pesyJIbTarTBl. 


Koro4ueBble C10Ba: NIOJIMSTHIICH, COUNYTCTBYFOMLAA BOa, TBCpAble MaTepHaJIbl, COJIM, HHCTHJVIMpOBaHHad BOTA. 
BuaarogqapHocrn. ABTop BbIpaKaeT NPH3HATCJIBHOCTh PCHCH3eHTaM H peakiMu *KypHaJla 3a WCHHbIC 3aMedaHHA. 


Aaa uuTupopannsa. Autuoac U.P. Ucnonp30Bpanue MeTOa KOHC4HBIX IICMCHTOB JIA MOJeMpoBaHHA pesepByapa 
BBICOKOTO JaBJICHHA H3 MOIMMepa, ADMMPOBAHHOTO yIJIepOAHBIM BOJIOKHOM. Advanced Engineering Research (Rostov- 
on-Don). 2023;23(1):55—65. https://doi.org/10.23947/2687-1653-2023-23-1-55-65 


Introduction. Polyethylene is one of the main elements of the olefin family, which is usually characterized by a 
waxy appearance, chemical inertia and physical properties deterioration at high temperatures. Since polyethylene is a 
partially crystallized material and does not go into secondary chemical reactions with any liquids, it has no solvents, and 
its water absorption at ambient temperature is relatively low. 

Polyethylene does not dissolve in organic solutions, but it expands at a temperature of more than 70 °C and 
dissolves in granular aromatic hydrocarbons, it is not affected by oil, fats, kerosene and other hydrocarbons. It is also 
resistant to aqueous solutions of acids, salts and alkalis [1—5]. 

Polyethylene is oxidized by nitric acid, which causes a deterioration in the characteristics of the products obtained 
from it. Oxidation occurs in the presence of ultraviolet rays at ambient temperature, causing final brittleness and 
deterioration of the mechanical and thermal insulation properties [6—8]. The scientific interest is focused at the practical 
interaction study of polyethylene and oil, which is a mixture of hydrocarbon materials consisting mainly of compounds 
containing carbon and hydrogen as well as some elements mixed with hydrocarbons, such as sulfur, nitrogen and 
oxygen, in addition to metal elements in simple proportions, such as iron, nickel, arsenic and vanadium. Along with 
them, chlorides salts of land-based minerals as well as their sulfates are present in oil. Moreover, oil contains water in 
certain percentages, depending on its type and production conditions, and this percentage can increase up to 30 % 
relative to the total volume of oil, since it is in the form of saturated solutions of sodium chloride, magnesium sulfate 
and calcium. Carbon dioxide and hydrogen sulfide gases dissolving in water increase its volume. Acid and oxygen, 
reacting with iron and manganese dissolved in the associated water, form insoluble products. The associated water 
contains a small percentage of organic chloride salts, not exceeding 50 parts per million [9, 10]. 

The impurities contained in crude oil (solid materials — salts — water — chemical additives) cause serious risks to the 
oil refining equipment, its operating time and efficiency. Polymer materials are known to be destroyed in a completely 
different way than metal materials; therefore, when producing various products from them, it is necessary to take into 


account their physico-chemical properties and electrochemical reactions. 
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Plastics cannot dissolve like metals, but they are damaged or destroyed due to swelling, that is, the so-called 
embrittlement occurs, in which they soften, harden and change colour, which leads to their mechanical properties 
deterioration. 

There are several mechanisms of polymer materials destruction: 

— structure swelling; 

— dissolution; 

— bonds are destroyed under heat, chemical reactions or free radical reactions. The surface layer embrittlement can 
occur by any one mechanism or by a combination of mechanisms [1 1—13]. 

Dissolution and swelling without breaking chemical bonds are the main causes of embrittlement upon contact with 
impurities, and this process is complex, since small amounts of polymer solution penetrate inside [14], forming a 


seeping layer with altered physical properties (Fig. 1). 


Solvent 


Liquid layer 


Gel layer 


Solvent 


Solid expanded layer 
Polymer 


Permeable layer 


Polymer 


Fig. 1. Layer-structure of embrittlement process 


The research objective is to study the polyethylene resistance to oil and associated water at room temperature, both 


at rest and in motion. 


Materials and research methods. Table | shows the substances dissolved in oil after its extraction. 


Table | 
Substances soluble and insoluble in oil after its extraction 
Soluble substances Insoluble substances 

solid substances gases metal substances organic substances 

positive electrolytes: plants and O, HDS, Fe(OR)3, Compounds obtained 
Mg", Fe*, K*, Na*, Mn* animals N2,CO;z FeS), as a result of the 
negative electrolytes: SCO? addition of 

Br, Cl, F, CO;°, sand with clay substance, demulsifiers, bacteria, 
SO,°, HCO; sulfate compounds, aquatic plants as well 

carbonates, alkaline land- as plant and animal 

based minerals. residues. 


Table 2 shows the composition of the associated water used in the experiments. 
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Composition of associated water used in experiments 


Table 2 


Analysis Measured value Measurement units 
Density 1.015 gr/cem* 
Cr 10,100 
Na* 2,427 
Mg", 370 
Ca’* 1,270 
Total stiffness 4,243 ppm 
CaCO; 
Alkalinity 530 
(CaCO3)m 
SO,* 1,093 
NO; 100 
NH4 13 
AS 130 
PH 8.5 


Petroleum alcohol was used as a cleaning agent. The experiments were carried out using granules and pipes made of 


high-pressure polyethylene HDPE of low density (Table 3). 


Table 3 


Technical characteristics of polyethylene granules and pipes used in the experiments 


Granules Pipes 
Type Low-density high-pressure polyethylene 
Density 0.914 0.917 
Shore hardness 44.6 46.7 
Elongation-to-cut ratio absent 4.42 
The degree of softening onset 109-122 110-112 


Laboratory instruments. Oil and associated water passed through polyethylene 


assembled by the author (Fig. 2). The complex of equipment included: 


— sealed glass containers; 


— scale with a sensitivity of 4 


t+ 0.1 g,%; 


— electronic micrometer with an accuracy of 1 micrometer; 


electronic microscope conne 


cted to a computer 


— An electron microscope connected to a computer. 


pipes in a laboratory device 
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Fig. 2. Laboratory device for passing oil and associated water in polyethylene pipes 


Experiments. In calm mode (immersion method), the sequence of actions when setting up the experiment was as 
follows: 


— weight W, and volume V, were determined for several polyethylene granules samples, the characteristics of which 
are given in Table 3; 

— the samples were completely immersed in oil, associated and distilled water in the sealed containers for certain 
periods of time; 

— the taken samples were cleaned; 

— the sample weight was determined as W2 and its volume as V> during the experiment; 

— the weight change was calculated as AW=W, — W, and the volume change — as AV = V2 -V, 

—the relative weight change as AW 460 and the relative volume change as oY 00 were calculated after each 

1 1 

sample testing. 

In dynamic mode: 

— oil and associated water flowed in the polyethylene pipes at a constant speed at the ambient temperature, their 
weight as W, and the wall thickness were determined for specific periods of time; 

— the pipes were cleaned after each experiment; 

— their weight as W, was determined after the experiment W>, and the wall thickness — after the experiment as X>; 

— the weight change was determined as AW = W,— W,, the volume change was as AV = V2 — V,, the relative weight 
change was as a -100 and thickness change was as AX after the experiment. 

\ 

Results and discussion 

The weight and volume changes test results of particles immersed in oil, associated and distilled water were 
summarized in Tables (4—9). A graphical dependence was found for the relative weight and size changes of oil granules 


immersed in water, associated and distilled water (Fig. 3 and 4). 
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Table 4 
Weight change of granules immersed in oil 
No. 1 2 3 
Immersion W), WwW AW, W), W AW, W), W AW, 
time, days = a % - er % gr x % 
10 27.022 27.673 2.240 26.9611 27.5542 | 2.200 25.384 | 25.921 | 2.261 
30 26.652 | 28.6376 | 7.450 26.4489 | 28.4589 | 7.600 26.577 | 28.613 | 7.661 
60 25.839 | 28.2497 | 9.330 23.486 25.7346 | 9.574 25.266 | 25.266 | 9.396 
90 26.863 | 29.6979 10.55 27.3432 | 30.1732 | 10.350 | 25.683 | 25.683 10.28 
Table 5 
Volume change of oil granules 
No. 1 2 3 
Immersion time, Vi, V>, AV, Vi, V>, AV, Vi V>, AV, 
days a ar %o ue sal %o ae sae %o 
10 50.0 51.0 2.0 50.0 51.0 2.0 50.0 51.0 2.0 
30 50.0 52.1 4.2 50.0 52.4 4.8 50.0 52.3 4.6 
60 50.0 53.0 6.0 50.0 53.0 6.0 50.0 53.0 6.0 
90 50.0 53.7 74 50.0 53.7 TA 50.0 53.9 78 
Table 6 
Weight change of granules immersed in the associated water 
No. 1 2 3 
Immersion W), Wr, AW, Wi, Wo, AW, Wi, Wo, AW, 
time, days a & % ee oe % e ee % 
10 55.045 55.430 0.70 58.765 58.765 0.70 28.514 28.695 0.63 
30 55.333 55.948 1.11 26.400 26.400 1.10 26.962 27.25 1.07 
60 45.831 55.629 1.45 27.332 27.332 1.60 26.614 26.99 1.44 
90 26.226 | 26.6452 | 1.595 27.0853 | 27.0853 | 1.650 | 26.8433 27.254 1.530 
Table 7 
Volume change of particles immersed in the associated water 
No. 1 2 3 
Immersion Vi, V>, AV, Vi V>, AV, Vi V>, AV, 
time, days sae “ee Yo i ea %o ae a %o 
10 50.0 50.0 0.0 50.0 50.0 0.0 50.0 50.0 0.0 
30 50.0 50.5 1.0 50.0 50.5 1.0 50.0 50.5 1.0 
60 50.0 51.0 2.0 50.0 51.0 2.0 50.0 51.0 2.0 
90 50.0 51.52 2.5 50.0 51.1 2.2 50.0 51.1 2.2 
Table 8 
Weight change of granules immersed in distilled water 
No. 1 2 3 
Immersion Wi, Wo, AW, Wi, W2 AW, Wi, W2 AW, 
time, days a a % x x % gr gr % 
10 56.003 56.003 | 0.000 26.761 26.762 0.00 26.723 26.723 0.00 
30 54.40 54.40 0.000 27.058 27.058 0.12 26.599 26.599 0.00 
60 55.48 55.489 | 0.000 26.303 26.303 0.00 26.084 26.084 0.00 
90 55.37 55.37 0.000 26.8300 | 26.8300 | 0.000 | 27.1665 | 27.1665 | 0.000 
Table 9 


Volume change of granules immersed in distilled water 
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No. 1 2 3 
Immersion time, Vi, V2, AV, Vi, V>, AV, Vi V>, AV, 
days om on % om om % om ont % 
10 50.0 50.0 0.0 50.0 50.0 0.0 50.0 50.0 0.0 
30 50.0 50.1 0.2 50.0 50.1 0.2 50.0 50.1 0.2 
60 50.0 50.0 0.0 50.0 50.0 0.0 50.0 50.0 0.0 
90 50.0 50.0 0.0 50.0 50.0 0.0 50.0 50.0 0.0 


=@— Distilled water 
—f Associated water 
—* Oil 


Relative weight, % 


0 20 40 60 80 100 
Time, day 


Fig. 3. Relative weight change of the immersed granules by time 
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Fig. 4. Relative volume change of immersed granules by time 


Figures 3 and 4 show that the time dependence of the relative weight and size changes of polyethylene granules has 
the character of a quadratic dependence and follows Fick's law in diffusion in the cases of three studied liquids (distilled 
water, associated water and oil), which indicates that the diffusion speed of liquid through polyethylene is the most 
important factor. 

The relative weight change of granules immersed in the distilled water is not mentioned, but it is rarely noticeable in 
the associated water and increases more in the oil, as shown in Figure 4. This is due to the fact that salt diffusion 
significantly effects on the relative weight and volume between the granules immersed in the distilled water, on the one 
hand, and the granules immersed in the oil and related water, on the other hand. 

The test results of the pipe weight and thickness changes, in which oil and associated water flow, are systematized in 
Tables 10 and 11. A graphical dependence of the pipe relative weight and thickness changes, through which the oil 
passes with the associated water flow shown in Figure 5 and 6, was obtained. 

The sample cross section micrographs were taken before and after the tests with a 510-time increase and shown in 
Figure (7). 

Table 10 
Pipe weight and thickness changes through which oil flows 
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Immersion time, days 30 60 90 
Pipe W, Xi, W, AX AW, AX W, AX 
number (er) mm % %o %o 
1 18.019 2.483 0.60 0.024 1.07 0.05 3.02 0.17 
2 18.6440 2.485 0.90 0.032 1.12 0.07 3.01 0.19 
3 18.6009 2.490 0.58 0.022 1.14 0.07 3.12 0.21 
4 18.6441 2.483 0.91 0.032 1.12 0.096 3.00 0.19 
5 18.5999 2.487 0.37 0.022 1.13 0.08 3.11 0.22 
Table 11 
Pipe weight and thickness changes through which the associated water flows 
Immersion time, days 30 60 90 
Pipe W, Xi, W, AX W, AX W, AX 
number (er) mm % % %o 
1 15.22 2.483 0.399 0.012 0.521 0.025 1.35 0.087 
2 19.134 2.484 0.404 0.015 0.478 0.029 1.42 0.086 
3 15.905 2.483 0.409 0.011 0.459 0.028 1.38 0.075 
4 18.6009 2.400 0.39 0.013 0.469 0.027 1.40 0.080 
5 18.019 2.483 0.38 0.012 0.498 0.027 1.39 0.079 
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Fig. 5. Pipe relative weight change through which oil and associated water flow, over time 
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Fig. 6. Pipe relative thickness change in which oil and associated water flow, depending on time period 
Due to Figures (6) and Table (12), it is noted that the relative change in the weight and thickness of the polyethylene 
pipe walls through which oil passes is greater than those through which the associated water passes, which indicates that 


polyethylene is more affected by oil. 
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8) h) 
Fig. 7. Sample cross section micrographs before and after tests with a 510-time increase: a — cross section of a clean granule; 


b — cross section of an oil granule; c — cross section of a pipe with associated water in oil; d — cross section of a pipe immersed in 
oil; e — inner surface of the pipe with associated water; g — cross section of the pipe immersed in oil after 200 hours; 
h — inner surface of the pipe immersed in oil 


Figure 7 shows cuts of pipe and granules samples before and after the tests, reflecting the surface colour and 
smoothness of the cross sections subjected to the tests, compared to the clean cross sections, which confirms the 


previous results obtained due to the basis of the graphs. 


Conclusions 

The experiments and the results obtained have shown that the liquid diffusion speed through polyethylene is the 
most significant factor. Immersion in oil has a greater impact than immersion in associated water and distilled water 
according to the presence of salts. 

It has been found out that the relative change in the weight and thickness of the polyethylene pipe walls through 
which oil passes is greater than those through which the associated water passes, and the microscopic cross sections 


images in the samples before and after the tests have confirmed the obtained results. 
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Abstract 

Introduction. Electronic devices capable of collecting individual telemetry data have opened up prospects for 
preclinical detection of COVID-19 signs. Known solutions involve the analysis of information that is difficult to obtain 
at the moment. We are talking, specifically, about the blood condition or a PCR test. This significantly limits the 
possibility of integrating algorithms with wrist gadgets. At the same time, the cardiovascular system as an object of 
observation is quite informative, the data collection is well developed. The article describes the problem of detecting 
covid anomalies in rhythm strips. The work aims at creating a mathematical model based on machine learning 
algorithms to automate the process of detecting covid abnormalities in the heart rhythm. The possibility of integrating 
the results obtained with fitness bracelets and smart watches is shown. 

Materials and Methods. The work involved an open technology stack: Python, Scikit-learn, Lightgbm. When assessing 
the quality of models for binary classification, metric F'; was used. 229 cardiac rhythm strips (cardiointervalographies) 
of patients with COVID-19 were studied. The presence or absence of signs of an anomaly was determined taking into 
account the time of the rhythm strip and the intervals between heartbeats. Deviations that could indicate infection were 
shown graphically. Based on the exploratory analysis results, a list of signs indicating an anomaly was made. 

Results. As a result of the work done, a mathematical model was obtained that detected heart rate abnormalities specific 
to COVID-19 with an accuracy of 83 %. The basic features determining the predictive ability of the model were 
identified and ranked. They included the current value of the interval between heartbeats, the derivatives at the 
subsequent and previous points of measuring the duration of the heartbeat, the first derivative at the current point, and 
the deviation of the current value of the duration of the RR-interval from the median. The first indicator in this list was 
recognized as the most significant, the last — the least. For machine learning purposes, the potential of five algorithms 
was evaluated: IsolationForest, LGBMClassifier, RandomForestClassifier, ExtraTreesClassifier, SGDOneClassS VM. 
The normal and abnormal results of observations in isolation trees were visualized. A parameter was set that 
corresponded to the probability of regular observation outside the norm, and its value was selected — 0.11. Taking into 
account this indicator, a graph was constructed for the SGDOneClassSVM model. Based on the data set, using the 
cross-validation technique, the quality metric was calculated. The case in hand was a rhythm strip with a time series of 
observations taken in one continuous time interval from one person. A step-by-step process of obtaining averaged 
metric values for each model was described. In comparison, the highest indicator was recorded for the LGBMClassifier 
model, the lowest — for SGDOneClassSVM and IsolationForest. 

Discussion and Conclusions. The resulting mathematical model takes up little space in the memory of a mobile device, 
ie., it does not impose significant requirements on computing resources. The solution has an acceptable detection 
quality for pre clinical screening of COVID-19-related cardiovascular disorders. The algorithm detects anomalies in 
83 % of cases. Four minutes is enough to record a rhythm strip. The proposed scenario for using an integrated solution 
is concise and easy to implement. Widespread use of the development can contribute to the detection of COVID-19 at 
an early stage. 
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AHHOTalna 

Beedenue. QiexTpoHuHble yctTpolictBa, cimocoOuble coOupaTb aHHble MO TeseMeTpHH MHAMBUyyMa, OTKPbIIN 
TI€pCMeKTHBbI JOKIMHM4eECKOTO BbIABICHHA IWpH3HaKoB COVID-19. Vs3Bectubie pelieHua peanomararoT aHasn3 
MHOpMallHu, KOTOPYt0 CJIO2%HO NOJYYNTb B MOMEeHTe. Peusb ueT, HallpuMep, 0 cocTosHuN KpoBu uu IIL[P-tecte. Sto 
CYIeCTBeCHHO OrpaHH4nBaeT BO3MOXHOCTH MHTerpallMu asITOpUTMOB C HapydHbiIMu raspKeTamu. IIpu sToM cepye4Ho- 
cocyqHcTad cucTemMa Kak OObeKT HaOIOAeCHHA TOCTaTOUHO HH(POPMaTHBHA, CbeM JaHHBbIX Xopouio mpopabotau. B 
cTaTbe OlMcaHa 3ajjaua J[eTeKIMH KOBMHBIX aHOMAaJIMM B pHTMorpammax. Ilenb padoTbl — co3qaHHe MaTeMaTH4eckKol 
MOJleJIM Ha 6a3e arOpHTMOB MalUMHHOrO OOy4eHHA JIA ABTOMATH3AIHH TIpolecca BbIABIICHHA KOBHAHBIX aHOMAJINE B 
putTme cepgua. Hoka3aHa BO3MO%KHOCTb MHTerpallHu MOJYYCHHIX pe3yIbTaTOB c (uTHecc-OpacueTaMH WM YMHbIMU 
yqacaMu. 

Mamepuanoi u memoovi. B padote 3aqevicTBOBaIM OTKPbITbIM CcTeK TexHONOru: Python, Scikit-learn, Lightgbm. [pu 
olleHKe KayecTBa MOJeNei Td OuHapHOM Kiaccuukallun UcCOMb30BaIach MeTpuka F|. U3yaensr 229 purmorpamMM 
cepmia (KapqMouHTepBalorpadui) mamuentop c COVID-19. Hanwume unm otcyTcTBue pH3HaKOB aHOMaJIHu 
olpeyesaocb C yYeTOM BPpeMeHH PHTMOrpaMMBI HM HHTepBaIOB MexKTy cepaieOueHuamu. Ipadwuecku noKka3aHbl 
OTKJIOHCHHA, KOTOPhIe MOTYT CBU JeTeIbCTBOBATb O 3apaxKeHHH. Ilo uToraM pa3BeO4HOTO aHaiHM3a coOpaH MepeyeHb 
IIPH3HAKOB, YKa3bIBarOIlHx Ha AHOMAJIMIO. 

Pezyivmamol ucciedosanua. B pe3yibTate WposewaHHo padoTEI NWollyuyeHa MaTeMaTHYecKad MOJesIb, KOTOpat 
WeTeKTupyet cnemwu@uyupie pzua COVID-19 aHomManuu cepfeaHoro puTMa c TouHocTbIO 83%. BplaBieHbl u 
PaHKUpPOBaHb! OCHOBHEIe IIPH3HAaKH, OlpeseAIOMHe MPOrHOCTHYeCKyt0 CIHOCOOHOCTh MOA eM. ITO TeKylilee 3HadeHHe 
WHTepBalla MeKLY yapaMu cepa, Mpov3BOAHbIe B TocHenywmeit U Mpesbiqyujei TouKax M3MepeHHA 
MIpOAOMKUTeIbHOCTH cepAlleOveHHA, TepBad MNpOM3BOTHAaA B TeKyleH TOUKe HU OTKIOHEHMe OT Me{HaHbI TeKyllero 
3HayeHHA YMTembHOcTH RR-vunTepBasa. IlepBeli oKa3aTeib B 9TOM TepedHe pH3HaH HavOoee 3HaYHMBIM, 
TocueqHuii — HaMMeHee. Jia Weneli MallIMHHOrTO OOyyeHHA OeHHBaICA TOTeHWMal MATH asIrOpHTMOB: 
IsolationForest, LGBMClassifier, RandomForestClassifier, ExtraTreesClassifier, SGDOneClassS VM. 
Bu3yalv3HpoBaHbl HOPMAaJIbHbIc HW AHOMAJIBHbIC pe3yiIbTaTbI HaOMFONCHHM B W30JMpyIOUIMx WepeBbax. YcTaHOBJIeH 
HapaMeTp, KOTOPHIM COOTBETCTBYeT BCPOATHOCTH peryAPHOTO HaOIOAeHHA 3a MpeyeaMH HOPMBI, HW BbIOpaHO ero 
3HaueHHe — 0,11. C yueTom waHHoro oKa3aTesa WOcTpoeH rpaduk wa Mogemu SGDOneClassSVM. Ho Habopy 
TaHHbIx C IpHMeHeHHeM TeXHHKH TepekpecTHoH IpoBepKU paccuuTaHa MeTpHkKa KayecTBa. Peub ueT O pPHTMOrpamMMe 
C BPCMeHHBIM pAOM HaOJIOAeCHHH, CHATHIX 3a OHH HelipepbiBHbIi MHTepBasl BPeMeHU y OFHOrO YesOBeKa. OnucaH 
MOWAaroBbI Mmpowecc MowyYeHuA ycpeHeHHBIX 3HadYeHu MeTPHKU WIA KaxoOM Moyen. IIpu cpaBHenuu cambiii 
BbICOKHH WOKa3aTeIb 3a=~uKcHpoBaH y MogemH LGBMClassifier, Haumenbumme — y SGDOneClassSVM u 
IsolationForest. 

O@6cystcoenue u 3akro4yvenua. TlonyaeHHad MaTeMaTHYeCcKad MOJ[elIb 3AHMMaeCT MaJIO MeCTa B TlaMATH MOOMJIBHOTO 
yCTpolicTBa, TO eCTb He MIpebABACT 3HAYMMBIX TpeOOBaHH K BEIYHCIMTeIbHbIM pecypcaM. PerieHwe oOmaqaeT 
TIPHeMJICMbIM KaYeCTBOM J[eTeKIMU JIA MOKIMHMYeCKOTO CKPHHHHTa cBA3aHHbIX C COVID-19 cepyeuno-cocyqucTBIXx 


Information Technology, Computer Science, and Management 


67 


http://vestnik-donstu.ru 


Advanced Engineering Research (Rostov-on-Don). 2023;23(1):66—75. eISSN 2687-1653 


Hapymenuii. AnropHTM oOHapy2KHBaeT aHOMaJINH B 83 % crry4daes. [Ia 3almlMCH PHTMOrpaMMBI JOCTaTOUHO 4 MHHYT. 
IIpeqmaraempli cleHapuii MCHONb30BaHHA MHTerPHpOBaHHOTO pellieHHA JIaKOHMYeH HM JlerKo peamu3yem. LIupoxoe 
MCIOIb30BaHHEe pa3xpaOOTKH MO2XKeT CIOCOOCTBOBATS BbIABICHHIO COVID-19 Ha paHHeli cTaqun. 


Ksoyespre copa: COVID-19, mpw4unHbl cMepTH KOBHJ-MOOXKMTCIbHBIX WalMeHTOB, OCJIOXKHeEHHA B padoTe 
cepeuHo-cocyyuctoli cuctembl, III[P-rect, 1oKmmHMueCKH KOHTPOsIb Cepe4HO-COCyAMCTOM CHCTeMbI, BCTPOeCHHbIe 
WaTYHKH UYacTOTbI Tybca, puTMorpammMa, RR-nuTepBas, 9eKTpOKapqMorpaMMa cepa, aHOMaJIbHOe 110 
MIpOAOMMKUTeIBHOCTH cepyeOueHHe, cepauleOveHHe C aHOMAJIbHbIM PHTMOM, MallIMHHOe OOyyeHHe, aJIropHT™M 
LGBMClassifier. 


Baarogapnoctu. ABTopbl BbIpaxaloT OaroqapHOCTb PYyKOBOJCTBY HM MOsepaTopaM OTKpPbITOrO BCepocculicKoro 
copeBHOBaHHaA TpodeccnouasoB B ctbepe WHMppoBolt 9KoHOMUKH «LInpoBol NpopLlB» 3a MpeqocTaBJIeHHbIe TaHHbIe 
IIA UCCIeOBaHHA. 


Aaa wuTupospannsa. Mexos M.C., Ko3suyun B.O., Kayep 10.]]. Mogens MauimHHoro oOyyeHna Aid OOHapyxKeHHA 
COVID-19 ua paHHei craquu no aHoMasMaM B puTMe cepyua. Advanced Engineering Research (Rostov-on-Don). 
2023;23(1):66—75. https://doi.org/10.23947/2687-1653-2023-23-1-66-75, 


Introduction. Investigation of the impact of COVID-19 on humans remains a challenge. Thus, in 2021—2022, more 
than 16,000 scientific papers were published on this topic. One of the main causes of death of covid-positive patients 
was complications in the cardiovascular system (hereinafter referred to as CVS) caused by exposure to coronavirus [1]. 
Two methods are mainly used for preclinical diagnosis of COVID-19: biochemical method based on polymerase chain 
reaction (PCR test) and blood analysis. Contacts with medical staff needed in this case (including visits to medical 
institutions) complicate regular operational control and increase the burden on the healthcare system. Thus, it seems 
relevant to use modern technologies of preclinical control of CVS for early detection of COVID-19 signs. 

Wearable electronic devices can provide regular monitoring. The most common of them are fitness bracelets and 
smart watches with built-in heart rate sensors and the ability to perform measurements with high discreteness [2]. This 
approach opens up opportunities for analyzing data flows based on machine learning’ [3]. 

The presented study aims at creating a trainable model capable of detecting covid anomalies based only on data on 
heart rhythm. A number of papers [4-6] consider similar problems, but the solutions are based on additional 
information about the state of the blood and other characteristics”. This significantly limits the possibilities of their 
integration with wearable devices, because at the moment, it is impossible to enter the results of a blood test or a smear 
for a PCR test into the model. The novelty of the proposed solution is in the fact that only heart rate data is used, which 
can be taken with a high frequency in a way convenient for a person and interpret the indicators in real time. 

Materials and Methods 

Data characteristics. 229 impersonal rhythm strips (cardiointervalographies) of patients with COVID-19 were used 
in the research. The information was obtained in 2021 as part of the open All-Russian competition “Digital 
Breakthrough” for professionals in the digital economy. A data fragment is presented in Table 1. 


Table | 
A fragment of the data set 
Number of Time RR interval between heartbeats Sign of covid 
rhythm strip in milliseconds in milliseconds anomaly* 
81 0 576 0 
81 568 568 0 
81 1,140 572 0 
176 44,332 568 0 
176 44,968 636 1 
176 45,596 628 0 
*0 — no anomaly, | - there is an anomaly. 


Figure 1 shows the relationship of the rhythm strip (RR interval) and the electrocardiogram of the heart (ECG). 


' Permyakov SA, et al. Ehndogennye anomalii kardioritma u patsientov s COVID-19. In: Proc. VII All-Russian Conf. “Nelineinaya dinamika v 
kognitivnykh issledovaniyakh — 2021”. Nizhny Novgorod: Institute of Applied Physics of RAS; 2021. P. 109-110. (In Russ.) 

> Diagnosis of COVID-19 and Its Clinical Spectrum. Kaggle Inc. URL: https://www.kaggle.com/datasets/einsteindata4u/covidl19 (accessed: 
10.09.2022). 
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Fig. 1. Comparison of electrocardiogram and heart rhythm strip: horizontal axis shows the time in seconds, 
vertical axis for the ECG — microvolts 


In all rhythm strips from this set, there are marked abnormal areas. In Figure 2, abnormal areas are highlighted with 
a red dotted line. The x-axis shows the duration of one measurement of the rhythm strip in milliseconds, the y-axis — 
the interval between adjacent heartbeats in milliseconds. 
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Fig. 2. Chart of rhythm strip No. 69: abnormal areas are highlighted with red dotted lines, 
abnormal points — with black bullet points 


Each rhythm strip is presented by its own identifier. The duration of rhythm strips in the studied data set is different: 
4 minutes on average, 31 minutes maximum. Each measurement inside one rhythm strip has a timestamp in 
milliseconds from the start of recording. The duration of the RR interval is also presented in milliseconds. Each specific 
value in the rhythm strip allows us to talk about the signs of an anomaly (0 — no, | — there is). 2.53 % of observations 
are marked with number 1. Thus, the data set has a strong class imbalance, which is typical for anomaly detection tasks. 

In the data markup, there are various approaches to the allocation of abnormal areas. Groups of points in the vicinity 
of a characteristic peak and fall in the duration of the heart rhythm were distinguished as abnormal: 3rd, 4th, 6th 
measurements (Fig. 2). Not always the number of points in the neighborhood is marked the same — there may be a 
different number of abnormal points to the left and right of the peak. Moreover, rhythm strips with noisy indications 
were detected. This was the case when the connection with the gadget was lost, and measurements were taken when 
installing or removing the device. Sixteen rhythm strips with incorrect data had to be excluded from consideration, and 
the markup was redone: 

— only one point stands out in the anomalous section, characterizing the anomalous fragment; 

— abnormal points are indicated by black bullets (Fig. 2). 

Feature extraction. In its pure form, only one signal was presented — the value of intervals between heartbeats. 
Therefore, to refine the model, additional features were prepared based on the available signal: deviation from the 
median value and derivatives in neighboring rhythm measurements. This list of features was selected after an 
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exploratory data analysis and visual identification of the pattern in places corresponding to abnormal areas. In Figure 2, 
they were marked with a red dotted line. 

Research Results 

Metric for evaluating the quality of anomaly detection. To assess the quality of the model in the binary 
classification problem, due to the imbalance of classes, metric F, [7] (1) was used. It provided evaluating how well a 
constructed model detected a rare class. In that context, a rare class referred to abnormal heartbeats in duration — 


heartbeats with an abnormal rhythm: 
accuracyx completeness 


F, =2x Geen erica” (1) 
Here: 

— accuracy — the proportion of abnormal heartbeats correctly detected by the model from the total number of 
heartbeats that the model identified as abnormal; 

— completeness (or in other words, sensitivity) — the proportion of heartbeats that the model correctly detected as 
abnormal from the total number of abnormal heartbeats in the entire data set. 

Machine learning algorithms. As part of the study, five machine learning algorithms described below were 
applied. 

1. IsolationForest — an algorithm with uncontrolled self-learning based on extremely randomized decision trees [8]. 

2. Light Gradient Boosting Machine Classifier (LGBMClassifier) — an algorithm for gradient boosting over 
decision trees [9]. To increase the operation speed, two techniques were used: Gradient-based One-Side Sampling and 
Exclusive Feature Bundling’. 

3. RandomForestClassifier is based on decision trees and implements multiple selection of a random subset of 
features. They are used to build simpler estimators — decision trees. The results are aggregated to obtain a final 
prediction [10]. 

4. ExtraTreesClassifier is similar to RandomForestClassifier, however, it additionally implements a random 
selection of the boundary along which nodes branch in decision trees [11]. 

5. SGDOneClassSVM* — a linear version of One-Class Support Vector Machine using Stochastic Gradient 
Descent. 

IsolationForest and SGDOneClassSVM were chosen due to their wide use in anomaly detection tasks [12, 13]. 
LGBMClassifier, RandomForestClassifier and ExtraTreesClassifier perform well enough in different tasks, therefore, 
they were also used to compare the results [12, 13]. 

The specific feature of the IsolationForest and SGDOneClassSVM algorithms is that they do not require a clear 
marking of anomalous observations at the input, while it is mandatory for the rest of the algorithms used in the study. 

IsolationForest is based on the assumption that when constructing isolating trees, abnormal observations can be 
isolated (separated) in fewer operations than normal observation instances. For each observation, the algorithm 


calculates the anomaly score by the formula: 
_E(a@)) 
s(xn)=2 @™, (2) 


where h(x) — number of edges up to instance x in each isolating decision tree; FE (h(x)) — average value h(x) on the 


entire set of isolating trees; c(n) — normalizing constant for a data set of size n (3). 
c(n) = 2H(m- 1) ~ 3) 
H(k) = In(k) + ¥. (4) 
In equation (4) y — Euler’s constant equal to 0.57721... 
If observation x has an anomaly estimation value s, close to 1, then it is considered anomalous. If s is close to 0.5, 
then the observation has no obvious signs of an anomaly. If s is close to 0, then the observation can be considered 


normal (Fig. 3). 


> LightGBM: A Highly Efficient Gradient Boosting Decision Tree. www.microsoft.com URL: https:/Avww.microsoft.com/en-us/research/wp- 
content/uploads/2017/1 1/lightgbm.pdf (accessed: 10.09.2022). 
“ Online One-Class SVM. Scikit-learn developers (BSD License). scikit-learn.org URL: https://scikit-learn.org/stable/modules/sgd.html#online-one- 


class-svm (accessed: 10.09.2022). 
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Fig. 3. Normal and anomalous observations in isolating trees (the authors’ figure) 


SGDOneClassSVM is based on the opposite approach to IsolationForest. The algorithm determines the boundaries 
of normal observations and compares all new observations to the boundaries of this norm to identify an anomaly. 


Feature Significance. An assessment of the degree of impact of features on the predictive ability of the model is 


- 


shown in Figure 4. 


Q next_diff —| 

: So 
S 

1 

o 

oy 


median_deviation -} 


200 400 600 800 1,000 


Feature significance, ea 


Fig. 4. Feature significance diagram: x — current value of the interval; next_diff — derivative at the next point of measuring the 
heartbeat duration; prev_diff — derivative at the previous point of measuring the heartbeat duration; diff — the first derivative at the 
current point; median_deviation — deviation of the current value of the RR interval duration from the median 

To calculate the numerical significance estimate, a mechanism built into LGBMClassifier was used, which returns 
an array of numerical estimates for each feature via the feature_importances_ property of the trained model. 
Significance in models based on gradient boosting over decision trees is usually calculated on the Gini-impurity 
Index? [14], used in the process of determining the branching points when training the model: 

Gini(d) = 1- YEP}. (5) 
Here, d — a set of observations that match the conditions at the considered branching point, d € D; k — number of 
classes presented in the entire training dataset D; p; — probability of observations belonging to class i at the considered 
branching point of the decision tree. 

The following features were the most significant: the current value of interval (x), the derivative at the next 
(next_diff) and previous (prev_diff) points of measuring the heartbeat duration (Fig. 4). A complete list of the features 


used is given in Table 2. 


° Karabiber F. Gini Impurity. learndatasci.com URL: https://www.learndatasci.com/glossary/gini-impurity/ (accessed: 10.09.2022). 
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Table 2 
List of features used 
No. Feature Description 

1 x RR interval at the current measuring point 
2 next_diff First derivative at the next point 
3 prev_diff First derivative at the previous point 
4 diff First derivative at the current point 
Be - aeatandemahoa pos of ne current value of ae RR interval duration from 

the median within one rhythm strip 


Comparison of models. For effectiveness of SGDOneClassSVMmodel, it is important to select parameter nu, 
which corresponds to the probability of detecting regular observation outside the norm. In other words, nu determines 
the upper bound of the error rate when training the model, and the lower bound of the support vector fraction®. To 
select nu taking into account the available data nature, the quality metric was additionally assessed at different values of 


the specified parameter (Fig. 5). As a result, nu equal to 0.11 was selected. 
F,, units 
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Fig. 5. Estimate of parameter nu (on the horizontal axis) for SGDOneClassSVM model. 
On the vertical axis — values of metric F, 


To calculate the quality metric on various models, the entire data set was used through the cross-validation 
technique. Within one rhythm strip, we had a time series of observations taken in one continuous period of time from 
one person, therefore, they should be considered as dependent [15]. The following strategy was used to divide the data 
into training and test sets. The selected data set consisted of 213 rhythm strips marked with a unique identifier (id). This 
made it possible to allocate rhythm strips for training and testing models. A set of rhythm strips for the test could be 
randomly selected by identifiers. The approach used in the presented work is described below. 

I. Five actions were performed in the data partitioning cycle. 

1. The initial number for generating pseudo-random numbers was fixed (seed) —np.random.seed(fold), where fold 
— number of the current data partition. 

2.42 random integer values were generated in the range from 1 to 213. This was how we got random numbers of 
rhythm strip identifiers for the test data set. 

3. The numbers of rhythm strip identifiers that remained after the selection of identifiers for the test were entered in 
a separate list. They were used for a training set. 

4. Models were trained on rhythm strips from the training set, and prediction quality metrics were evaluated on 
rhythm strips from the test set. 

5. The quality metric value was recorded for each model calculated on the test set of rhythm strips at the current data 


split. 


° SGDOneClassS VM documentation. Scikit-learn developers (BSD License). scikit-learn.org URL: https://scikit- 
learn.org/stable/modules/generated/sklearn.linear_model.SGDOneClassS VM.html#sklearn.linear_model.SGDOneClassS VM (accessed: 10.09.2022). 
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II. Steps 1-5 were repeated for each data split number. 
Il. The obtained values of the quality metric were averaged for each of the models. 


A comparative assessment of the average prediction quality metric for each model is given in Table 3. 


Table 3 
Evaluation of quality metric Fy 
Model Metric F,* 
LGBMClassifier 0.8328 
RandomForestClassifier 0.7638 
ExtraTreesClassifier 0.7369 
SGDOneClassSVM 0.0169 
IsolationForest < le-4 
* Average value for the selected cross-validation 
strategy on five partitions. 


Discussion and Conclusions. A mathematical model for detecting anomalies in the heart rhythm with an accuracy 
of 83 % has been developed. According to quality metric F), the model based on LGBMClassifier algorithm turned out 
to be the best. IsolationForest and SGDOneClassS VM showed weak results on current data. 

The proposed model can be implemented as a component of the software part of wearable personal smart devices. 
The proposed scenario for using the solution is as follows: 

— the recording of the rhythm strip is activated on a personal wearable device through the user interface; 

— upon completion, the record is submitted to the developed model for analysis; 

— based on the results of data analysis, the mathematical model issues a notification about the presence or absence of 
anomalies on the screen of the wearable device. 

Note that an average of 4 minutes is probably enough to record one rhythm strip. During this time, it is possible to 
detect covid anomalies in the heart rhythm. 

The model occupies 493 kilobytes in the memory of the wearable device, which is quite suitable for practical use. 
The solution relies only on information about the heart rate and does not involve factors inaccessible to mobile personal 
gadgets. 

Improving the accuracy of anomaly detection involves additional research. They should focus on the development of 
unique features that are detected by the initial heart rate signal. However, the current solution already makes it possible 
to quickly and easily assess the probability of COVID-19 at an early stage. This, along with the implementation of 
medical recommendations, can further contribute to reducing the risk of mortality from the negative impact of 


coronavirus infection on the cardiovascular system. 
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Abstract 

Introduction. The banking sector assigns high priority to data storage, as it is a critical aspect of business operations. 
The volume of data in this area is steadily growing. With the increasing volume of data that needs to be stored, 
processed and analyzed, it is critically important to select a suitable data storage solution and develop the required 
architecture. The presented research is aimed at filling the gap in the existing knowledge of the data base management 
system (DBMS) suitable for the banking sector, as well as to suggest ways for a fault-tolerant data storage cluster. The 
purpose of the work is to analyze the key DBMS for analytical queries, determine the priorities of the DBMS for the 
banking sector, and develop a fault-tolerant data storage cluster. To meet the performance and scalability requirements, 
a data storage solution with a fault-tolerant architecture that meets the requirements of the banking sector has been 
proposed. 

Materials and Methods. Domain analysis allowed us to create a set of characteristics that a DBMS for analytical 
queries (OnLine Analytical processing —- OLAP) should correspond to, compare some popular DBMS OLAP, and 
offer a fault-tolerant cluster configuration written in xml, supported by the ClickHouse DBMS. Automation was done 
using Ansible Playbook. It was integrated with the Gitlab version control system and Jinja templates. Thus, rapid 
deployment of the configuration on all nodes of the cluster was achieved. 

Results. For OLAP databases, criteria were developed and several popular systems were compared. As a result, a 
reliable cluster configuration that met the requirements of analytical queries has been proposed for the banking industry. 
To increase the reliability and scalability of the DBMS, the deployment process was automated. Detailed diagrams of 
the cluster configuration were also provided. 

Discussion and Conclusions. The compiled criteria for the DBMS OLAP allowed us to determine the need for this 
solution in the organization. Comparison of popular DBMS can be used by organizations to minimize costs when 
selecting a solution. The proposed configuration of the data warehouse cluster for analytical queries in the banking 
sector will improve the reliability of the DBMS and meet the requirements for subsequent scalability. Automation of 
cluster deployment by the mechanism of templating configuration files in Ansible Playbook provides configuring a 


ready-made cluster on new servers in minutes. 
Keywords: DBMS, OLAP, data warehouse, ClickHouse, failover cluster. 
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AnHOoTanHa 

Beedenue. bankosckul ceKTop mpusaeT OombuIoe 3HAaYeHHe XpaHeHHIO TaHHBbIX, MOCKOJIbKY 9TO KPHTHYCCKH BaxKHBIM 
acliekT Ou3Hec-onepalHi. OObeM aHHbIX B WaHHOM ciepe HeyKIOHHO pacTeT. C yBesH4YeHHeM OOBeEMOB JaHHBIX, 
KOTOpble HeEOOXOAMMO XpaHUTb, OOpabaTHIBaTbh M aHaJIM3HpOBaTb, KpaliHe BaxKHO BbIOpaTbh NOAXOAAWee pelieHHe DA 
XpaHeHHa JJaHHbIX HM pa3spadoTaTb HeoOxogzuMyro apxuTeKtypy. IIpeyctaBieHHoe HMcceqoBaHve HalpaBeHo Ha TO, 
YTOOBI 3aN0JHHTb upoben B cylecTByIomjMx 3HaHHAx CYB], nogzxopsawHx Id OaHKOBCKOrO ceKTopa, a TakoKe 
TIpeAIOKUTh CHOCOObI AI OTKA30YCTOMYMBOTO KJlacTepa xpaHeHHaA aHHBIX. Llerb paOoTbl — aHasIM3 KJIIOUCBBIX 
CYB] wa anamuTwueckux 3anpocos, onmpesemenve mpHoputetos CYB] ana OaHKoBcKoro ceKTopa MH pa3pa0oTKa 
OTKa30ycTOM4MBOrO KylacTepa xpaHeHHa anHbIx. Jia BbINOMHeHHA TpeOoBaHHH K MpOM3BOAHTebHOCTH U 
MaclUTaOMpyeMOCTH TipeANOXKeHO pellieHHe JIT XpaHeHHA TaHHbIX C OTKa30yCTOMYMBOM apXHTeKTypoH, OTBeYAaTOLee 
TpeOoBaHHAM OaHKOBCKOrO CeKTopa. 

Mamepuaavi u memoooi. Ananu3 ipeqMeTHO OOacTH NO3BONMI CO31aTb Hadop XapaKTepHCTHK, KOTOPbIM JODKHAa 
cooTBeTcTBoBaTb CYbJ] aa avamuTM4eckux 3armpocos (OLAP), BsinoHHTbh cpaBHeHHe HeKOTOPBIX MOMYyApHbIX 
OLAP CYB], u mpeqnoxuTs oTKa30ycTOMuMBylO KacTepHy!O KOH®urypalMio, HallucaHHyIo Ha s3bIKe xml, 
nosepxuBaemyro CYB] ClickHouse. ABToMaTu3alHA BbIMOHeHa c HOoMoubIO Ansible Playbooks. On mHTerpupoBaH 
c CHcTeMOM yupaBleHua Bepcnamu Gitlab u miaOnoHamn Jinja. Takum oOpa30M ocTuraeTca ObICTpoe pa3sBepTEIBaHHe 
KOH(urypalluu Ha BCexX HOax KlacTepa. 

Pezyismamot uccnedosanua. na 0a3 QanHerx OLAP Osim pa3spaboTaHbI KpHTepHu, MpoBeyeH cpaBHUTeJIbHbIi 
aHaJIM3 HECKOJIBKMX MOMMyIAPHbIX cHcTeM. B pesynbTaTe Oba WpewoxwKeHa HayexKHad KlacTepHad KOHurypalHA B 
OaHKOBCKOM MHAYCTpHu, KOTOpad yOBeTBOpseT TpeOoBaHHAM aHaIMTMYeCKHX 3anpocos. Jina yBerM4eHHa 
HajjexKHOCTH HM MacuITaOupyemoctu CYB] mporecc pa3BepTpiBaHua OBI aBTOMaTH3NpoBaH. Takxke IpHBeeHbI 
WeTasIbHble CXeMbI KOHUTypalHn KslacTepa. 

O6écyorcdenue u 3aKnio“uenua. CocTaBileHHble KpuTepun ya OLAP CYB/J] nosBonsaroT ompeseuTh HeOOXO,HMOCTb 
JaHHOrO pelieHHA B OpraHH3auun. CpaBHenve nonyapHErx CYB] MoxerT OpiTb UciOb30BaHO OpraHH3alHAMH JIA 
MHHMMV3aHH 3aTpaT upH BbIOope pewienua. IIpeqnaraemaa KOHuTypalua KacTepa xpaHvIMUla WaHHbIX JA 
aHaJIMTHYeCKHX 3alIpOCoB B OaHKOBCKOl ciepe MO3BOJIMT MOBbICHTb HayexHocTs CYB] u yqoBleTBopuTs TpedoBaHHa 
K TlocleyyrouleH MaciuTaOupyemocTu. ABTOMaTH3al[HA pa3BepTbIBAHHA KjlacTepa MyTeM MexaHH3Ma WadOHH3alHu 
KOH@UTypalHOHHBIx dalioBn B Ansible Playbooks no3BomaeT HaCTpOHTb TOTOBbI KIacTep Ha HOBbIX CepBepax 3a 


MHHYTBI. 
K.unovespre copa: CYB, OLAP, xpanumnure qanHbrx, ClickHouse, orka30ycTOH4MBbIit KacTep. 


Baarogapuoctu. Astop BbipaxaeT OnarogapHocts B.A. boraTsipepy, OKTOpy TexXHW4eCKHX HayK, Mpodeccopy 
Kaepbl BbIMHCIMTeIbHOM TexHukH YuuBepcuteta MTMO, mouetHomy padoTHuky HaykKH MM TexHuKH PO, 


TIPOBOAMBINCMY SKCIICPTHbIC HHTCPBbIO COBMECCTHO C aBTOPOM CTaTbH. 
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Introduction. Data storage in the banking sector is one of the key business factors. To ensure the security of 
customer information and transactions, it is required to take measures of protection, distribution and creation of 
backups. For operational analysis, employees should be able to make operational analytical requests to the data 
warehouse, while not interfering with the work of other processes within the organization and without causing a heavy 
load on the storage itself. Databases and Data Warehouse are information systems in which data is stored, but they are 
also used to solve various tasks. The article describes what such systems do, what the main differences between them 
are, and why their effective use is essential for business development. 

Many organizations make mistakes in designing the architecture of databases and data warehouses, losing sight of 
aspects of information security, scalability and fault tolerance. The urgency of this problem is due to the intensive 
development of systems in banks, the expansion of their fields of application and the increase in the amount of data in 
need of constant analysis. For operational analysis of a large amount of data, a storage is needed that must meet all 
reliability and security requirements. 

Effective decision-making processes in business depend on high-quality information. In today's competitive 
business environment, flexible access to a data warehouse is required, organized in such a way as to increase business 
productivity, provide fast, accurate and up-to-date data understanding. The data warehouse architecture is designed to 
meet such requirements and is the basis of these processes [1—5]. 

The objective of the work is to determine the priority DBMS for performing analytical queries in the banking sector 
and design a fault-tolerant data warehouse cluster. This solution will significantly increase the speed of execution of 
analytical queries, solve problems with scalability and reliability of the data warehouse. 

Materials and Methods. The database stores real-time information about one specific part of the business. Its main 
task is to process daily transactions. Databases use Online Transaction Processing (OLTP) to quickly delete, insert, 
replace and update a large number of short online transactions. 

Data warehouse is a system that collects data from lots of different sources within an organization for reporting and 
analysis, using operational analytical processing (OLAP) to quickly analyze large amounts of data. This system focuses 
on reading, rather than changing historical data from lots of different sources, therefore, compliance with ACID 
(Atomic, Consistent, Isolated and Durable) requirements is less strict. Data warehouses perform complex functions of 
aggregation, analysis and comparison of data to support management decision-making in companies. 

A warehouse in the banking sector may contain: 

— user account information (personal data, addresses, phone numbers); 

— information about banking products and services (loans, deposits, plastic cards, mobile banking, etc.); 

— data on transactions (including card transactions) in minimal detail for the last three years; 

— information about accounts, balances on them, etc. 

To meet the needs for OLAP, there are separate types of database management systems (DBMS) [3-6]. Each of the 
systems has its own characteristics in the construction of architecture. 

To perform an effective analysis of compliance with these requirements, warehouses must: 

— have a high capacity capable of accommodating huge amounts of data (billions or trillions of rows); 

— be organized as wide tables with multiple columns; 


— perform queries with a small number of columns; 
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— have a high query execution speed (in milliseconds or seconds); 

— provide for most of the read-only requests; 

— support fast bulk data loading when updating (more than 1,000 rows at a time) and adding, but without changing 
them; 

— have high throughput to process a single request (up to billions of rows); 

— have high reliability; 

— ensure data security and consistency. 

For the OLAP scenario of work in the banking sector, it is preferable to use column-based analytical databases, 
since they can store a lot of columns in a table, which will not affect the speed of reading data. Column-based DBMS 
provide strong compression of data in columns, since data in one column of the table is usually of the same type, which 
cannot be said about a row. They also enable to get a tenfold increase in query execution speed on lower-power 
equipment. At the same time, thanks to compression, the data will occupy 5-10 times less space on the disk than in the 
case of traditional DBMS [7-11]. 

During the requirements analysis, the following column DBMS were selected: ClickHouse, Vertica, Amazon 
Redshift. 

ClickHouse is the preferred solution due to the following advantages: open source; it is possible to define some or 
all structures that will be stored only in memory; high speed; good data compression; http and command line interface; 
cluster can be scaled horizontally; high availability; ease of installation and configuration. Installation is carried out on 
the organization's servers in an isolated segment, which meets the security requirements for sensitive data in the 
banking sector. The DBMS is also included in the register of domestic software; therefore, it provides implementing 
this software product in state-owned companies. 

Amazon Redshift solution is provided only as a cloud service. For organizations from the banking sector that cannot 
place their data in the clouds for a number of security-related reasons, this product loses its appeal. 

Vertica is an alternative version of ClickHouse with a paid license for large clusters and the installability on the 
company's local servers. 

The implementation of the distributed data warehouse architecture is presented below. To increase fault tolerance 
and performance, the implementation of a distributed ClickHouse failover cluster with three shards and two replicas is 
proposed. 

Sharding (horizontal scaling) makes it possible to write and store parts of data in a distributed cluster, process and 
read them in parallel on all nodes of the cluster, increasing data throughput. 

Replication is copying data to multiple servers; thus, each bit of data can be found on multiple nodes. 

Scalability is determined by sharding or segmentation of data. The reliability of the data warehouse is determined by 
data replication [12-16]. 

Sharding and replication are completely independent, different processes are responsible for them. It is required to 
localize small data sets on one shard and ensure a fairly even distribution across different shards in the cluster. To do 
this, it is recommended to take the hash function value from a field in the table as a sharding key. 

Sharding and replication are completely independent, different processes are responsible for them. It is required to 
localize small data sets on one shard and ensure a fairly even distribution across different shards in the cluster. To do 
this, it is recommended to take the hash function value from a field in the table as a sharding key. 

Depending on the number of available resources and servers, it is proposed to implement this configuration on 3 or 6 


nodes. For a production environment, it is recommended to use a cluster of 6 nodes. It should be noted that replication 
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does not depend on sharding mechanisms and works at the level of individual tables, and also, since the replication 
coefficient is 2, each shard is represented in 2 nodes [17—19]. Configuration options are described below. 

The logical topology diagram is as follows: 

3(Shard) x 2(Replicas) = Clickhouse Cluster of 6 nodes. 
The probability of trouble-free operation of a system with 2 replicas and 3 shards on 6 nodes is equal to: 
P, =[1-G-p)’]’. 

The probability of trouble-free operation is an objective possibility that the system will work for time ¢ without 
restorations [7, 13]. 

Thus, a table containing 30 million rows will be distributed evenly across 3 nodes of the cluster. The remaining 
3 nodes will store replicas of the data. When one of the cluster nodes is disabled, data will be taken from another 


available node that contains its replica, thereby achieving reliability [20]. A cluster of 6 nodes is shown in Figure 1. 


node O01 node 02 node 03 


shard 01 shard 02 shard 02 
replica 01 replica 01 replica 02 


node 04 node 05 node 06 


shard 03 shard 03 shard 01 
replica 01 replica 02 replica 02 


Fig. 1. Fault-tolerant cluster of 6 nodes (the authors’ figure) 
To replicate data and execute distributed DDL queries, we need to use +1 node with ZooKeeper installed. You can 
also use ClickHouse Keeper, compatible with ZooKeeper, which does not require installation on a separate server. 
An example of a fragment of the configuration file is shown in Figure 2, from which it can be seen that the shard has 


replication configured for the Ist and 6th nodes. 


<yandex> 
<remote_servers> 
<cluster_1> 


<shard> 
<weight>1</weight> 
<internal_replication>true</internal_replication> 
<replica> 
<host>{{ node1 }}</host> 
<port>9000</port> 
</replica> 
<replica> 
<host>{{ node6 }}</host> 
<port>9000</port> 
</replica> 
</shard> 


Fig. 2. Fragment of the configuration file for 6 nodes (the authors’ figure) 
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An option of the cluster configuration of 3 nodes with cyclic replication is shown in Figure 3. 


node O01 node 02 node 03 


shard 01 ina shard 03 
replica 01 Rsaceeha replica 01 
replica 01 


shard 03 shard 01 shard 02 
replica 02 replica 02 replica 02 


Fig. 3. Fault-tolerant cluster of 3 nodes (the authors’ figure) 


This implementation requires two different segments located on each node. The main problem arises due to the fact 
that each shard has the same table name, ClickHouse cannot distinguish one shard/replica from another when they are 
located on the same server. 

To solve this problem, it is needed: 

— to place each shard in a separate database (schema); 

— to set default_database for each shard; 

— to set internal_replication parameter of each shard to true; 

— to use an empty database parameter in a distributed table DDL script. 

For this topology in an industrial environment, 6 server nodes are required, where each server stores data from only 
one segment, a security trapdoor for a separate database is not required. To save resources in the development or testing 
area, a configuration with 3 nodes can be used. 

Automation is performed using Ansible Playbooks and integrated with Gitlab version control system. Thus, rapid 
deployment of the configuration on all nodes of the cluster is provided. When changing the configuration, it can be 
applied to all nodes with a single command or deploy a new DBMS cluster in a few minutes [21]. 

Research Results. The fault-tolerant cluster of the analytical DBMS provides redundancy for important system 
components, which allows for continuous operation even in case of errors in individual cluster nodes. This is done 
through load balancing, data replication between cluster nodes, and high reliability of the components used in the 
cluster. The result is an increase in the availability and reliability of the analytical DBMS, which is business-critical 
when analytical queries play a key role. The fault-tolerant cluster configuration of the data warehouse for analytical 
queries in the banking sector, taking into account the automation of the deployment process, enables to increase the 
reliability of the analytical data warehouse and meet the requirements for scalability. The developed task of automating 
cluster deployment using the mechanism of templating configuration files in Ansible Playbooks provides for the 
configuration of a ready-made cluster on new servers in a few minutes. The tasks of the template include operations to 
install the required packages, create the needed configuration and launch the cluster. 

An example of configuration files for automatic deployment of a DBMS cluster is shown in Figure 4. The j2 
extension says that they are created using the Jinja template engine. Purpose-built placeholders in the template provide 
writing code similar to Python syntax. Parameters are passed to the template for automatic insertion into the final 
document, thereby achieving automatic assembly into development, testing and industrial operation zones, which does 


not require manual modification of configuration files. 
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’ templates 
clickhouse_config.xml.j2 
clickhouse_keeper.xml.j2 
clickhouse_Idap_auth.xml.j2 
clickhouse_ldap_user_directory.xml.j2 
clickhouse_macro_n1.xml.j2 


clickhouse_macro_n2.xml.j2 


clickhouse_macro_n3.xml.j2 


clickhouse_macro_n4.xml.j2 
clickhouse_macro_n5.xml.j2 
clickhouse_macro_n6.xml.j2 
clickhouse_users.xml.j2 


cluster.xml.j2 


Fig. 4. Configuration files 


Description of configuration files: 

clickhouse_config.xml.j2 — general cluster configuration; 

clickhouse_keeper.xml.j2 —- zookeeper configuration, which is responsible for node synchronization and 
replication; 

clickhouse_Idap_auth.xml.j2 — LDAP connection configuration for data security; 

clickhouse_ldap_user_directory.xml.j2 — role-based configuration by access groups to ensure data security; 

clickhouse_macro_n1(6).xml.j2 — macro files (each node has its own); 

clickhouse_users.xml.j2 — configuration file for creating local users needed for administration; 

cluster.xml.j2 — cluster configuration file. 

To test the reliability of this configuration, an experiment was conducted during which data was loaded into a 
DBMS cluster with a replication factor equal to 2. The dwh schemas and cluster_test_data tables were created on each 
of the nodes of the DBMS cluster, and a distributed table was created on the dwh cluster.cluster_test_data_distributed. 
The rows of the dwh.test_data_distributed table distributed across the cluster were 27,547,855. The rows of the 
dwh.cluster_test_data table with each of the cluster nodes are listed below: 

9,186,544 rows — Ist node; 

9,182,959 rows — 2nd node; 

9,182,959 rows — 3rd node; 

9,178,352 rows — 4th node; 

9,178,352 rows — 5th node; 

9,186,544 rows — 6th node. 

Conspicuously, the table was distributed over the entire cluster. According to the configuration shown in Figure 1, 
the replication factor was 2, which means that each data block would be presented on 2 nodes. This can be seen from 
the number of rows on the nodes: the sixth node stored a copy of the first, the third — a copy of the second, the fifth — 
a copy of the fourth. 

The fault tolerance of this configuration can be checked by alternately disabling nodes in the cluster. To do this, you 
can turn off the node or stop services on one of the nodes with the systemctl stop clickhouse-server command. During 
the experiment, DBMS services were stopped on cluster nodes. 

With simultaneous disconnection of the 3rd, 4th, 6th or Ist, 2nd, 5th nodes that contained replicas, users continued 
to receive data from the dwh.cluster_test_data_distributed table, and the number of rows was equal to 27,547,855. 
When one of the nodes was disabled, the data continued to be displayed, and the number of rows was equal to 
27,547,855. When the nodes containing the replica and the original data were disconnected at the same time, data loss 
occurred. This configuration can be scaled to 12 nodes, then the replication coefficient will be 3, and the sharding 
coefficient will be 6. 

Discussion and Conclusions. The proposed solution can increase the speed of execution of analytical queries, solve 
problems with the scalability and reliability of data storage in banking organizations. The authors have automated 


cluster deployment by using templates in Ansible Playbooks, which provides setting up a ready-made cluster on new 
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servers in minutes. This configuration can be scaled by increasing the number of nodes and adding them to the 
configuration files. 

A set of characteristics that DBMS OLAP should correspond to was indicated, DBMS comparison was performed, a 
fault-tolerant cluster configuration of a data warehouse for analytical queries in the banking sector was proposed, 
automation of the configuration deployment process was performed. A similar solution is applicable for deployment on 
FreeBSD, Linux, macOS. The cluster configuration diagrams are given. This configuration can solve the problem of 


reliability and scalability, which is often found in organizations. 
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Abstract 
Introduction. At present, the concepts of fog and edge computing are used in a wide range of applications of various 


kinds. One of the key problems in the organization of computing in groups of mobile devices that make up the edge/fog 
layer is the mission assurance based on battery power availability. In this context, a lot of developments aimed at energy 
saving of device systems have been presented to date. However, one important aspect remains beyond the consideration 
of the problem of resource saving, namely, the issue of saving the residual resource of a computing device. The aim of 
this research is to formalize the workload distribution problem as two-criteria optimization problem, and to develop the 
basic solution technique. 

Materials and Methods. Within the framework of this article, an approach to resource saving is proposed. It is based on 
the evaluation of two device criteria: battery life and residual resource of a computing device. The residual resource of a 
computing device can be estimated using the probability of failure-free operation of the device, or as the reciprocal of 
the failure rate, taking into account that the exponential law of failure distribution is used in the simulation. From this, a 
model of the problem of two-criteria optimization is formulated, taking into account the dynamics of the network 
topology in the process of performing a user mission. The topology dynamics is reflected in the model as a sequence of 
topologies, each of which corresponds to a certain period of time of the system operation. 

Results. Based on the proposed model of the two-criteria optimization problem, a method was proposed for resource 
saving in the edge and foggy layers of the network. It reflected the specifics of the dynamic layers of the network, and 
also took into account the importance of the criteria for estimating the consumption of device resources. An experiment 
was conducted to evaluate the impact of the method of distributing tasks over a network cluster on the probability of 
failure-free operation of devices and on the average residual resource. 

Discussion and Conclusions. The conducted experiment has demonstrated the feasibility of using the developed 
method, since the distribution of tasks among executing devices had a significant impact (up to 25 % according to the 


results of the experimen t) on the average residual resource of a computing device. 
Keywords: resource saving, calculation planning, fog computing, edge computing, optimization. 
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AnHOoTanHna 

Beedenue. B nactosjee BpeMA KOHIeMIMH TYMaHHbIX HM KpaeBbIX BBIYHCIICHHM UCHON3YIOTCA WIMPOKHM KpyroM 
TIpHIOKeHHH CcaMOM pa3IH4HOM HanpaBNeHHOCTH. OHO U3 KIIOYeBLIX MpoOseM opraHH3allHH BbIYMCIeHHM B 
rpyumax MOOWJIBHBIX YCTpOMCTB, COCTaBJIAIONIMX KpaeBoli/TyMaHHbIi CIO, ABIAeTCA OOecHeyeHHe BbITIOHeHHA 
MHCCHH Ha OCHOBe HasIM4unA 3apaya Oarapen. B cBa3H Cc 3THM K HaCTOALIeMy BPeMeHH MpecTaBIeHO HeMaJIO 
pa3pa0oToK, HalipaBJIeHHbIX Ha 9HeprocOepexeHve CHCTeM ycTpolicTB. OHaKO OYeCHb BaxKHBIM acCIeKT OCTaeTcA 3a 
paMKaMu paccmMoTpeHHaA TpoOeMbI pecypcocbepexeHHa, a MMeCHHO — BOTIpoc cOepexKeHHA OCTaTOUHOrO pecypca 
BBIYMCIMTeIbHOTO ycTpolicTBa. Llenb1o WaHHoro McceqOBaHuA ABIAeTCA COPMaIM3alMA 3aqayu pacripezenenua 
Harpy3Ku Kak JBYXKpHTepHasIbHOH 3aa4n ONTHMH3alMH MU BbIOOp Oa30BOro MeTOa ee PeleHHA. 

Mamepuanei u memoooi. B pamkax JaHHO cTaTbu UpefyaraeTca MOAXoA K pecypcocbepexeHMIo Ha OCHOBe 
OWCHHBaHHA ByX KpuTepveB ycTpolicTB: pecypca OaTapeH M OCTaTOUHOTO pecypca BBIYMCMTebHOTO YCTPOMcTBa. 
OctaTouHbIii pecypc BBINHCIIMTeIbHOrO YCTpOlMCTBa MO2KeT ObITb OL|CHeH MPH MOMOLIM 3HaYeHHi BeposATHOCTU 
Oe30TKa3HOM paOoTEI ycTpolicTBa WIM Kak BeIMYMHa, OOpaTHad HHTCHCHBHOCTH OTKa30B C y4eTOM TOTO, 4TO TIpH 
MOJeJIMpOBaHHH CTIOb3yeTCA 9KCHOHCHIMAIbHbIM 3aKOH paciipefzemeHua oTka30B. Ha ocHoBe 9Toro 
CPopMysIMpoBaHa MOselIb 3aa4 JByXKpHTepHasIbHOM ONTHMU3allH C YHeTOM JHHAMUKH TOMOJOFHH CeTH B IIpolecce 
BBINOJIHCHHA TOb30BaTeIbCKOM MUccHH. J[MHaMuKa TOMONOrHH OTpaxkKeHa B MOJeJIM KaK MOCIIeOBATeJILHOCTb 
TOMOOIMH, KaxKad U3 KOTOPbIX COOTBETCTBYeT OlIpeeJICHHOMy OTpe3Ky BpeMeHH (PYHKUMOHMpOBaHHA CHCTeEMBI. 
Pezyismamoi ucciedoeanua. Ha ochopaHuu mpescTaBseHHOH MOJeM 3aa4H AByXKpHTepHabHOM ONTHMN3alHH 
Tipeso%KeH MeTOA ObecrieyeHHA pecypcocbepexeHHA B KpaeBOM HM TYMAHHOM CJI0e CeTH, OTpaKaloWIMi cneunpuKy 
J[MHAMMUYeCKHX CJIOeB CeTH, a TaKXKe YAMTHIBAIOMMM BaKHOCTh KPHTepHeB OMCHHBAHHA pacxoa pecypcoB ycTporicts. 
IIpopeyeH 9KCHepHMeHT, NOSBOAIOMIMH OLCHHTb BIIMAHMEe ciocoba pacripeyeueHuA 3aay MO CeTeBOMy KJIacTepy Ha 
BCPOATHOCTL Oe30TKa3HONM paooTs! YCTPOMCTB H Ha CpeqHHM OCTATOYHBIM pecypc. 

O6cyscdenue u 3akiio4uenua. VI[poBeqeHHbIii 93KCIIepHMeHT eMOHCTpupyeT WeecooOpa3sHocTb MmpHMeHeHHA 
pa3paOoTaHHoro MeTOJ{a, HOCKOIbKy paciipeyeueHve 3aja4 10 HCTOJHAIONIMM YCTpOMCTBaM OKa3bIBaeT CYLeCTBeHHOe 


BJIMAHHC (a0 25 % 110 UTOraM 9KCIIepHMeHTAa) Ha cpeaHnii OCTaTOUHbIN pecypc BbIYHCIINTeJIBHOTO yCTpovictBa. 


Ku1roueBble C10Ba: pecypcocbepexeHue, TIJaHHpOBaHHe BbIUMCJICHHH, TYMaHHble BEIYHCJICHHA, KLPacBbile BhIYHCJICHHA, 


ONTHUMU3alnaA. 


Baarogapuoctu. ABtop BbipaxkaeT IpH3HaTeJIbHOCTb pyKoBoscTBy MucTuTyTa MHHPOpMalMOHHBIX HaykK HM TeXHOJIOrHi 
Oe30nacHoctTuH Poccuiickoro rocyqapcTBeHHOTO TyMaHUTapHOrO YHUBepcuTeTa 3a NOMOLIb, OKa3aHHYy!O B Tpowecce 


TIOMrOTOBKH IIpOeKTa. 


Aisa waTupopanua. Knumeuxo A.b. J[ByxkputepuasbHbIi MeTo obectieueHua pecypcocbepexeHua B KpaeBoM U 
TyYMaHHOM CIIOAX ceTH. Advanced Engineering —_ Research (Rostov-on-Don). 2023;23(1):85—94. 
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Introduction. Currently, applications using edge and fog network segments are widely used [1-5]. Distributed 


computing performed on the nodes of these segments has significant difference from distributed computing performed 
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in cloud structures, namely, relatively high topology dynamics and, in addition, the presence of certain distances 
between the nodes used. This results in the need to take this parameter into account when modeling the time of 
information exchanges between executable tasks. The existence of topology dynamics is explained by the fact that 
network edge devices are, as a rule, user devices (including all kinds of sensors, smartphones, laptops, etc.) that belong 
to certain people and can be moved, turned off, access to data/resources may be prohibited. At the same time, all of 
them depend on the presence of a battery charge (power source) [6, 7]. The fog layer of the network is less dynamic. 
However, it is also characterized by all of the above to one degree or another: the occurrence of a certain distance 
between nodes with the presence of transit sections of the network, the possible mobility of the nodes themselves, e.g., 
in case of data processing by groups of mobile devices. An interesting example in this sense are groups of low-orbit 
satellites that are currently used for data processing, but the nodes are displaced relative to the Earth's surface at a fairly 
high speed (access time is on the order of several minutes) and, therefore, the question arises about the routing and 
transmission of processing results to ground stations [8]. 

As part of the issues of resource saving of devices used in the dynamic layers of the network, a number of works 
have been published [9-11], in which attention is focused on energy saving. As a rule, the model assumes the 
dependence of the consumed energy of the device on its workload [12]. At the same time, such an important factor as 
saving the residual computing resource of the device remained outside the focus of research [13]. 

The residual computing resource of a device is a value that is closely related to such reliability characteristics as the 
probability of failure-free operation and gamma-percentage time to failure. Also, the average residual life, taking into 


account the use of the exponential law of failure rate distribution, is characterized by the following expression: 
R(t) =f P(xydv a) 
P(t). 


where P(x) — probability of faultless operation (PFO) of the object during time 4 
For the exponential resource distribution law, the PFO is defined as: 
P(x)=e™, (2) 


where A — hazard rate, A>O. 

Then, according to formula (1), the average residual life is determined from the following formula: 

RQ)=R=~. 3) 
Xr 

At the same time, there are works devoted to the relationship between the values of PFO, gamma-percentage time to 
failure, and the average residual resource on the temperature of the computing element of the device, which, in turn, is 
described as a function of the device load [14—16]. The problem of device resource saving can also be posed as the 
problem of maximizing the average residual resource of the device. 

Thus, the problem of resource saving is considered in recent works either as a task of minimizing the power 
consumption of devices, or as a task of saving the residual computing resource. 

However, under present-day conditions, it is reasonable to talk about two-criteria optimization of resource 
consumption of devices for the following reasons: 

— electricity is a renewable but critical resource for the mission of the device, whether it is a sensor of an 
information and control system or a controlled mobile device; 

— residual computing resource determines the duration of the device's expedient operation, and the operation time, 


respectively, is the sum of the missions performed by the device. 
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Therefore, the objective of this study was to develop a resource saving method for devices of dynamic network 
layers, within which the assessment would be carried out on the basis of two criteria —- power consumption and 
computing resource consumption of the device. 

Materials and Methods. Let us define the concept of “mission of a group of devices” as a predetermined sequence 
of solving by a group of devices a complex of computational tasks related to information exchanges. The mission of a 
group of devices can be described by an acyclic graph, whose vertices will be weighted by the complexity of the tasks, 
and the arcs, respectively, will determine the limits of the sequence and the amount of data transmitted by the tasks. Let 
the mission be described by graph W = {w, z, I }, where w, — complexity of the task, z; — proportion of the task 
completed by the time of the topology change, J — matrix of data volumes transferred between tasks. A group of 
devices in a dynamic network layer is described by a set of graphs, each of which defines the topology of the network at 
time ¢;. We assume that the mission execution time includes a certain time interval with discrete instants of time /to, t), 
...f,], each of which corresponds to the topology graph G;=/{<p,j, ej, Ry>, H}, j=1...m, pi — performance of the j-th 
node of the i-th topology, e;; — energy resource of the j-th node of the i-th topology, Rj; — average residual computing 
resource of the j-th node of the i-th topology, m— number of devices in a group, H — matrix of network connections 
between nodes. We also assume that at each subsequent time ¢;,;, for graphs G;,; and G;, the following is true: at least 
one of the vertices belonging to G;, will be contained in G;,), which is required to continue the mission. 

We also assume that if the topology graph changes at time ¢;,;, it is possible to reassign tasks that were performed in 
the i- th topology on the nodes that ceased to exist at i+/ time, to any of the available topology nodes i+/. At the same 
time, we make the assumption that the tasks transferred to new nodes start to be executed at the same place where their 
execution was interrupted. 

Estimation of the residual computational resource will be implemented as follows: to estimate the PFO of the device 


of the fog layer, we will use the expressions proposed in [17, 18]: 


A 
Kay 2%, (4) 
where D — calculator load (in fractions). 
D= eo ; (5) 
P; *Econstraint 


where W, — complexity of the work performed. 


Accordingly, PFO (t) can be estimated, as well as the average residual resource: 


Apt: gkw/10 ptconstraint 


P(t) =e™ =e? =e ; (6) 


where teonstrainn — time limit for which the task must be completed by the node to meet the mission time limit f,,. 
When estimating the residual energy resource, we also assume that the energy consumption is proportional to the 


workload of the computing node, and there are no more energy costs: 
sore = En =e) ? (7) 


where EF jo — initial level of energy available, E jm — ftesidual level, & — coefficient expressing the relationship 


between the load of the computing element and its energy consumption. 
The result of the operation of the two-criteria method for ensuring resource saving will be a sequence of 
distributions of still unsolved mission tasks over time /fo, t;.;/ by topology G; taking into account the fact that some of 


the tasks that are in the process of being solved on the nodes inherited from topology G;.;, remain at their places and are 


not transferred: A = {A, ie 
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Aj=|.. ty -+/- (8) 


Matrix A describes the distribution of unsolved problems for which z,;</. At the same time, the number of nodes 
corresponds to new topology G;. Nodes inherited from G;.;, are renumbered and occupy indices starting from 1. 

Since the mission may involve more than one topology, we will decompose the problem as follows: we will use the 
“greedy” policy and assume that a satisfying solution (minimization of resource consumption, estimated by two chosen 
criteria) can be obtained, if for each topology, we choose the best solution. 


That is, the task model will have the following form: Vg,,z,<1t,>0 to get such Aj, that R(t,) 
R,(t,,) > max, e(t,,) > max , mpu max(A,,) <1, . 


Consider the possible priority of the formulated objective functions. For each individual topology, it is required to 
minimize both the power consumption and the consumption of computing resources. In some special cases, the 
coincidence of optimization goals for two criteria can be achieved. Specifically, the task is distributed to the node in 
such a way that it will provide it with the minimum workload of the computational element and at the same time the 
minimum power consumption. 

However, this situation is not always possible. For example, there is a node loaded less than other calculations, but 
at the same time its energy resource is almost exhausted. By assigning a task to it, we jeopardize the mission (or it is 
necessary to change topologies as the supply of electricity runs out), and by assigning a task to a more loaded node, we 
worsen the value of the residual computing resource criterion. The opposite situation can also occur: a node with a large 
residual power resource may have a low residual computing resource as a result of participating in previous device 
missions. 

Therefore, it is advisable to rank the objective function (OF) in terms of importance, on the basis of which to select a 
further solution method. 

For critical missions, the battery life of the device plays a leading role [19-21]. Accordingly, the OF that determines 
the residual computing resource can be reduced to an additional constraint. Uncompleted mission tasks will be 
redistributed depending on the remaining battery charge. 

It should be noted that the above is also true for missions that do not have critical execution. In the event that the 
device loses battery power without completing the mission, it may be lost. Then the optimization of the computing 
resource will become meaningless. 

On the other hand, by reducing the residual computational resource to a limitation, the solution may not be obtained 
at all due to the lack of available nodes whose characteristics would satisfy it. 

Research Results. Multicriteria optimization problems, as a rule, are reduced to solving one or more single-criteria 
problems. In this case, it is more expedient to leave one OF, and present the second one as a constraint. 

For the system under consideration, the fulfillment of the mission is critical for each node; therefore, we transform 
the OF for energy consumption into a restriction, the implementation of which will cut off unacceptable options for 
distributing tasks among the nodes. Further, within the framework of the resulting set of nodes, the distribution will 
occur with minimization of the consumption of the residual computing resource. 

Let us describe the main stages of the two-criteria method for ensuring resource saving as the basic result of the 
study: 

—on the existing G,-topology, select nodes that have sufficient energy resources to complete the mission without 


changing the topology; 
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— distribute mission tasks among selected nodes in such a way as to minimize the consumption of computing 
resources for idle nodes, to which unsolved tasks will be attached. 

It is possible to rank the nodes obtained at the first step of the method in descending order of the availability of an 
energy resource. But then, by distributing tasks primarily among nodes with the maximum amount of energy resource, 
there may be a deterioration in the values of the residual computing resource criterion. 

Let us illustrate the dependence of the PFO values of devices on the selection of a node for solving the problem (in 
this experiment, we compare the states of the nodes under the transmission of incoming data and during their 
processing). 


The topology is shown in Figure 1. 


Fog node 
100 


Edge device 
100 Fog Layer Cloud 
Broker 1,000 
500 
Fog node 
300 


Fig. 1. Experimental fragment of network topology 


The experiment was carried out for 10-time intervals, each of which was 100 hours with the following parameter 
values: Wreceive (amount of data received by the task) — 500 units, Wsend (amount of data sent by the task) — 
100 units, Wprocess (data processing complexity) — 150 units, Tdecl (time allotted for solving the task) — 50 units. 

The following are graphs of PFO nodes depending on whether the node processes data or only transmits it (Fig. 2—5). 
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Fig. 2. Comparison of PFO values for node 1 under data processing (Pp1) 
and when transferring data (Ptr1) 
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Further, Figure 6 shows comparison of differences in the values of the average residual computing resource. 
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3. Comparison of PFO values for node 2 under data processing (Pp2) and when transferring data (Ptr2) 
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. 4. Comparison of PFO values for node 3 under data processing (Pp3) and when transferring data (Ptr3) 
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5. Comparison of PFO values for node 4 under data processing (Pp4) and when transferring data (Ptr4) 
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Fig. 6. Comparison of average residual resource values for nodes 1—4 during data transit 
(Remain_transit) and under data processing (Remain_perform) 


Discussion and Conclusions. The article proposes a two-criteria method for ensuring resource saving in the edge 
and fog layers of the network. It differs from those previously presented in the subject area of organizing distributed 
computing in the dynamic layers of the network in the following. The proposed method takes into account both the state 
of the battery of the device (which is critical for the mission) and the residual resource of the computing device (what 
affects the duration of the appropriate use of the device outside the mission). 

The developed method is based on the formulated problem of two-criteria optimization, which is then proposed to 
be reduced to a single-criteria one through converting the objective function of the energy resource balance into a 
constraint. 

The nodes that are not suitable for the battery level are cut off in advance. Then there is a distribution of tasks 
according to the criterion of the residual computing resource. 

The conducted experiment makes it possible to compare the results of task distribution over nodes in the presence of 
information interactions between tasks. It is clearly shown that the method of task distribution affects significantly the 
value of the residual resource of a computing device (up to 25 %). Thus, it is possible to improve the resource 


consumption indicators of devices involved in calculations in the edge and fog layers of the network. 
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Abstract 

Introduction. Essential nonlinearity of the chemical reactions of acids and bases determines the control algorithms in 
the mode of acidification or alkalization, that is, periodic dosing of a minimum volume of acid or alkali. Such regulation 
may be ineffective, specifically, it allows insufficient or excessive concentration of the controlled substance. The article 
discusses the problem of precise regulation of the hydrogen index in mini-bioreactors. It is proposed to use a digital 
model of the acidity control system to select the concentrations of topped-up solutions, determine the regulation 
methodology, and improve accuracy. The objective of the work is the assurance of required accuracy of pH regulation 
in an in vitro mini-model of the gastrointestinal tract of a static type. 

Materials and Methods. The initial block diagram of the model included accumulators and flows. It was the base for 
the main differential equations characterizing the change in volume and acidity. To correct the acidity readings of the 
resulting solution by temperature, a static model based on the polynomial approximation of experimental data using the 
least squares method was created. The structural elements of the mathematical model were investigated in the Matlab 
Simulink application package. To validate the adequacy of the mathematical model, transient characteristics were 
determined on a real system of in vitro modeling of the artificial gastrointestinal tract of poultry. 

Results. Within the framework of this work, the authors created and analyzed a nonlinear mathematical model of pH 
changes in a bioreactor taking into account external control actions. The flows of hydrochloric acid solution, alkali 
solution and drain from the reactor were presented as elements of a differential equation describing the accumulation of 
liquid in the reactor. To improve the accuracy, the solution was modified taking into account the temperature 
dependence of the hydrogen index. A dosing mathematical model based on a regulator with alkali and acid channels 
was proposed. The data obtained made it possible to generate a combined model of the pH regulation process in the 
bioreactor. The adequacy of the solution was confirmed empirically. The models of pH regulator, regulation of the 
volume of contents in the reactor and chemical reactions were shown in the form of structural diagrams. The transients 
of a mathematical model and a real control system were compared. It was established that the transient characteristics of 
the mathematical model and the real system were identical in terms of regulation time. The relative error of regulation 
of the real system was 0.35 %, and the mathematical model — 0.1 %, which corresponded to the required accuracy of 
regulation + 0.1 pH. The influence of the studied flows on the neutralization reaction was shown in the form of graphs. 
Discussion and Conclusions. The proposed mathematical model will provide selecting optimal methods and algorithms 
for regulating acidity, which will accelerate the creation of a regulator for the nonlinear process of regulating the 


hydrogen index. In the future, these developments can be integrated into a comprehensive digital model of the entire 
artificial gastrointestinal tract of poultry to optimize control algorithms (dosing, mixing, periodicity, etc.), as well as 
approximation to objects in vivo. 


Keywords: mathematical modeling, acidity, pH, in vitro modeling, in vivo modeling, control system, acidity regulation 
algorithm, digital model. 
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AHHOTalna 

Beedenue. CyjecTBenHad HesIMHeMHOCTh XMMMYeCKHX peakUMii KMCIIOT HM OCHOBaHHM oMpeyemaeT alropHTMbl 
yiipaBsIeHHA B peyKHMe NOAKMCICHHA WIM MOAMWeIaMBaHHA, TO CCTb MepHOAMYeCKON O3alHH MHHMMasIbHOrO OOBeMa 
KMCJIOTbI HIM Weroun. Takoe perymupoBaHuve MOxeT OBITh Maslo9dPeKTHBHBIM, T. €. JOMycKaeT HeOCTaTOYHy!IO WH 
W30bITOUHYIO KOHUeCHTpalluio KOHTposMpyemoro BelljecTBa. CTaTbA MocBAuleHa MpoOsemMe TOYHOTO perysIMpoBaHHA 
BOJOPOAHOTO MoKa3aTeIa B MMHH-OuopeakTopax. IIpeqnaraetcd MCHONL30BaTb WMPpOByIO MOJ[eJIb CHCTeMBI 
ylpaBsIeHHa KHCIOTHOCTbIO JIA WOAOopa KOHUeHTpalMii MONMBaeMbIX pacTBOpOB, OlIpeseyeHHA MeTOAMKH 
peryIMpoBaHHA MW MOBbILeHHA TOUHOCTH. Liens padoTEr — obecteyeHHe TpeOyemMol TOuHOCTH perymupoBaHua pH 
B in vitro MMHH-MOJeJIM 2KeJTy [OWHO-KHMIe4YHOrO TpakTa CTaTHYeCKOrO TUM. 

Mamepuanoi u memoooi. Ucxoyuas crpykTypHad CXeMa MOJeJIM BKIKOUAeT HAKONMTeIM MW MOTOKH. Ona pescTaBiaerT 
coOoi Oa3y AIA OCHOBHBIX WudpdepeHuUMaIbHbIX ypaBHeEHHH, XapakTepv3yIOWIMX H3MeHeHHe OObBeMa HM KHCIIOTHOCTH. 
Jia KOppeKTHPOBKH TOKa3aHHii KMCMOTHOCTH pe3yIbTHpyrOllero pacTBopa MO TeMIlepaType co3qaHa cTaTHuecKkaad 
MOJesIb, OCHOBaHHad Ha MOJMHOMMAIbHOM allMpOKCHMalM 9KCIePHMeHTAJIbHbIX JJaHHbIX Me€TOJOM HaMMeHBbIWIMX 
KBagpatoB. B mpuknaqHom makete MATLAB Simulink uccneqoBaHbl CTpyKTypHbIe 3JIEMeCHTbI MaTeMaTHYeCKoH 
Moyes. Ha peambHoi cuctTeme in vitro MojjeuMpoBaHuaA UCKYCCTBCHHOMO 2KeslyJOUHO-KMINedHOTO TpakTa JOMalHel 
IITHUBI Ope eeHbI MepexOMHble XapaKTePHCTHKH JIA NOATBepxK CHUA aeKBAaTHOCTH MaTeMaTHYeCKOH MOEN. 
Pezsyismamot ucciedosanua. B pamMKkax jaHHo paOoTbI aBTOpbI co3qaiu MW MpoanaM3HpoBalM HesIMHelHyt10 
MaTeMaTHYeCKy!O MOJeIb H3MeHeHHA PH B OMOpeakTope c yYeTOM BHEINUHHX yupaBIAIOWMx BosfelicTBui. [loroKu 
pacTBopa COJAHO KUHCIOTHI, pacTBOpa WesOUH U CIMBa pacTBOpa M3 peakTOpa Ipe{CTaBJIeHbI KaK 93JIEMCHTHI 
TupPepeHunalbHOro ypaBHeHHA, OMMCbIBAaIOWerO HaKOMIeHHe %KHKOCTH B peakTope. JIIa NOBbILIeHHA TOUHOCTU 
pewieHHe opadoTamm c y4eTOM TeMMepaTypHOM 3aBHCHMOCTH BOZOpogzHoro noKa3atena. Ipeqnoxena 
MaTeMaTHUYeCKad MOJeJIb NOSHPOBAHHA Ha OCHOBE perysaATopa C KaHalaMM Weoun WU KucMOTHI. lomyaeHuble JaHHbIe 
TO3BOIMIM TeHepHpoBaTb OObeAHHEHHYIO MOJeIb Mpolecca perymupoBanua pH B Onopeaxtope. AyeKBaTHOCTB 
pelleHHua NOATBepAWIM OMBbITHDIM TyTeM. B Bue CTPyKTypHbIX CXeM [MOKa3aHbI MOJeuIM: perysatopa pH, 
perysIMpoBaHua oObeMa cofepxKHMOrO B peakTOpe H XHMHYeCKHX peakuMit. CpaBHuBaloTcaA MepexOHble IMpOLlecchl 
MaTeMaTHYecKOH MOJeM UW pealibHOM CHCTeMbI yiipaBleHus. YcTaHOBJIeHO, YTO TepexoHble XapakTepHCTUKU 
MaTeMaTHYuecKOH MOeH HU peasbHOM CHCTeMbI UJCHTHYHbI TO BpeMeHH perysMpoBaHua. OTHOCHTeIbHaAd 
MOrpeLIHOCTb perysIMpOBaHHA peasbHOM cucTembI cocTaBusia 0,35 %, a MaTemaTu4eckoH Moyenu — 0,1 %, 4To 
COOTBeETCTByeT TpebyemMoli TOUHOCTH perymupoBaHua +0,1 pH. B Buge rpapukosB oKa3aHO BIMAHMe MCCIeOBAaHHBIX 


MIOTOKOB Ha peakIM1O HeliTpalu3allun. 

O@6cystcoenue u 3aknrouenua. IIpeqnaraemad MaTeMaTHYeCcKad MOJesIb MO3BONMT MOLOOpaTh ONTHMAJIbHbIC MCTO]bI 
asIFOPHTMB! peryIMPOBaHHA KUCIIOTHOCTH, YTO YCKOPUT Co3aHUe peryiaATOpa HeMHeMHOTO Mpowecca peryiMpoBaHHuAa 
BOOpozHoro HoKa3atesa. B OyayWieM 3TH HapaOOTKH MO2%KHO HHTerpHpoBaTb B KOMIMICKCHY!O IMpoBy!lo MOJeJIb 
BCerO HCKYCCTBCHHOTO Key IOUHO-KHMICYHOTO TpakTa JOMAIIHeH UTHIbI Td ONTHMU3aIIMM aNrOpHTMOB yipaBseHHa 
(03allMu, NepeMelIMBaHHA, TepHOAMUYHOCTH U T. I1.), a TAKKe AallIIpOKCHMAallHH K OOBEKTaM in Vivo. 


Koro4ueBble CsI0Ba: MaTeMaTHYecKkoe MOJCJIMNPOBaHHe, KHCJIOTHOCTH, pH, MOJCJINPOBaHHe in vitro, MOJCJIMpOBaHHe in 
vivo, CHCTeMa yUpaBJIeHHsA, asITOPpHTM peryJINpOBaHHA KHCJIOTHOCTH, mudposasa MOJICJIb. 
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Baarogapuoctu. ABTopb! BbIpaxkaloT OarofapHOcTb 3a (PHHAaHCOBY!IO NoOfWepxXKy UCcieqOBaHHit B paMKax rpaHTa 
MuuuctepctBa HaykKH HM BBICHIero oOOpa30BaHHa Poccniicxoi Wegepaunun Ne 075-15-2022-285 or 9 urona 
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Introduction. Since the early 1980s, the scientific community studying the physiology and functional features of 
human internal organs has come to understand that in the near future, the focus of experimental techniques will shift 
from “in vivo” practices to “in vitro” systems [1]. The development of systems for hardware simulation of physiological 
processes will be of primary relevance, including for the following reasons: reduction of costs and time for conducting 
research; possibility of simultaneous conducting of the same type of experiments; and, what is not unimportant — 
increasing the ethics of research, reducing the need for experiments on living beings. By the beginning of the XXI 
century, a whole scientific direction was formed, engaged in the development of techniques and equipment for “in 
vitro” modeling and research [1]. 

Currently, in vitro models of the gastrointestinal tract are divided into two classes: static and dynamic. Static 
installations allow you to simulate the processes of digestion when using liquid nutrient media [2]. Dynamic models, 
unlike static ones, allow modeling peristaltic movements in the gastrointestinal tract, and are used when working with 
content nutrient media. 

From the point of view of the simulated volume, the gastrointestinal tract models (regardless of the previous 
classification) are divided into micro- (up to 50 ml), mini- (from 50 to 400 ml) and macro-systems (over 400 ml) [3]. 
Microfluidic systems should be mentioned separately [4], but they go far beyond the scope of this work. 

Patent analysis makes it possible to determine that the first is a patent for an invention “In vitro model of an in vivo 
digestive tract”'. This system is dynamic and precedes the market entry of the first TIM-2-type systems [5]. The 
development of this approach occurred in the models DGM [6], TIM-2 [5], HGS [7], etc. Static models, such as: DIDGI 
[8], SIMGI [9], SHIME [10], ARCOL [11], etc., moved away from complex modeling of peristalsis, and focused on 
more accurate modeling of the actual digestive processes. One of the reasons for this is the technical complexity and 
high cost of implementing systems with dynamic content compression. Another reason was the development of 
mathematical modeling methods and microcontroller control systems, which made it possible to control the medium in 
a test tube (reactor) with high accuracy and stability [4]. 

In this paper, we consider the solution to a particular problem of mathematical modeling of the pH regulation 
process for an in vitro mini-model of the gastrointestinal tract (gastrointestinal tract) of a static type [12]. The 
importance of precise control of the acidity of the chyme is determined by rather strict requirements for the conditions 


of in vitro experiments. Ensuring plausible conditions for the course of digestive processes, as well as the requirements 


for repeatability of experiments, require automatic pH control in the bioreactor with an accuracy not worse than +0.1 
units.” At the same time, the pH range can be from slightly alkaline (7.9 pH, oral cavity) up to strongly acidic (1.3 pH, 
stomach). 

The objective of this work is to ensure the required accuracy of pH adjustment in an in vitro mini-model of the 
gastrointestinal tract (gastrointestinal tract) of a static type. 

To achieve this goal, the following tasks will be solved in the work: 

— development of mathematical model of nonlinear pH change in a bioreactor taking into account external control 
actions; 


— study of a mathematical model by simulation methods; 


'US5525305A. In vitro model of an in vivo digestive tract. URL: https://patents.z0ogle.com/patent/US5525305A/en 
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“No. 075-15-2019-1880. Veterinarnye probioticheskie preparaty napraviennogo modulirovaniya zdorov'ya zhivotnykh. URL: 
https://probioticdonstu.com/itno-2021-ru/ (In Russ.) 
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— refinement of the mathematical model in order to increase its accuracy by taking into account the temperature 
dependence of the hydrogen index; 

— development of a mathematical model of the dosing process based on a two-channel competitive regulator (via 
alkali and acid channels); 

— development of a unified mathematical model of the pH regulation process in a bioreactor; 

— conducting experimental studies to confirm the adequacy of the combined mathematical model and determine the 


achieved accuracy of regulation. 


Materials and Methods. To simulate the neutralization reaction, we will use the method of accumulators and 
flows [13]. The method is based on the use of fundamental or local conservation laws (accumulators) and mathematical 
description of patterns that change the amount of accumulated feature (flows). 

The block diagram of the model is described in the figure, and contains two drives: 


(1) Volume 
acid 
atts (3) Drain 
(3) Volume 
alkali 
Substance 
(4) Flow 
siciid Q> (6) Digestive 
(5) Flow pubstaliee processes 
HC1 


Fig. 1. Block diagram of the model 
The reactor as a physical object simultaneously acts as a storage of the amount of liquid (Q1) and a storage of the 
amount of hydrogen ions (Q2). This is a container in which a controlled digestion process takes place, accompanied by 
a change in pH and a change in fluid volume. 
There are three threads acting on the first drive: 
— the first stream is the volume of acid topped up; 
— the second stream is the volume of added alkali; 
— the third stream is the volume poured out of the reactor. 
The second drive is also affected by three threads: 
— the fourth flow is the supply of sodium hydroxide by a pump [2]; 
— the fifth stream is the intake of hydrochloric acid. [2]; 
— the sixth stream is a potential change in pH during digestion. However, from a methodological point of view, it is 
a disturbing effect, and is not discussed in this paper. 
Now, let us consider the mathematical implementation of the model for describing physico-chemical processes. 
The volume of liquid in the reactor is described by the following equation: 
V=S-h. (1) 
Therefore, the change in the volume of liquid in the reactor: 
dV =S-dh. (2) 
For convenience, let us call the streams: 
— the hydrochloric acid solution flow: q1=G1; 
— the flow of the alkali solution (sodium hydroxide): q2= G2; 


— the discharge flow of the solution from the reactor: q2 = G3. 
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Thus, the differential equation describing the accumulation of liquid in the reactor: 
dh_1 
To determine the dynamics of changes in the concentration of the solution in the reactor, it is required to find the 
differential of the amount of substance according to the following formula: 
v=V-C, (4) 
where V is the volume of the substance, 1; C is the molar concentration, mol/l. 


From (4) and (1), we obtain: 


# -£(s-h-c)=s(cBrnZ), (5) 
dt dt dt dt 


The change in the concentration of the solution is influenced by the flows of alkali and acid with a certain 


concentration, then: 
dh dC 
s(c-Ben E) = (6,-Cyou -Gs-Gra) (6) 


From equality (6), we express the change in concentration and simplify the expression: 


dC 1 
dt = rare (Cynon =) 4G (C=Caa)) . (7) 

Thus, the system of differential equations describing the dynamics of changes in the volume and concentration of 
the solution in the bioreactor is presented in formula 8: 


dh_ 1 


a = 5G +G, —G,), 

(8) 
dC 1 
FG Cul Canon ~C)+ G:C Cc). 


To obtain the concentration of H ((H+]) cations from the solution of the system of equations (8), let us recall the 
basics of chemistry. The concentrations of hydrogen ions ([H+]) and hydroxide ions ([OH-]) in distilled water at 25 °C 
are equal and amount to 10-7 mol/l, this directly follows from the definition of the ionic product of water, which is 
presented below [11]: 

K = [H*] x [OH] =10'mol’/?? (at 25 °C). (9) 

It is also worth considering that: 

T =[OH ]-[H"]=[Na]{Cr]. (10) 

Now we find [H+] (denoted as Cy) through the square root of the previous expressions, taking into account the 


concentration of the solution according to equation (4) of our mathematical model [14]: 


C,, =VK,T =0, (1) 


Where ions i = 4. 
To create a control system with feedback on acidity in the form of a hydrogen index (pH), it is necessary to calculate 


a negative decimal logarithm from the concentration calculated according to equations (7). Thus, we get the acidity of 


the solution under normal conditions, that is, at 25 °C. 
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But the contents of the bioreactor are not always in normal conditions. For example, to create an artificial 
environment characteristic of the in vivo environment of the gastrointestinal tract of poultry, a temperature of 42°C is 
maintained in the reactor’. Therefore, experimentally and on the basis of available research* a matrix of changes in the 
pH of buffer solutions depending on temperature was obtained on high-precision pH meters, and the readings were 
obtained using a measuring system with thermal compensation for the electrode used? [13]. 


Based on these data, a static model was constructed, it was obtained by polynomial approximation by the least 
squares method. A visual representation of the pH dependence on temperature is shown in Figure 2: [14-16]. 


Real pH 
14 


12 


10 


2 
pH 25 °C 


Fig. 2. A surface describing the dependence of pH on temperature (12) 


F(x) = p00+ pl0-x+ p0O1- y+ p20-x* + pll-x- y+ p02: y? + 12) 
+p30-x° + p21-x?- y+ pl2-x-y’. 


The coefficients of the equation are presented in Table 1. 
Table | 


Coefficients of the equation (11) 


Coefficient Value 
poo 0.154 
plo -0.00511 
pol 0.899 
p20 -2.998 - 10° 
pli 0.00333 
po2 0.0130 
p30 -3.226 - 107 
p21 1.588 - 107 
pl2 -4.98. 107 


To determine the transients in the system, we will build a mathematical model in the MATLAB Simulink 
application package (Fig. 3): 


> Automatic Temperature Compensation in pH Measurement. _ https://www.horiba.com/esp/water-quality/support/technical-tips/bench- 
meters/automatic-temperature-compensation-in-ph-measurement/?utm_source=uhw&utm_medium=301&utm_campaign=uhw-redirect 

Toa 

Ibid. 
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Gnaiger E. pH measurement and temperature dependence of pH. URL: 
https://www. bioblast.at/images/archive/b/be/20190329092613!MiPNet08.16_pH-Calibration.pdf 
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V_ pH(25C) 


Static model pf pH 
adjustment from 
solution temperature 


Model of the dosing 
control system 


Volume 
change model 


Neutralization 
reaction model 


Fig. 3. General structure of the mathematical model in MATLAB Simulink 


Let us consider the created system by blocks. The first block is a block of a two-channel acidity regulator with 
feedback. Noise has been introduced into the feedback circuit to simulate random events of a real system for measuring 
acidity and measuring instrument errors. 


kite 


Switch H Discrete Controller H 


Regulation error 


Real pH 


Real pH GC) 


Fig. 4. Structure of the pH regulator, block “Model of the dosing control system” 


The second unit is a system for monitoring the volume of solution in the reactor. The “Drain volume controller” 
function simulates the activation of pumping. If the volume of the solution approaches the maximum volume of the 
reactor, the function gradually resets the solution to the required volume. 


x 1/S (reactor) 


Volume to h 


Fig. 5. Structure of the reactor content volume control model, block “Volume change model” 


Information Technology, Computer Science and Management 


101 


http://vestnik-donstu.ru 


102 


Advanced Engineering Research (Rostov-on-Don). 2023;23(1):95—106. eISSN 2687-1653 


The third block is the main one. It implements the calculation of the chemical reaction of acid and alkali when added 
to the reactor and calculates the pH of the final solution under normal conditions. The function called “Neutralization” 
calculates the concentration of hydrogen ions in a solution [12]. And the “Concentration in pH” block implements the 


conversion of concentration into a hydrogen pH indicator. 


pH (25 C) 


Concentration in pH 


-( Neutralization) 


Fig. 6. Structure of the model describing the neutralization reaction 


The fourth block enables to take into account the temperature changes of the solution with the function of the 


statistical model in the “Temperature correction” block. 


Temperature correction 


Fig. 7. Structure of the model describing the change in acidity from temperature, block “Static model of pH adjustment from solution 


temperature” 


To determine the adequacy of the transient characteristics of the resulting mathematical model, experiments were 
conducted on a simulation model. The real model is represented by an automated control system for the artificial 


gastrointestinal tract of poultry with one control unit and a bioreactor. Figure 8 shows the experimental equipment and 


its description. 
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a) b) 
Fig. 8. Research equipment and its description: 
a — bioreactor environment control units; b — elements of bioreactor. 1 — controlled agitators; 2 — heaters; 
3 — ports with sensors and a system for topping/draining solutions 


The methodology of the experiment is as follows. The mathematical model takes into account the characteristics of a 
real system, such as: 

— set acidity; 

— reactor parameters (base area and height of the bioreactor); 

— concentrations of added substances (acids and alkalis); 

— initial concentration of the solution; 

— initial volume of the solution; 

— the temperature of the solution throughout the experiment; 

— discreteness of the control process [14]. 

The mathematical model is an ideal system in which mixing of solutions in the reactor occurs instantly and 
uniformly. In real conditions, this is impossible, therefore, forced mixing is needed for the uniform dissolution of the 
added acid or alkali. In the simulation model, this problem is solved by a magnetic stirrer controlled via a SCADA 
system. 

Unfortunately, this type of mixing device induces electromagnetic interference on the pH electrode. To eliminate 
interference, the mixing device was switched off for a short time, and the discreteness of measurements and the 
operation of the regulator was 20 seconds, which is due to the experimentally established time of normalization of the 
pH readings of the electrode in the reactor (according to the documentation of the electrode, sampling can reach 1-2 
minutes).°’ The regulator of the real system is limited in the maximum volume of the dosed acid or alkali in order to 
ensure safety, which was taken into account in the model regulator. 

Experimental Results. Figure 10 shows a comparison of the transition process of a mathematical model and a real 
control system. The pH electrode is a sensitive element that is susceptible to receiving external electromagnetic 
interference. 

One of the experiments was carried out at a given acidity of the solution in the range of 2.15 + 0.05 pH. The acidity 
of the initial solution was 6.2 pH. The temperature of the solution in the reactor: 22.25 °C. Solutions of 0.1 mol/liter 
NaOH and 0.1 mol/liter HCl were used to change the acidity. The dosing algorithm consisted in discrete proportional 
regulation. 


° ph-4502c pH meter calibration notes. URL: https://tlfong01.blog/2019/04/26/ph-4502c-ph-meter-calibration-notes/ (accessed: 30.11.2022). 
7 EZO™ pH Circuit Datasheet. URL:https://files.atlas-scientific.com/pH_EZO_Datasheet.pdf (accessed: 30.11.2022). 
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Fig. 8. Transient characteristics of a mathematical model and a real system when regulating from 6.2 to 2.15 pH 


Figure 11 shows the flow integrals from differential equation (7). 


6-10° T T T T x T T T T , — GOH] 
~~» G[H] 
4 4-108 + —— 4 G, [fluid] 
bp an --- G, [fluid] 
E 
3 
a + 
amy 
4 4 = 


350 400 450 ‘Time, s 
Fig. 9. The effect of each flow on the neutralization reaction during experiments 
Discussion and Conclusions. Analysis of the simulation results have shown that flows with a solution 
concentration in the reactor do not make a significant contribution to the regulation of the entire system, therefore, they 
can be neglected for simplification. Then the system of differential equations of the digital model (8) will take the 


following form: 


h 1 
TE eG 26.26); 
d Ss : 
dC 1 aa 
rs = a *Cyaon ~ Gr *Cuer ) 
According to the results of the experiments, taking into account the concentration of the solution in the reactor, the 
relative error of the steady-state acidity of the solution for the mathematical model was 0.1 %, and for a real system 


with the same algorithm 0.35 %. 


Thus, taking into account the error of the electrode readings and the neglect of the mixing process of solutions in the 


mathematical model, the absolute error of the control results is within acceptable limits and does not exceed + 0.1 pH of 
the set value. This means that the developed mathematical model adequately describes the neutralization reaction curve, 
and in the future, it will allow choosing optimal control algorithms, which would be a laborious process on a real 


system. 
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The developed model solves the problem and is a mathematical tool for the development of methods and algorithms 
for pH regulation, taking into account the temperature of the contents of the mini-reactor, the initial concentration of the 
main solution and the concentrations of top-up solutions, which will improve existing control algorithms. 
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